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CVI. The Sensitivity of Electron Diffraction as a Means of 
Detecting Thin Surface Films : I 
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ABSTRACT 


Very thin layers (0-204) of silver bromide and of copper have been 
deposited on to smooth (111) surfaces of silver single crystals. The 
thickness of these layers has been determined by means of a radioactive 
technique. The electron diffraction patterns from the layers have been 
studied as a function of their thickness. The layers do not grow as 
uniform films but as nucleated deposits. It is found that the thinnest 
detectable layer is 0-44 in mean thickness for the silver bromide and 
0-8 A for the copper. 


§ 1. INTRODUCTION 
Evectron diffraction provides a powerful tool for the study of thin 
surface films when the reflection method is used. Although information 
concerning crystal size, structure, orientation and lattice spacing is 
readily obtainable, it is extremely difficult to estimate the average 
thickness of a very thin layer. In particular, the minimum thickness of 
a detectable surface film is not known. 

Various types of surface deposits, of 100A or less in mean thickness, 
are commonly studied. These layers are formed by electrodeposition, 
chemical attack, condensation from the vapour in vacuo and growth 
from solution. The thickness of the layers may be estimated by such 
methods as weighing, multiple beam interferometry and from Faraday’s 
law in the case of electrodeposits. For film thicknesses of less than 104, 
these methods become insensitive or unreliable. 

The direct estimation of film thickness from the electron diffraction 
pattern is limited by : 

(a) the lack of a detailed theory which predicts the intensities of the 
diffracted beams relative to the background intensity, 

(b) the marked influence on the pattern of the substrate surface 
topography, 

(c) the effect on the pattern of the mode of growth of the deposit. 


Raether (1949, 1950, 1951) has considered the effect of (b) on the 
assumption that the deposit grows as a uniformly coherent film, which 
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follows the topography of the substrate. He considered three special 
cases : 
(1) electron optically smooth surfaces (e.g. crystal cleavage surfaces), 
(2) slightly undulating surfaces (e.g. electropolished surfaces), 
(3) rough surfaces (e.g. etched or abraded surfaces). 


The sensitivity for observation of the thinnest deposits decreases as 
the surface. becomes progressively rougher. In the case of a smooth 
surface, the deposit will cause the diffraction pattern from the substrate 
to be obscured below a certain value of the angle between the electron 
beam and the surface (hereafter denoted by «). If this critical value of « 
can be determined, and if the relevant mean free path for electron 
scattering by the deposit material is known, then the thickness of the 
deposit may be estimated. Raether gives an approximate method of 
determination of this mean free path. The method has been used by 
both Raether (1950) and Lucas (1951) to estimate the thickness of zine 
oxide on a cleavage surface of zinc, but it is limited and not very 
accurate. 

Schulz (1951, 1952 a, b, c) has systematically studied the growth of 
alkali halides on various substrates, and he estimated the thickness of 
the layers, to an accuracy of -+-20°%, from the radiation characteristics 
of the source. He calibrated this by depositing a thick layer, the 
thickness of which was determined either by weighing or by an optical 
method. He claimed that with a cleavage face as substrate an observable 
diffraction pattern could be obtained from a deposit only one or two 
Angstroms in average thickness. Seifert (see Gomer and Smith 1953, 
pp. 360 and 383) has expressed doubts as to whether such thin layers 
can be detected in this way. 

An accurate and reliable method of thickness measurement is therefore 
required to test the validity of Schulz’s conclusions concerning the 
sensitivity of the electron diffraction method. Now that a large number 
of radioactive isotopes are available, it is possible to use tracer methods 
for thickness measurements. Deposits have therefore been formed on 
flat substrate surfaces and the change of the diffraction pattern has 
been studied as a function of the deposit thickness in the range 0-20, 
this thickness being measured by a radioactive method. In particular, 
the thinnest deposit layer which would give rise to an observable 
diffraction pattern has been determined. The substrate consisted of 
atomically smooth (111) faces of silver crystals, and the deposits consisted — 
of (i) silver bromide grown by chemical attack of the silver and (ii) metallic 
copper deposited from the vapour phase. 


§ 2. SPECIMEN PREPARATIONS 


The silver substrates were prepared by condensation on to a cleavage 
surface of mica heated to about 270°c. <A detailed account of the 
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preparation and properties of these specimens will be given later. The 
main advantages for the purpose of these experiments are : 


(1) the specimens are easily prepared in large numbers, 


(2) the silver is well oriented, in single positioning, with its (111) face 
parallel to the mica surface, 

(3) the surfaces are ‘ atomically smooth’ over areas large compared 
with atomic dimensions. 

The specimens were mounted on small brass blocks, which were 
clamped directly in the specimen holder of the camera. The edges of 
the mica were painted with ‘ Aquadag’ which prevented charging up 
of the specimens in the electron beam, and allowed fairly low values of « 
to be used. 


(a) Bromide Layers 

The bromine was prepared by the oxidation of ammonium bromide, 
which had been exposed to neutrons in the Harwell pile. The bromine 
gas was dissolved in water to give a stock solution, from which subsequent 
working solutions were obtained by dilution. The silver specimens were 
exposed to the bromine vapour, which was in equilibrium with a given 
strength of solution, by suspending them above the surface of the solution. 
By trial and error, a strength of solution was found which gave a convenient 
rate of growth of the bromide of about 3A/min. The exposure times were 
varied so that specimens of different thickness were obtained ; these 
specimens could be stored indefinitely without any loss in the quality of 
their diffraction patterns. 


(0) Copper Layers 

The copper layers were prepared in the evacuated electron diffraction 
camera, by evaporation of metallic copper from a tungsten filament. 
The ‘ specpure ’ copper was cut into suitably sized pieces, for subsequent 
insertion in the tungsten spiral, before being sent for exposure to neutrons 
in the pile. The evaporation was carried out with the substrate at room 
temperature, and the diffraction pattern from the surfaces was simul- 
taneously observed on the fluorescent screen. The arrangement of the 
specimen chamber has already been described (Kehoe, Newman and 
Pashley 1954); the vacuum was about 10~-*mm Hg, as measured by a 
hot filament ionization gauge. 

The rate of deposition was adjusted to a convenient value of about 
0-54/minute by variation of the current through the tungsten filament. 
Since the copper oxidizes when exposed to air, all the plates required 
from a given deposit were taken before the specimen was removed from 
the camera. The reloading of the plate holder was carried out while the 
vacuum was maintained immediately around the specimen by means of 
a special vacuum lock (Kehoe, Newman and Pashley 1954). 
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§ 3. THICKNESS MEASUREMENT 


Both isotopes used (Br 82 and Cu 64) are £8, y emitters, and in the 
present experiments the f-rays were counted. The relevant data for 
these isotopes are given in table 1. 


Table 1 A 
Isotope | Br 82 | Cu 64 
Maximum f energies (Mev) 0-45 0-66 
Half-life (hours) 35°5 12:8 
Initial specific activity (me/g) 110 50 
Initial counting sensitivity 
(counts/min/A above background) 110 180 


The counting apparatus consisted of a G.E.C. type GM4 counter, 
used in conjunction with a standard power pack and scaler. The counter 
was housed in a lead container which halved the background counting 
rate. The specimen was mounted immediately behind a circular aperture 
+ in. in diameter, which was placed directly in front of the counter window. 
In this way, all specimens were placed in exactly the same position, and 
their effective areas were equal. Since the back scattering corrections 
were also the same for all specimens, an accurate comparison between 
the thickness of any two specimens was possible. 

The absolute values of thickness were determined from a calibration, 
which was carried out after the isotopes had decayed by a factor of 
about a thousand. For the purpose of this calibration the ammonium 
bromide was converted into silver bromide, by precipitation from a 
silver nitrate solution. Thick layers (about 1000 A) of the silver bromide 
were deposited by evaporation on to clean glass slides (1 in. x 1 in.), and 
the amount of bromide in these layers was determined by means of a 
semi-microbalance. An identical thickness of bromide was simultaneously 
evaporated on to one of the silver specimens ; this specimen was used 
for the determination of the counting rate corresponding to the measured 
film thickness. The counting sensitivity could then be calculated for 
any time during the experiment. The same procedure was used for the 
copper, which was evaporated as the metal. 

The absolute accuracy was about +15°%. The relative accuracy was 
determined only by the number of counts registered and was better than 
+5%, except for layers less than 2 4 in thickness for which measurements 
were accurate to +0-1 4. 


§ 4. THe Errect or Deposit THICKNESS ON THE DIFFRACTION PATTERNS 


All diffraction photographs were recorded on Ilford thin film half tone 
plates (N50), which were developed for 2 minutes in ID2. The accelerating 
potential for the electrons was maintained at about 45 ky. 
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(a) Silver Bromide Layers 


A large number of specimens was prepared, ranging in thickness from 
about 0-24 to 204. After the thickness of the bromide layer had been 
measured, the specimen was examined in the diffraction camera and 
photographs were taken at various values of « (defined in §1). In all 
cases where the silver bromide was observed, it was oriented with its 
(111) planes parallel to the silver (111) surface. The orientation was 
partly fibrous with a preference for the [110] and [110] silver bromide 
directions to be parallel to the silver [110] direction. 

‘ Specimens covered by layers of mean thickness less than 0:44 gave 
rise to diffraction patterns which were indistinguishable from those 
obtained from the freshly prepared silver specimens (see fig. 1, Plate 15). 
Layers of mean thickness between 0-44 and 1-04 were observable by 
means of secondary diffraction patterns, such as that shown in fig. 2 
(Plate 15), the interpretation of which is given in fig. 3. Primary 
diffraction spots were observed from the bromide when the thickness 


Fig. 3 
Aue B A B 


Interpretation of fig. 2 (Plate 15). 


@ Spots due to the substrate silver; (111) surface, [110] azimuth. These 
spots are drawn out into almost continuous streaks due to the flatness 
of the surface, which restricts the penetration of the electrons to the 
top one or two atomic layers, when small values of « are used. 

x Spots which would appear due to the oriented silver bromide ; streaks A 
and B are due to secondary diffraction, the silver streaks acting as 
sources for diffraction by the bromide. 


was increased beyond 1:04. Figure 4 (Plate 15) shows such a pattern 
from a layer 2-6A in thickness, at a low value of «. As « was increased, 
the cross-grating (primary) spots became weaker and the secondary 
diffraction pattern started to appear (see fig. 5, Plate 16). At higher 
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Table 2. Results for the Silver Bromide Layers 


X X Intensity of silver pattern 
Ae Y Intensity of primary silver bromide pattern 
Z Z Intensity of secondary pattern 


R. C. NEWMAN & D. W. PASHLEY Phil. Mag. Ser. 7, Vol. 46, Pl. 15. 


Freshly prepared silver substrate. Silver substrate plus 0-6 A of silver bromide. 
(111) surface, [110] azimuth. a=1-8°, silver [110] azimuth. 


Silver substrate plus 2-6 A of silver bromide. 
a=1-0°, silver [110] azimuth. 


R. C. NEWMAN & D. W. PASHLEY Phil. Mag. Ser. 7, Vol. 46, PI. 


Fig. 5 


As for fig. 4, except a=2-2°. 


As for fig. 4, except a=3-4°. 


R. C. NEWMAN & D. W. PASHLEY Phils Magecers 7. Vol. 46, Pl 1/7. 


Fig. 7 Fig. 8 


> substrate plus 7-0 A of silver bromide. Silver substrate plus 15 A of copper. 
a—=1-2°, silver [110] azimuth. a=0-9°, silver [110] azimuth. 


Silver substrate plus 0-9 A of copper. 
a—0-9°, silver [110] azimuth. 


R. C. NEWMAN & D. W. PASHLEY Phil. Mag. Ser. 7, Vol. 46, Pl. 18. 


Silver substrate plus 3-3 A of copper. 
a=1°, silver [110] azimuth. 


As for fig. 11, except a=2-6°. 
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values of «, as shown in fig. 6 (Plate 16), the secondary diffraction pattern 
predominated over the primary pattern, and a strong silver pattern 
was observed. With further increase of thickness, the primary diffraction 
pattern from the bromide became quite strong (see fig. 7, Plate 17), and 
was visible at all values of « up to 5°. These results are presented in 
table 2. The intensities of the components of the pattern are estimated 
qualitatively by visual inspection, and are denoted in order of decreasing 
intensity by means of the letters A-H. O indicates the absence of a 
particular component. The table is intended to show only the character 
of the variations of the diffraction pattern with thickness and «, and is 
not quantitative. 

The silver substrates give rise to quite prominent Kikuchi patterns, 
which become slightly less distinct as the bromide is grown on the 
specimens. Even with the thickest layers (204), these Kikuchi patterns 
are still easily distinguished. 

Since prominent cross-grating spots are observed from bromide layers 
as thin as 1-24, the deposits do not grow as a uniform film, but in the 
form of a large number of nuclei. Hence only a fraction of the surface 
is covered. Attempts have been made to estimate the size and shape 
of these nuclei. The diffraction patterns show no definite refraction 
effects, so that the nuclei do not possess atomically smooth surfaces. 
The measurement of the size of nuclei is always uncertain, particularly 
since the size of the diffraction spots varies considerably with photographic 
exposure. The size of the focused central spot on the photographic plate 
was small compared with the size of the diffraction spots, so that the 
length of the diffraction spots in any direction can be used for the 
measurement of the thickness of the nucleus in that direction. It has 
been assumed that there is no other cause of appreciable spot broadening, 
such as ‘ tailing ’, for which there was no evidence. Measurements gave 
for the average size a height of about 404 and a width of about 704 ; 
it was not possible to detect any change in the size of the nuclei as the 
thickness of the layer was increased. 

The lattice spacing of the bromide layers has been measured relative 
to that of the substrate silver. Because of the nature of the patterns, 
it was possible to measure spacings only of planes perpendicular to the 
_ surface. From measurements of the primary pattern on a large number 
of plates, the value of 5-737--0-018A was obtained for a). This is 0-5% 
lower than the x-ray value of 5-7684. In this case the measurements 
from the various plates were quite consistent, within the above accuracy 
which was limited by the size and shape of the diffraction spots from the 
silver bromide. When the secondary diffraction patterns were measured, 
however, satisfactory consistency was not obtained. Most specimens 
gave values which were in agreement with the above value for the 
primary measurements, but some gave much higher values (by as much 
as 2°/) which varied from plate to plate taken from the same specimen. 
The reason for these variations is not known. 
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(b) Copper Layers 

A number of specimens was prepared ranging in mean thickness from 
0-44 to 204. Initially the deposited copper was always oriented with 
its (111) planes parallel to the silver surface, and with the [110] and 
[110] copper directions parallel to the [110] silver direction. As the 
thickness was increased some randomly oriented copper was simultaneously 
formed, being just visible as broad rings at a mean thickness of about 5A. 
When the deposit thickness was increased to about 15 A, the rings became 
more distinct and had approximately the same intensity as that of the 
spot pattern (see fig. 8, Plate 17). 

Layers of less than 0-74 in mean thickness gave rise to no extra 
features in the diffraction patterns from the substrate silver. At about 
0:84 elongated spots due to primary diffraction by copper nuclei appear. 
These spots increase rapidly in intensity for further small increases in 
thickness. A typical example is shown in fig. 9 (Plate 17); the inter- 
pretation is given in fig. 10. At about 5A, the silver pattern is almost 


Fig. 10 
4 ' e @ 
¢ é : : 
r) i e @ 
s 


Interpretation of fig. 9 (Plate 17). 


e Spots due to the substrate silver; (111) surface, [110] azimuth. These 
spots are drawn out into almost continuous streaks because of the 
flatness of the surface, as for fig. 3. 


§ Spots due to copper ; (111) surface, [110] and [110] azimuths. 


completely obscured for all values of « less than about 2°. Most 
specimens of thickness less than 5A gave rise to secondary diffraction 
effects ; they were much less prominent than in the case of the bromide 
layers. These results are given in table 3, where the same scheme is 


Means of Detecting Thin Surface Films : I 935 


used as in table 2. At «=1°, quite a thin layer of copper causes the 
diffraction pattern from the silver substrate to be almost obscured (see 
fig. 11, Plate 18). The copper pattern consists of a cross-grating of strong 


Table 3. Results for the Copper Layers 


a O'S 0-7 1-0 1:3 1-6 2-0 2-5 3:0 
| 
Thick- 
ness (A) 
aa A A A A Fiea ems marry 
0-8 A A A B B C D F 
O O L®) H H H G G 
B A A x A A eae tens 
1-0 A A A B B C C D 
O O O H H G G F 
E D D C B A A A 
1-3 A A A A B B B B 
O H G G G G G E 
F E D C B B A A 
2-0 A A A A A A B B 
H G G F F D D D 
G F F E E D B A 
2-5 A A A A A A A B 
sl H H G G F E D 
HH G F in E E C C 
3:3 A A A A A A A A 
H ist H H G G F D 
O O O O O O O G 
Tet A A A A A A A A 
O O O O O O O H 
x X Intensity of silver pattern 
Na Y Intensity of primary copper pattern 
Z Z Intensity of secondary pattern 


elongated spots. Some of these spots are seen to be split into pairs in 
a direction parallel to the shadow edge. This fine structure fits neither 
with the splitting expected from the occurrence of reciprocal lattice 
spikes (Laue 1936, 1937) arising from crystal shape effects, nor with 
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effects due to refraction. No satisfactory explanation has been found. 
A possible explanation seems to be that dislocations are present in the 
copper nuclei. A more detailed examination would be necessary to 
check this point ; it is beyond the scope of the present work. An increase 
of « causes the silver pattern to become stronger, as shown in fig. 12 
(Plate 18). 

The Kikuchi patterns from the silver substrate became gradually 
less distinct as the thickness of the copper was increased. At a thickness 
of about 204 the patterns were completely obscured for «<2°. 

The diffraction spots from the copper nuclei were not in the normal 
Laue positions, but were displaced towards the shadow edge of the 
patterns. This shift is due to refraction of the electrons during entry 
and exit at the crystallite surfaces. The magnitude of the shift was 
less than the full value to be expected for an inner potential of 14 volts. 
which is the value calculated by Tull (1951). This would result from an 
averaging effect, due to the possibility of the electrons entering the 
nuclei by one facet and leaving by another. A displacement occurs only 
as the electrons pass through the top surface, which is therefore very 
smooth on an atomic scale. From the size and shape of the diffraction 
spots it is deduced that the nuclei are plate-like, with a height of about 
124A and lateral dimensions of about 604. The measurement of the 
height must be taken as a lower limit, because part of the elongation 
of the spots might arise from variable refraction effects. No variation 
of spot length was observed as « was changed, so that it was deduced 
that the electrons were penetrating the whole depth of the nuclei. 
Halliday, Rymer and Wright (1954) have recently studied thin films of 
copper prepared by vacuum evaporation. They deduce that the films 
are composed of nuclei which are plate-like, with a height of 15 A, and a 
diameter of 90 A. 

A nuclear thickness of 124 implies that, for the thinnest detectable 
layer of 0-8A, the proportion of the surface covered by the nuclei is 
about 7°%. Measurements showed that the size of the nuclei did not 
increase appreciably as further copper was deposited, but that the 
number of nuclei increased. It is interesting to note that the non- 
oriented copper started to appear only after a large fraction (say >40°) 
of the silver surface had become covered by nuclei. 

The lattice spacing of the copper nuclei was measured relative to 
that of the silver substrate. A mean of several measurements gave 
Aqy=3°635-+-0-015 A, which is about 0-75°% greater than the x-ray 
value of 3-608. A similar expansion of the copper lattice in thin films 
has previously been observed (see Shishakov 1952). 

All these results were found to be independent of the rate of deposition 
of the copper, over the range of 4A/hour to 80A/hour, the normal 
deposition rate being about 25A/hour. It was found that the rate of 
deposition could be maintained at a constant value for a period of at 
least one hour. 
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§ 5. Discussion 


The present experiments have shown that the radioactive method, 
which requires no elaborate or laborious technique, is highly accurate 
for the measurement of the thickness of surface layers as thin as 0-14. 
In addition, the method allows independent measurements to be made 
on each specimen as distinct from methods which involve only the 
calibration of the rate of deposition. As a general method, however, 
it has the following limitations : 

(1) many elements do not have a suitable isotope, 

(2) an isotope with a half-life of a few hours is required where 
permanent apparatus is contaminated, 

(3) if deposits as thin as 0-14 are to be measured on specimens of the 
size normally used in electron diffraction experiments, the specific 
activity of the material should be at least 5 mec/g; this ensures that 
the counting time for each specimen is not prohibitively long. It is 
preferable to count B- or y-particles, rather than «-particles, in order to 
minimize effects due to particle absorption in the deposit. The dis- 
advantage of the use of «-particles is illustrated by the inconsistency of 
the counting rate, for a given film thickness, found by Antal and Weber 
(1952). 

By means of the radioactive method described in §3 it has been 
conclusively shown that layers with an average thickness of less than 1 A 
can give rise to observable electron diffraction patterns. These results 
apply to electrons with energy of about 45 kev; the sensitivity of 
detection is expected to vary with changes in this energy, but it is 
difficult to predict the form of this variation. 

The lower limit of detection is 0-84 for the copper layers, and only 
0-44 for the bromide layers. A difference of this order of magnitude 
is to be expected because the scattering power per unit volume of the 
bromide is about one and a half times that of the copper, over the 
employed range of scattering angles. This high degree of sensitivity 
agrees with the results obtained by Schulz (1951, 1952 a) and thus 
Seifert’s doubts (see Gomer and Smith 1953, pp. 360 and 383) are not 
justified. 

An interesting difference occurs between the diffraction patterns from 
the copper and the silver bromide layers. From tables 2 and 3 it is 
seen that for a deposit thickness of 14, and for small values of «, the 
bromide layer gives no primary diffraction, whereas the copper gives a 
very strong primary pattern. This is a surprising result, since it would 
be expected that with equal volumes of the two deposits the bromide 
would give rise to stronger primary diffraction patterns than the copper, 
provided that (a) all of the nuclei have their normal crystal structure 
and orientation, (b) the nuclei are sufficiently small for multiple scattering 
to be negligible, and (c) none of the nuclei is shielded from the incident 
electron beam. Since conditions (a) and (b) are considered to be 
satisfied, the above difference appears to be due to condition (c) not 
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being -fulfilled. A likely explanation is that in the initial stages the 
silver bromide nuclei are partially embedded in the silver substrate. 
This could occur as a result of the mechanism of the chemical reaction, 
which might be expected to lead to the diffusion of bromine into the 
surface layers of the silver. It might also be due to the occurrence of 
growth in holes, grain boundaries, steps or cracks in the surface. Under 
such conditions there would be a high probability of the electrons being 
diffracted by metallic silver before or after entering the silver bromide 
nuclei. This would cause the intensity of the primary diffraction pattern 
to be reduced and that of the secondary diffraction pattern to be 
increased. This interpretation would, therefore, be consistent with the 
observation of very strong secondary patterns from the bromide layers, 
compared with those from the copper deposits (cf. tables 2 and 3). 

These results will now be compared with those of Schulz (1951, 1952 a). 
Schulz grew alkali halides, by condensation in vacuo, on to cleavage faces 
of various single crystals, and according to his interpretation the layers 
were of two types: (i) monolayers, (ii) nucleated deposits. The mono- 
layers were most easily detected at low values of «, the patterns being 
attributed to primary diffraction. He found, in the case of the nucleated 
deposits, that a more sensitive method of detection was to use a strong 
diffracted beam from the substrate crystal as a source for diffraction by 
the deposit nuclei. This latter method (Schulz 1950) required the use 
of a somewhat higher value of «. In the present experiments, however, 
although the deposits grew as nuclei, it was found in all cases that the 
greatest sensitivity of detection was obtained at the smallest possible 
values of «. 

In the case of copper layers, the first indication of the presence of a 
deposit consisted of definite primary diffraction spots. The smaller the 
value of «, the more intense were these primary spots, and secondary 
diffraction spots occurred comparatively weakly only at higher values 
of « (cf. table 3). It has been shown in § 4 (b) that the copper nuclei 
are plate-like in shape, and hence the copper deposits resemble Schulz’s 
monolayers more closely than his nucleated deposits. This is consistent 
with the primary diffraction pattern being predominant. 

In the case of the bromide layers, the first indication of the presence 
of the overgrowth consisted of extra streaks on the diffraction pattern 
(see figs. 2 (Plate 15) and 3). Streaks of type B arise only as a result 
of secondary diffraction, as illustrated in fig. 3, while streaks of type A 
could be due to a superposition of both primary and secondary diffractions. 
Since type B streaks were usually of the same intensity as those of type A, 
it is reasonable to conclude that type A streaks, and hence all of the 
streaks, are of secondary origin. As shown in §4(a) the bromide 
grows in the form of nuclei, and streaks arise in the secondary diffraction 
pattern simply because the silver spots acting as the sources are them- 
selves extremely elongated. The bromide layers are similar to Schulz’s 
nucleated deposits, and so their initial appearance by means of secondary 
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diffraction is to be expected from his results. The fact that Schulz used 
larger values of « than those found most suitable in the present work 
might be due to the difference in the diffraction patterns from the 
substrate ; the silver specimens give very good diffraction patterns even 
at the lowest possible values of «. 

Since, in Schulz’s experiments, it is reasonable to suppose that none 
of his deposits was embedded in the substrate, it follows that the 
appearance of secondary diffraction effects from the silver bromide layers 
is not sufficient evidence of the embedding of the nuclei. It appears 
that the secondary diffraction will appear anyway, provided the nuclei 
fulfil certain conditions on size and shape. However, it is considered, in 
view of other reasons given above, that the embedding of the silver 
bromide nuclei does occur and that there is more than one condition 
favourable to the appearance of the secondary diffraction. 

To summarize, it is concluded that the copper nuclei are plate-like 
and form on the surface of the silver. The measurement of the nuclear 
size is uncertain because of refraction effects. The silver bromide nuclei 
have no well defined shape, and are partially embedded in the metallic 
silver. The measurement of their size is not reliable because of the 
embedding. Since the primary diffraction arises mainly from the pro- 
truding part of the nucleus, the measurement of the spot size will include 
only that part of the nucleus. It will not measure the part of the nucleus 
embedded in the silver. 

In the case of the thicker deposits, a layer of 7-7 A of copper completely 
obscures the silver substrate pattern for values of «<3°. However, such 
an effective cut-out is not produced even by 194 of silver bromide, 
unless « is reduced to less than 14° (cf. tables 2 and 3). This large 
difference proves that the shape and size of the deposit nuclei are 
important factors in the determination of the conditions for the cut-out 
of the substrate pattern. These factors have not been considered by 
Raether (1949, 1950, 1951). His method of thickness measurement 
applies to deposits in the form of uniform layers, but very few deposits 
grow in this manner, and hence the method has very limited application. 
A theoretical treatment of the diffraction patterns to be expected from 
deposits on various types of surface would be extremely useful. Such 
a theory would have to consider the formation of both primary and 
secondary diffraction patterns, and also the effect of nuclear size and 
shape. 

The measurement of the lattice spacing of the deposit during the 
earliest stages of growth is of importance in connection with theories 
of orientation. By considering the initial formation of an oriented 
deposit to be in the form of a monolayer, Frank and van der Merwe 
(1949, 1949-50) deduced that the lattice spacing of this first layer will 
be modified to fit that of the substrate. They maintained that the 
formation of such a monolayer was an essential first stage for the 
occurrence of an oriented overgrowth. In the present experiments, 
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however, the deposits have not grown as extensive monolayers, but as 
oriented three-dimensional nuclei. As the deposit thickness is increased, 
the size of the nuclei remains approximately unchanged, and the number 
of nuclei increases. Thus the various stages of growth of a single nucleus 
are not observed, even for deposits of very small mean thickness. There 
is, therefore, no direct evidence concerning the lattice spacing of a single 
nucleus during the earliest stages of its growth. All that can be stated 
is that the lattice spacing of the completed nuclei, although slightly 
modified from the value for the bulk material, is not constrained to fit 
that of the substrate. If, however, a deposit grows in the form of a 
uniform oriented ‘monolayer’ (see Schulz 1951, 1952a, Gomer and 
Smith 1953, p. 372), the electron diffraction method is sufficiently sensitive 
to allow the first atomic layer to be detected, and its lattice spacing 
could then be measured. 
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ABSTRACT 


The graphical rules of Feynman and Dyson are stated in a very 
general form, which depends only on the field Lagrangian, and does 
not explicitly require either the use of field operators or the construction 
of a Hamiltonian scheme. It is shown that by choosing these rules as 
a starting-point one is able to quantize the electromagnetic field in a 
simpler way, since there is no longer any need to modify the classical 
Lagrangian {fF ,, (which does not lead to a Hamiltonian as it stands) ; 
to impose a supplementary condition ; or to employ an indefinite metric. 
Dyson’s theorem on photon propagation functions is established, and the 
gauge-invariance of the S-matrix is discussed. 


§ 1. LyrRopucTION 


General Statement of the Feynman—Dyson Rules 


Ir seems to be true that nearly all the techniques so far developed 
in quantum field theory that are not purely formal, but likely to 
be of direct interest in practical calculations, can be obtained in a 
straightforward way from the Feynman—Dyson rules of S-matrix theory 
(Feynman 1949 a, b, Dyson 1949). (These techniques would include, 
for instance, the one-particle Dirac equation with radiative corrections, 
the Bethe-Salpeter equation, and probably also the various Tamm-— 
Dancoff equations.) In other words, a knowledge of these rules is all 
that is really required to enable any physical quantity to be calculated. 

It is also true that the Feynman—Dyson rules can be written down 
at once in a very simple form, for any particular system of interacting 
fields, as soon as the Lagrangian # is known. Suppose that # (expressed 
as a function of the field operators ¢%, their conjugates ¢,, and derivatives) 
has been divided into two parts - 


Gi+ Gi, ah seh eee ae nd 0 (1) 


where # is bilinear and represents the free fields, while Y’ describes the 
interactions. In order to avoid any explicit consideration of the form 


* Communicated by the Author. 
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of Y, it is convenient to rewrite each part of the corresponding action 
function J in such a way that 


=| 2) de=55|| $,(x)Lg(x, y)p*(y) dwdy . . . (2) 


(one term for each free field), and 


1 
=| 2 ie=z (1175) NPEAL de 
(x) ; d;! (i. | Bale )ba(%2 
x V8G (a, % a+ 65 Yrs Yar JP(Yr)P' (Yo). -- dey darg...dy,dyg-.. . (3) 


(one term for each interaction). This defines a set of kernels L,*, V35-: 
containing 5-functions and their derivatives, together with Dirac eee 
and fundamental constants. (For a typical example cf. eqn. (20) below.) 
The numerical factors d, d;! are designed to take into account the 
possibility that a given field operator may be raised to a higher power 
than the first: thus d=2 for self-conjugate fields (photon, neutral 
meson) ; d=1 for non-self-conjugate fields (representing charged particles, 
neutrons etc.) ; d; is the degree to which the jth field operator (together 
with its derivatives) occurs in the corresponding term of #’. The 
transformation (2) allows the variational equations for each of the 
free fields to be expressed in a standard ‘ integral ’ form 


| Lee DPW) dy=o, | $.(0) (x)Lgt(e, y)da=0, . . (4) 


which has the advantage of being valid for differential equations of 
arbitrary order (and also for integral equations, of course). Equation (4) 
may be written more concisely as 


Ld—0, “¢L==0,. “ J ae 


thus defining a linear operator L. 

Since the Feynman—Dyson rules are already well known, in order to 
bring out their dependence on the Lagrangian it should be sufficient to 
remark that the correct propagation function Q, (say) for each line can 
be defined by the formal limiting process 


—ith 
Q,= 
re) ae L—te’ 


(6) 


where «>0, and therefore satisfies the inhomogeneous equations associated 
with (5), 
LQ,=Q,L=—izh, » Aly et ee 


while the ee for each vertex is simply (ih)-1 times the corresponding 
operator V%4---, which acts on the propagation functions or wave-functions 
belonging te the lines just as on the field operators in (3). The definition (6) 
is equivalent to Feynman’s introduction of a negative-imaginary mass 
increment (Feynman 1949 a, § 6), which ensures that positive- frequency 
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waves are always radiated into the future, while negative-frequency 
waves are radiated into the past. (For further details about these rules 
the reader is asked to consult the papers of Dyson (1949). The fact 
that the vertices really depend on #’, rather than on the interaction 
Hamiltonian H* (which may contain extra terms), first. became apparent 
from the work of Matthews (1949).) 

In using this Lagrangian version of the rules one is able to circumvent 
nearly all the abstract formalism of quantum field theory (such as 
canonical coordinates and momenta, commutation rules, state-vector, 
Hamiltonian, Heisenberg and Schrédinger equations, definition of the 
vacuum, etc.); all that matters in practice is contained in the 
Lagrangian # and in the definition (6). Furthermore, even the field 
operators themselves do not play any essential part, for what one really 
needs to extract from # is simply the set of kernels L,*, V%:.. Of 
course the Lagrangian version can only be justified at present by deducing 
it from the usual Hamiltonian scheme, but such a simple and apparently 
universal formulation of relativistic quantum theory must have a deeper 
significance than has yet become clear. In this connection it is very 
instructive to examine a fairly complicated type of field, such as that 
of the vector meson, where the construction of a Hamiltonian scheme 
encounters the problem of the ‘ missing momentum ’, conjugate to the 
time component U° of the meson potential. In spite of the various 
manipulations that are necessary if the theory is to be expressed in a 
covariant Hamiltonian form (e.g. one might begin by eliminating the 
time-like component U“n,(c), normal to an arbitrary space-like surface ¢), 
the Lagrangian version of the rules is eventually found to be quite 
unchanged. This is a remarkable fact, and it suggests that the same 
rules might even be used when no Hamiltonian scheme exists at all. 


Application to the Electromagnetic Field 


A natural application for this idea is provided by the electromagnetic 
field. The most obvious choice of Lagrangian in classical theory is 


emi Ee (8) 

GA” dAt 
(ss = a a 9 
where F ame ni (9) 
because this leads directly to Maxwell’s equations.* The corresponding 
energy density is POs(HILE), . . .. . . + (10) 


Unfortunately, there is no Hamiltonian associated with (8), since the 
momentum conjugate to A® does not exist, and the usual procedure in 
setting up a canonical scheme for quantization has been to modify the 
original Lagrangian to 


eae ae nae 


Ont! 


* Upper indices are assumed to be contravariant, and. 91:=922=933—=!, 
Io=—l1; [J=—O?/Oxdx,, c=1. 
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so that the variational equations become 
[AF 0. Sh ee ad ae oe, eee 


Maxwell’s equations are then restored by imposing a supplementary 
condition, whose classical formulation would be 
one 
One 
However, the modification (11) introduces a number of serious diffi- 
culties, which have no analogue in the original version of the theory. 
The energy density is now 
4{(curl A)? +A2+( div A)? — gi (prad p)*}) 5. eee ee 
which is not positive definite, and the free field has four degrees of 
freedom instead of two; in quantum theory there are scalar and 
longitudinal as well as transverse photons, and the scalar photons appear 
to have negative energy. In particular this makes it difficult to define 
the free vacuum state in an invariant way ; normally one would expect 
to take 


(13) 


Att <0 Hees ig ne ne 


(where + denotes the positive-frequency part), but for negative-energy 
scalar photons A®+ would be an emission operator. This in turn reacts 
on Dyson’s definition of the photon propagation function as 

Dy P(x, y= PoP(A(x)A(y)) Po, - - - + + (16) 
where P is the chronological ordering operator (Dyson 1949). 

The most elegant treatment of these problems has been given by 
Gupta (1950) and Bleuler (1950), who avoid the occurrence of negative 
energies by introducing an indefinite metric in the Hilbert space of 
quantum states. This device allows a consistent definition of the 
vacuum, so that (16) can be given a rigorous meaning, and also makes 
it possible to prove a theorem of Dyson (1950), which states that the 
expression (16) may always be replaced by the simpler propagation 


function Dae (x; 9) = SPD e—y)e =. 
where CID (2) == — 1802). oe a ee 
These authors also show that the correct quantum analogue of the 
classical condition (13) is 

0AM + 

al Het A . : . . - . . (19) 


(in the Heisenberg picture). 

In spite of the formal success of the Gupta—Bleuler method, it is hard 
to give any physical meaning either to the difficulties that arise from 
the modification (11), or to the technique used for overcoming them, 
especially since the introduction of an indefinite metric leads to ‘ negative 
probabilities’. On the other hand it is known that the original 
Lagrangian (8) is perfectly satisfactory as it stands in the classical 
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theory, because the Hamiltonian scheme is never really needed. One 
would expect it to be equally satisfactory in a Lagrangian version of 
quantum theory. The main task of this paper will therefore be to apply 
our version of the Feynman—Dyson rules to the electromagnetic field in 
the most straightforward way, i.e. using (8) instead of (11), and to show 
that this procedure leads to consistent and intelligible results, while 
avoiding the difficulties encountered in the usual method. 

This programme will be carried out in § 3. However, the variational 
equations associated with (8) exhibit some unusual features, which it 
is necessary to discuss, and so these equations are dealt with in § 2. 


§ 2. STRUCTURE OF THE VARIATIONAL EQUATIONS 
The ‘ integral operator ’ belonging to (8) is 
0 0 0 0 
p py ES ASS cae 0 Cs Dae rs _ (20 
Lele, N= (s-BF— ast) (sede 38" )dle—y), - (20) 
leading to the equations 


- akg 
P S — Pp a ig ==(()). ae, al 
| Be NAW) dy = (BP+ ares) 4°) =0 (21) 
The associated inhomogeneous equations are 
oO? 
P eee eet be . . . . . . 22 
(ce. , ae 6 ee) 
which of course are equivalent to Maxwell’s equations 
CEs (23) 
Cate 


The pecularities of (22) are conveniently studied by making a Fourier 
transformation : 

LE(k)A(k) = (h75,P—k?k,)A(k)=7°(k), poy (24) 

FP (kh) = 1(kPA°(k)—k° AP(k)) ee tn oe) 

where P=k'k,. LS(k) is a 4X4 matrix, but it has rank 3, since the 
four rows and columns are linearly dependent : 

kL P(k)=0, Lek ee =O; eee ae (20) 

Equation (24) is therefore a set of inhomogeneous algebraic linear equations 

with vanishing determinant, and in accordance with the usual algebraic 


theory this set will be incompatible, and thus have no solution, unless 
the vector j°(k) lies in the 3-dimensional subspace V, defined by 


Peter) ee - (27) 
When (27) is satisfied (24) must possess an infinite set of solutions, since 
to any one solution A°(k) it is possible to add an arbitrary multiple of 


the vector k?: AP(k)--a(k)ke. Feet si oe | CAMs OMS (28) 


However, it is clear from (25) that this arbitrariness cannot affect Pel), 
and it may be removed by imposing the condition that A°(k) is to lie 
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in the same 3-dimensional subspace V, as j°(k), i.e. that it is to be 


orthogonal to k,, : re 
6 p k,A°(k)=0, eee 
so that «(k)=0. After reduction to V, (24) becomes 
LiP(k)A°(k) = (h?8,°)A%(ke)=7°(k), eres (::!) 


with the unique solution 
AP (Kk )e=f O()/RAy cay Al eee eo 
unless of course k?=0. 

It should be remarked that the orthogonality condition (29) also 
fails to define A°(k) uniquely when k?=0 (the Lorentz metric is of course 
indefinite), so that «(k)40, and in addition to the arbitrariness in F°’(/) 
there is a further arbitrariness in the relation between F°’(k) and A®(k). 

The algebraic method described here is obviously the exact analogue 
of the usual technique in which the differential eqn. (24) is simplified 
by imposing the Lorentz condition 0A“/dx“=0, but the latter seems at 
first sight to be a rather ad hoc procedure; it should become much 
clearer when related to the familiar theory of inhomogeneous algebraic 
equations. Thus (27) is equivalent to the conservation of charge ; 
(28) corresponds to the existence of unrestricted gauge transformations 


AP -> AP as : (¢ arbitrary) ; CTS Soest 


oes 
fe) 

(29) to the Lorentz condition ; (30) to the usual wave-equation [_]AP=7?, 

and the remaining arbitrariness in A®(k) to the existence of restricted 


gauge transformations of the form 
AP -> AP+ ad, (F}p==0). ep er 
Zc 


Finally, it should be pointed out that eqn. (22) is not as strange as 
it might seem; the differential operator on the left-hand side can be 
given a simple interpretation. A 4-vector A® would normally be expected 
to split up into two components on restriction to the 3-dimensional 
orthogonal group, belonging to representations with dimensions 1 (scalar) 
and 3 (vector). L,°(k) is in fact the projection operator on to the sub- 
space V, (multiplied by k?), and its purpose is to avoid the unwanted 
scalar representation automatically. By changing the Lagrangian (8) 
to (11) one replaces this projection operator by the identity, and so 
brings in a scalar component, which must then be removed by a suitable 
supplementary condition. 


§ 3. QuantuM THEORY 


A formal application of Feynman’s rule (6) to the operator (20), 
defined by the Lagrangian (8), leads to 


2 
(C+ aoe -id?) Do Ax, yy=—ihd PS(xa—y). . . (84) 
p 
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This equation can be solved by first differentiating with respect to «°, 
which gives 


Camner hoa 
Rage Del&, Y= = = 8(e—y). etre, JAPA 35) 


The second term on the left-hand side can then be transferred to the 
right, and one finds 


: & eke 
(-1D3.l2, =—i8 (B+ ara BU, BO) 
and therefore 
DE y= ee is D 
BNE J Oi ax,dy" R(L—Y), eee to) 


where Dy, is defined by (18). 


Dyson’s Theorem 


It is clear that Feynman’s introduction of a negative-imaginary mass 
increment has also succeeded in raising the rank of the operator L,°(a, y) 
to 4, so that it now has a unique inverse. It is true that the propagation 
function D§,(x, y) defined by (37) becomes singular in the limit «0, 
but as Dyson (1950) has remarked, in S-matrix theory this function will 
always act on currents J,”(y, x) satisfying 

OS eed 

0x, Oy" ney 

p 

(when summed over all graphs belonging to the same order in ¢). If 
one therefore makes the rule that elements of the S-matrix are to be 
evaluated by multiplying the propagation functions together and summing 
over all graphs before the limiting process «+0 is carried out, the second 
term on the right-hand side of (37) will give zero on integration by parts, 
and can be omitted altogether. Hence it is permissible to replace (37) 
by the modified propagation function D,6,? that is normally used in 
electrodynamics. This is just the content of Dyson’s theorem. 


(38) 


Photon Wave-Functions 
It is convenient to imagine that the initial and final photon wave- 
functions .V? associated with a given collision process are ‘ radiated ’ 
from sources in the past and future respectively, localized in the small 
regions of space-time where the initial photons were originally created, 
or the final photons will eventually be destroyed. These wave-functions 
will then be defined by 


r 1 Ge : 
A(a)= | Dios wy) dy = 2 | (3° 7 aa) Daley). (89) 
Since 0.J7/dy7=0, (39) is equivalent (after integration by parts) to 
b(a)= | Dyle—y)I*y) dy, (40) 


and from this it follows that 0.°7°/da°=0. 
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There can be no arbitrariness associated with J°, and so one may 
conclude that each individual photon has a wniquely defined ‘natural ’ 
wave-function or potential .»/?, and that this potential satisfies the 
Lorentz condition (as Feynman (1949 b, § 8) has pointed out). Feynman’s 
prescription (6) is therefore more powerful than the classical Lorentz 
condition by itself, which leaves the potential arbitrary to within a 
restricted gauge transformation of the form (33). (Of course, even in 
classical theory the arbitrariness (33) is removed if one insists on using 
retarded potentials, for example.) 


Gauge-Invariance of the S-Matrix 


Although each of the initial and final photons associated with a given 
collision has a unique natural potential which is defined by (40) (and 
therefore depends on conditions in some other region of space-time, 
where the photon is created or destroyed), in the neighbourhood of the 
collision itself each photon will propagate like a plane wave satisfying 
the homogeneous equation, and there will be no way of telling what its 
natural potential may be. In fact one can make a quite unrestricted 
gauge transformation of the form (32) on the wave-function for each 
photon separately, without altering the amplitude for any collision 
process, since in evaluating the S-matrix each photon wave-function is 
multiplied by a current J, satisfying 0J,/dx,=0, so that 


a, (a a7 ace 
[rege 2= | Fe, ob) de lace dx=0. 


The state of a free photon is therefore completely defined by its ‘ field ’ 
FP? (related to ? according to (9)), and any wave-function consistent 
with this field may be used in calculations. It follows from (40) 
(together with the conservation of charge) that °° satisfies Maxwell’s 
equations, so that at large distances from the source 6 and #& are at 
right-angles to one another and to the direction of propagation, as in 
the classical theory. 


Energy-Momentum Tensor 


Since the graphical method of quantization described in this paper 
does not use the concept of field operators A“(x), it is no longer possible 
to construct a canonical energy-momentum tensor by the usual varia- 
tional principle, and the simplest way of defining the energy and 
momentum associated with a free photon is to vary the metrical tensor g,,. 
It is well known that in classical field theory a symmetrical energy— 
momentum tensor can be defined by 

Ogee 
T(z) 72) eee ee (41) 
where J is the action function, and y,, is equated to its value for flat 
space-time after the variation has been carried out. In quantum theory 
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the effect of a change in GJur(z) will be to introduce a new type of 3-vertex , 
with kernal 

a Oy 2) NS ks) res ee AE a) 
say, describing the scattering of a photon by the new ‘external field ’ 
at the point z. Equation (42) may be evaluated by remarking that 
any small external interaction will modify the photon propagation 
function from D, to D,+5D,, where 


—h 
sD, =5 (=) =L1(thdL)L~1= D,[(ih)-15L]D,. . (48) 
This is equivalent to introducing a new vertex with kernal SL, so that 
) 
—TEe (ey Z)= amr LP(a, y). er, oer nee (44. 
(vy 3 2) Jal) (x, y) (44) 


(The interaction terms in the energy-momentum tensor can be obtained 
in a similar way by functional differentiation of the vertex kernels or 
V-functions, but these are not of interest for a free photon.) 

The result of the variation (44) is most simply stated in the form 


T(z) =| | Ayle\TQM(w, ys 2)f%(y) dee dy 


= Fe (F(z) + FF f (2) IFO (2)F sglz)g”, - (45) 


a 
which corresponds to the classical energy-momentum tensor, except for 
an extra factor 2. The expectation value of the energy density associated 
with a free photon is therefore 


PACH Oe Whe gees 2 (46) 


By integrating (45) over a space-like surface one can find the total flux 
of energy and momentum P“ associated with the photon. This enables 
the wave-function to be normalized correctly. 

Equation (45) is invariant under unrestricted gauge transformations 
of the form (32), since it depends only on the field 4°’. Hence, there 
is no energy or momentum associated with a potential d¢/éx,. This 
might be explained by saying that the energy-momentum operator 
T?""(k) acts only on the subspace V, of § 2. 


§ 4. DIscUSSION 


It has been shown that by starting from a Lagrangian version of the 
Feynman—Dyson rules (§ 1), one is able to develop a much simpler 
treatment of the quantized electromagnetic field than has hitherto been 
available within the framework of the Hamiltonian theory. The reason 
for this is that if one uses Feynman’s definition (6) for the propagation 
function, rather than Dyson’s definition (16), it is no longer necessary 
to construct field operators, commutation rules or the vacuum state 
explicitly, and so it becomes possible to work with a more general kind 
of Lagrangian, for which a Hamiltonian scheme would not exist. The 
classical Lagrangian }F”’F,, happens to be of this kind, and turns out 
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to be extremely well adapted for describing particles of spin 1 and 
mass 0 (§ 2). By avoiding the necessity for modifying this Lagrangian 
one greatly simplifies the gauge-invariance problem. 

On the other hand, since the Hamiltonian formalism of relativistic 
quantum mechanics has been dispensed with altogether in this paper, 
and a rather abstract method of Lagrangian quantization used instead, 
it is no longer possible to rely on general theorems in quantum mechanics 
for establishing either the physical interpretation of the self-consistency 
of our treatment of the electromagnetic field, and some explicit investi- 
gation ought to be made. Actually so far as the electromagnetic field 
itself is concerned, the new treatment can be justified for all practical 
purposes by comparing it with the more complicated Gupta—Bleuler 
method, which leads to the same numerical results. But it would 
obviously be much better to establish the general validity of Lagrangian 
quantization by some independent argument. 

From this point of view, it may be worth remarking that a new 
formulation of relativistic quantum theory can be developed, in which 
the Feynman—Dyson rules themselves are chosen as the axioms. I call 
this the Theory of Graphs (TG), since it starts from the basic idea of 
‘graphs’ rather than ‘fields’ or ‘elementary particles’. Roughly 
speaking, it might be said that TG represents an attempt to provide the 
simplest posssible logical framework for the intuitive graphical methods 
of Feynman (1949 a, b), so that these can be justified directly, without 
appealing either to the Hamiltonian formalism of quantum mechanics, 
or to the more classical concepts of ‘ continuous particle paths’ (Feynman 
1948) or ‘ continuous field distributions ’. TG is equivalent to quantum 
mechanics, but there are differences both in the mathematical formalism 
(which is Lagrangian rather than Hamiltonian), and also in the physical 
interpretation. It is found that a number of problems can be treated 
more simply with the help of TG ; in addition to gauge-invariance, these 
include time reversal, charge conjugation and the derivation of relativistic 
wave-equations. I hope to publish a paper about this theory shortly. 
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ABSTRACT 


The isotope °*Co has been studied by the method of nuclear alignment 
in cobalt rubidium sulphate. By determining the angular distribution of 
the y-radiation following its decay and its plane of polarization it is shown 
that the 123 kev transition is predominantly M1 with a small admixture 
of EK2 (6=-+0-19+0-02). 


§ 1. LyrropucTION 


In an earlier paper (Bleaney et al. 1954) a detailed account was given 
of a method for aligning cobalt nuclei by cooling a single crystal of a 
paramagnetic Tutton salt to temperatures of about 0-01°K. At these 
temperatures the lower h.f.s. levels are preferentially populated and an 
alignment of the nuclei is obtained. The spatial anisotropy of the 
y-radiation emitted in the decay of the aligned radioactive isotopes ®°Co 
(Bleaney e¢ al. 1954) and of *8Co (Daniels et al. 1953) was studied in 
previous experiments which have now been extended to °’Co. 

The general features of the decay scheme of °’Co are shown in fig. 1. 
The nuclear spin and moment of 5’Co which are almost identical with 
those of stable *°Co have been determined by Baker et al. (1953) and 
identify the state as f7/2. The most probable values for the spins and 
parities in the sequence J,, J,, J, are 5/2~, 3/2-, 1/2~, the prominent y-ray 
transition y, being mainly M1 in character (Lemmer et a/. 1955). Although 
the principal features of the decay scheme seemed to be well understood 
we thought it worthwhile to apply the nuclear alignment method to the 
study of this nucleus both in order to confirm these assignments and to 
obtain further information on the decay scheme. 


§ 2. EXPERIMENTAL EQUIPMENT 
Six single crystals of a Tutton salt of the composition (0-01 Co, 0-12 Cu, 
0-87 Zn) Rb,(SO,).6H,0, were grown from a radioactive solution containing 
some ®7Co which had been prepared about two years previously by 18 Mev 


* Communicated by the Authors, 
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deuteron bombardment ofiron. The total *7Co activity in these crystals was 
about 20 uCurie ; the *6Co activity which was the principal contaminant 
had, at the time of these experiments, decayed to a low value. These 
crystals were mounted in parallel crystallographic orientation in the low 
temperature cryostat. Adiabatic reduction of the magnetic field at 
0-9°K from 25 kilogauss to zero caused the temperature of the crystals to 
fall to between 0-01°K and 0-02°x. Scintillation counters were then 
placed round the cryostat and the dependence of the anisotropy («) of 


Fig. 1 
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Decay scheme of Cobalt 57. 


the y-radiation (y,+-y.) was observed as the crystal sample returned to 
0-9°K (e=1—J/(K1)/J(K2) where J(K1), I(K2) are the y-ray intensities 
observed along the magnetic axes K1, K2 of the crystal). The variation 
of « with 1/7’ is shown in fig. 2 where it is seen that even at the lowest 
temperature the anisotropy does not exceed 0-004. That this unexpectedly 
negative result was not due to failure to align the cobalt nuclei was 
proved by the fact that the hard background radiation from the residual 
Co activity did have an appreciable anisotropy. 
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§3. THEORETICAL DEPENDENCE OF ANISOTROPY ON TEMPERATURE 


The theoretical dependence of anisotropy on temperature for the 

decay scheme of fig. 1 has been evaluated using the treatment of Steenberg 
(1952) and is also shown in fig. 2. It is clear that there is a marked 
discrepancy between theory and experiment of at least a factor of 70. 
For the calculation four assumptions have been made : 
_ (i) The anisotropy arises from the y.-radiation only. Using NalI(TI) 
scintillation counters it was not possible to separate the contributions due 
to the 123 kev y-ray (y,) and the weak 137 kev y-ray (y,). From the 
known intensities of these y-rays it can be concluded that the detection of 
y1 Will not cause the expected anisotropy to be reduced by more than 
25°, at any temperature. 

(ii) That the h.t.s. is completely anisotropic. Since the slight departure 
from complete anisotropy of the h.f.s. will alter the expected y-ray 
anisotropy by only about 2°/ (for a discussion see Bleaney et al. 1954) 
this assumption is justifiable. 


Fig. 2 


Dependence of anisotropy («) on 1/7’ and the theroretical curve for a pure 
M1 transition. 

(iii) That the alignment of the state J, is not perturbed during its finite 
lifetime by the influence of extra-nuclear fields. The lifetime of this state 
can be estimated from the yield of Coulomb excited radiation (Temmer and 
Heydenburg 1954) and is found to be 52x 10-*%sec. For lifetimes of this 
order of magnitude perturbations of nuclear states are known to occur. 

(iv) That the y-radiation y, is pure dipole. 

Admixtures both of M2 and E1 and of E2 and M1 radiations have been 
identified in many y-ray transitions and it would not be surprising if such 
a mixture was present. 
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It is therefore clear that the explanation for the discrepancy between 
the theoretically predicted and the observed behaviour lies in incorrect- 
ness of assumptions (iii) or (iv), or both. It is possible to distinguish 
between the two possible explanations by determining the state of polariza- 
tion of the y-radiation from the oriented nuclei. The highest degree of 
polarization will in general be observed in the radiation emitted at right 
angles to the axis of alignment. This degree of polarization can be 
described by the intensities J, and J, of plane polarized components of the 
radiation in the directions parallel and perpendicular to the axis of 
alignment. Ifthe assumption (iii) is incorrect the perturbation will cause 
the ratio p=(I,—I,)/(Z,+1,) to be reduced in direct proportion to the 
reduction of the anisotropy (ec). If, however, it is assumption (iv) which 
is at fault, and it is assumed that there is an admixture of coherent 
quadrupole radiation which makes the y-ray intensity isotropic, it can be 
shown that the ratio » will remain almost unchanged. 


§ 4, POLARIZATION OF THE y-RADIATION 


The polarization of the radiation was measured using the apparatus 
described earlier (Bishop and Perez y Jorba 1955). This was set up so as 
to detect radiation emitted along the magnetic K3 axis of the crystal, for 
in this direction the polarization of the radiation should be greatest. The 
degree of polarization p was found by measuring the number of y-rays 
(NV, and N,) scattered in directions parallel and perpendicular to the 
axis of alignment, the scattering being a maximum in the direction at right 
angles to the plane of polarization. The degree of polarization p can be 
determined from the relationship 


N JN, =[+p)+ RU —p)]/(—p)+ R01 +p)] 
where RF is the analysing ratio for the polarimeter, determined from the 
geometrical arrangement and the differential Compton scattering cross 
section for the y-radiation in question. For 123 kev y-radiation the value 
of F is about 6. 

Figure 3 shows the dependence of N,/N,, on 1/7 and indicates a small 
but significant polarization at the lowest temperatures. Since V,/N, is 
greater than unity the above expression shows that p is negative. Since 
the measurements of Alburger and Grace (1954) indicate that this 
transition is predominantly dipole the sign of the polarization establishes 
that it must be magnetic (see for instance Steenberg 1953) thus confirming 
the suggested assignment of parities to these levels. The degree of 
polarization however is at least 6 times smaller than that to be expected 
at the temperatures reached and it can be concluded that this reduction is 
due to perturbation of the state /,. 


§5. Discussion 


If the transition is both mixed and follows a state which is perturbed 
the angular distribution of radiation intensity can be written, for dipole 
transitions, as: 


I(0,6)=1+-AB,U »(F'2+-892)Po(cos 6)-+ $0AB.U o( F .—8g4/3) cos 24P,2(cos 4) 
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(see Gray and Satchler 1955) where ¢ is the inclination of the electric 
vector of the radiation to the plane containing both its direction of 
emission and the axis of alignment, A is an attenuation coefficient (<1) 
describing the perturbation (this is assumed either to be random or to have 
symmetry about the symmetry axes for the ions), U,, F, and g, are 
_ coefficients involving the angular momenta of the transitions, B, is the 
population function of the initial state and depends on its magnetic 
moment and temperature, P,(cos @), P,?(cos @) are the Legendre and the 
associated Legendre polynomials for cos @; 8, which in the above expression 
is assumed to be small, is the ratio of the amplitudes of quadrupole and 
dipole in the mixed transition and o takes the values —1 and +1 for 
electric dipole and magnetic dipole transitions respectively. It can be 
seen clearly from the above expression why a vanishingly small anisotropy 


Fig. 3 
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y-ray polarimeter measurements. 
Dependence of V,/N,, on 1/7 and the theoretical curve for A=0-11, 6=+-0-19. 


is not necessarily accompanied by vanishing polarization. The reason for 
this is that the anisotropy, since it is measured with polarization insensi- 
tive detectors, is independent of the term containing ¢. This leads to the 
condition for zero anisotropy 5=—F,/g,; the polarization of the 
radiation for this value of 8 still remains finite and can be predicted for 
any particular value of A. Hence A can be determined since B,, U2, Fs 
and g, are known for this isotope. A and 6 may therefore be determined 
directly from the values of the polarization p and the anisotropy « obtained 
at the same temperature. 

Using the values of « and p at a temperature of 1/80°K it is found that 
§ is +0-19-+-0-02 and Ais 0-11. Account has been taken here of the small 
effect due to the detection of y,. In fig. 3 the theoretical curve using 
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these values of \ and 6 is shown as well as the experimental points. The 
agreement at temperatures other than 1/80°K gives one confidence not 
only that this small polarization is significant but also that the interpreta- 
tion of these results is correct. 


$6. CONCLUSIONS 


The residual polarization of the y-radiation and its nearly isotropic 
distribution supports the spin and parity assignments shown in fig. 1 
and requires an E2/M1 admixture with 6=+0-19+0-02. The value of A 
(0-11) indicates that the state 7, is strongly perturbed by extra-nuclear 
fields. These presumably arise as a result of the K-capture process. 
The large line widths of the x-ray lines indicate that the lifetimes of 
K and L shell vacancies are too short for the fields arising from them to 
influence this nuclear state. It is therefore likely that the perturbing 
fields originate in vacancies in the outer electron shells. 
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CLIX. On the Scattering Cross Section of an Obstacle 


By D. 8. Jonzs 
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SUMMARY 

Van de Hulst’s theorem giving the relation between the scattering 
cross section and the scattered ‘ amplitude ’ for an obstacle in an incident 
plane wave is derived by a calculation of the energy flow. In addition 
it is shown that, when the incident wave is due to a point source, the sum 
of the scattering and absorption cross sections is determined by the value 
of a certain field in the neighbourhood of the source. 

The argument, which is based on energy considerations, could be 
applied to many scattering problems, e.g. sound waves or atomic collisions 
but is given explicitly for the electromagnetic field. 


§ 1. IyTRoDUCTION 


Ir has been shown by van de Hulst (1949) that the scattering cross section 
of an obstacle in an incident plane wave is closely related to the 
“amplitude ’ of the wave scattered in the direction of the incident wave. 
(Various special cases of his theorem have been derived by, for example, 
van de Hulst (1946), Levine and Schwinger (1948, 1950) and by Storer 
and Sevick (1954). A proof for atomic collisions has been given by 
Feenberg (1932) and by Lax (1950).) The argument of van de Hulst is 
that the attenuation of the plane wave after it has passed the obstacle is 
a measure of the absorption and scattering of the obstacle. Since the 
scattered field and plane wave do not interfere at infinity except in the 
direction of the incident field the attenuation can be estimated by 
considering the neighbourhood of this direction only. The calculation, 
which involves a difficulty over the convergence of an integral, shows 
that the attenuation depends only on the ‘ amplitude’ of the scattered 
wave in the direction of the incident wave. 

In this paper we shall consider firstly the same problem using a some- 
what different method which avoids the convergence difficulty. The 
method, which is based on energy considerations will be given explicitly 
for the harmonic electromagnetic field, so that polarization effects may 
be observed but, apart from an alteration of the mathematical equations, 
most of the statements could be applied to other scattering problems, 
e.g. those concerned with small-amplitude sound waves or atomic 
particles. 


* Communicated by the Author. 
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Secondly we shall consider the problem which arises when the incident 
wave is due to a point source. In this problem the scattered field may 
be regarded as due to certain sources, situated within the obstacle, 
radiating into free space. Let the field which is produced by replacing 
all these sources by equal sinks (which give ingoing waves at infinity) 
be called the equivalent ingoing field. Then the sum of the scattering 
and absorption cross sections of the obstacle in the point source field is 
equal to average rate at which work is done in maintaining the source 
in the equivalent ingoing field only divided by the average rate at which 
incident energy is crossing unit area. 

Some readers may find this clearer when they see the appropriate 
mathematical expressions. These will now be given for two particular 
examples. 

Firstly, we consider scalar fields and suppose that the incident field, 
due to a source at r’, is derived from 


V'o+h*o=—gs(r—r’) 
where q is a constant. Then the sum of the absorption and scattering 


coefficients is 1672 
a SS is 1\ 7) 
ko 1 py iPalr )a*} 


where ¢, is the equivalent ingoing field, o the solid angle subtended by the 
obstacle at r’ and the star indicates a complex conjugate. 
Secondly, consider an electromagnetic field with the incident wave 


derived from curl Ey = —iwpgHy, 
tk? F 
curl Ho hy E,+po(r—r ) 


where p is a constant vector. Then the sum of the average rates at which 
energy is scattered and absorbed by the obstacle is equal to 


—PARIE,(r’) . p*} 
where E, is the electric intensity in the equivalent ingoing field. 


§ 2. THe INcipENT PLANE WAVE 


Let S be a sphere of radius R which completely encloses the obstacle. 
Then the average rate at which energy is entering S is the average rate 
at which energy is absorbed by the obstacle.t 

Let E, H and E,, Ho be the electric and magnetic intensities in the total 
field and incident field respectively, the time dependence exp (iwt) being 
understood. Denote by .7 the average rate at which energy is absorbed 
by the obstacle. Then the above statement about the rate of energy 
leaving S may be written 


re 


t{ (EAH*+E*,H).ndS=—v . ... . (1) 
Js 

eS ee eee 
} The corresponding statement for the scattering of a beam of particles by a 


Attn - that the number of particles entering S is equal to the number being 
absorbed. 
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where the star indicates a complex conjugate and n is the unit vector 
along the outward normal to S. 

Now the average rate at which the energy of the incident wave only 
enters S is zero. Hence (1) can be written 


| (Ess +Es"4H,) dS =—o—} | (Eyay*+E,aHyt+E,*aH, 
EEE AE nus eo, we... (2) 


where E, and H, are the electric and magnetic intensities in the scattered 
field. The integral on the left of the equation gives the average rate 
at which energy is scattered by the obstacle. The integral on the right 
gives a contribution only from a small portion of S because, over most 
of S, the incident and scattered waves carry energy independently. 

To calculate this contribution we observe that, as R+>«, 


E,—=A,(n) exp (—ikR)/R+0(1/R?), 
H,= = nrA,(n) exp (—ikR)/R+O(1/R?) Sey) 
res) 


J 
where n.A,=0, py is the permeability of free space and the wavelength 
is 27/k. The incident field is given by 


E,—e, exp (—7kRn.n,), 

H,=h, exp (—?kRn . ny) 
where the real constant vectors ey, hy, and no are related by 
@p-ep=1, My. e)=0, my. hy=0 and proohg—kngAxes- ——.-> (4) 
Also there is a lemma (Jones 1952) which states 


z) f(n) exp (—7kRn . ny) dS= ae {f(no) exp (—tkR) 
s 


—f(—np) exp (tkR)}+-0(1) 


as R->co. This lemma enables the integral on the right of (2) to -be 
calculated. After a use of (3) and (4) the limiting value of the integral 
is found to be 


It may be remarked that this is contributed by the portion of S (in the 
direction of the plane wave) which subtends a solid angle of order 
(kR)-°(1>B>%) at the centre of S and thus the attenuation can be 
measured by an instrument which subtends a solid angle of this order. 

The average rate at which energy is scattered by the obstacle is, 
consequently, 
4 : Mitt Fe, -Ay(ny))—L. 

OO 
On dividing by the average rate at which the energy of the incident 
waves crosses a unit ares. perpendicular to its direction of propagation 


2P'2 
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we find that the extinction, i.e. the sum of the scattering and absorption 
cross sections is equal to 


which is the result obtained by van de Hulst. 


§ 3. THE Point Source INcIDENT FIELD 


The incident field E,, Hy produces a scattered field E,, H,, which may 
be regarded as generated by sources, placed within the obstacle, radiating 
into free space. Let E,, H, be the equivalent ingoing field, i.e. the field 
produced when all the sources within the obstacle are replaced by equal 
sinks. As R->oo 

E,.—A,(n) exp (hR)/R+O/R*) 2 | Seas 
where n. A,(n)=0. 

In the calculation of A,(n) the only important quantity, apart from 
the source strength, is the phase difference between the position of < 
source and a plane perpendicular to n. A similar remark is true for 
A,(n) and clearly the phase difference occurring for A,(n) is the same 
as that occurring for A, in the direction —n. Hence 


A,(n\=A,(—n); ©) . ou ee 


ie. at infinity the amplitude of the equivalent ingoing field can be 
obtained by rotating the scattered field at infinity through 180°. One 
consequence of this is that the average rate at which energy is scattered 
by the obstacle is the same as the average rate at which energy is brought 
into S (a sphere which encloses the obstacle and the source of the incident 
field) by the equivalent ingoing field. 

Consider now the field E,—E,, H,—H,; the average rate at which 
energy in this field leaves S’, a surface enclosing the obstacle but not 
the source of the incident wave, is equal to the average rate at which 
work is done in maintaining the sources and sinks inside S’. In fact, 
there are only sinks namely those producing the equivalent ingoing field. 
The average rate at which work is done can be split into two parts 
(i) that done by maintaining the sinks in the field E,, H, and, (ii) that 
done by maintaining them in —E,), —Hpy. The first of these is just the 
average rate at which energy is taken out of S’ (or S) by the equivalent 
ingoing field alone and hence is equal to the average rate at which energy 
enters S’ due to the scattered field alone. This, again, is equal to the 
average rate at which work is done in maintaining the sources of the 
scattered field in —E,,—H,. In view of the relation between the sources 
of the scattered field and the sinks of the equivalent ingoing field the 
combination of this average rate and (ii) gives the average rate at which 
energy enters S’ in the field E)--E,, Ho+Hy,, and this is the average 


_ {The corresponding quantity for small amplitude sound waves, with incident 
field exp (—ikRn . ng) and scattered field 4 , (n) exp(—ikR)/R, is —(47/k) £A,(n,). 


Scattering Cross Section of an Obstacle 961 


rate at which energy is absorbed by the obstacle. Therefore, the average 
rate at which energy leaves S’ in the field E,—E,, H,—H, is equal to 7. 

Now the average rate at which energy leaves the interior of S’ in E,—E, 
is equal to the sum of the average rate at which work is done in maintaining 
the source of the incident wave in the field E,—E, and the average rate at 
which energy leaves S. Since the incident wave and equivalent ingoing 
field do not interfere the average rate at which energy leaves S is equal 
to the difference between the average rates at which energy leaves S due 
to the incident wave and scattered field separately. Hence the sum of 
the average rates at which energy is scattered and absorbed by the obstacle 
as equal to the average rate at which work is done in maintaining the source 
of the incident field in the equivalent ingoing field only. This is equally 
true if the incident wave is due to a number of point sources. 

Denote by o the area of the sphere, with centre the source of the 
incident field and unit radius, which lies within the tangent cone from 
the source to the obstacle. Define p, the average rate at which incident 
energy is crossing unit area, to be the average rate at which the energy 
of the incident field is crossing o divided by o. Then the scattering 
cross section is defined to be the ratio of the average rate at which energy 
is scattered and p. With a similar definition for the absorption cross 
section we find that the extinction is equal to the average rate at which 
work is done in maintaining the source of the incident wave in the presence 
of the equivalent ingoing field only divided by p. 

The cross section, defined in this way, does not become the plane wave 
cross section as the source tends to infinity. A more su table quantity, 
from this point of view, is the scattering coefficient which is defined as 
the ratio of the average rates at which energy is scattered and at which 
incident energy is crossing o. The value of the sum of the absorption 
and scattering coefficients may be deduced easily from the preceding 
paragraphs. 

When the incident wave is obtained from 


curl Ey>=—itwpu Hy, 
LHe Ue E d(r—r’) 
LO ie o+ po(r— 


where p is a constant vector and 6 is the Dirac delta function, the above 
results show that the sum of the average rates at which energy is scattered 
and absorbed by the obstacle is equal to 
—}{E,(r’). p*+E.*(r’) . p}-t 

If the source is some distance from the obstacle this reduces, on account 
of (5) and (6), to 
—i{A,(—n’) . p* exp ((kR')/R’+A,*(—n’) .pexp (—ikR’)/R’} 

BOM eerie oe (7) 
where r’=R’'n’. 


+ The corresponding formula for scalar waves is stated in the introduction. 
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The result of § 2 can be deduced from this. For, if we put n’=—n, 
and p=(47t/wp)) R’ exp (ikR’)ey, the incident wave becomes the plane 
wave of §2 in the limit as R’+00 and (7) becomes —(27/wp1)-%A4(Mo) « &p 
in agreement with the earlier result for the sum of the average rates 
at which energy is scattered and absorbed. a 


§ 4. THe HichH FREQUENCY EXTINCTION 


When the wavelength is very small compared with the linear dimensions 
of the obstacle the limiting case of geometrical optics is approached. 
Only that part of the incident wave which lies within the tangent cone 
to the obstacle is then affected by the obstacle. If there is no trans- 
mission through the obstacle the scattered field inside the tangent cone 
and on the opposite side of the obstacle from the source must annul the 
incident field. The incident field itself is either reflected or absorbed by 
the obstacle. Hence the sum of the average rates at which energy is 
scattered and absorbed by the obstacle is twice the average rate at which 
energy is incident inside the tangent cone. Consequently the extinction 
is twice the geometrical cross section, the geometrical cross section being 
the area of the unit sphere inside the tangent cone. The sum of the 
absorption and scattering coefficients is 2. 

It may be remarked that this simple result would not hold in general 
if p were defined by the incident energy over the whole unit sphere 
because the energy of the source is not necessarily radiated uniformly 
in all directions. 
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ABSTRACT 


A Nal(Tl) crystal was used to study the gamma radiations emitted by 
natural iron, *4Fe, and **Fe scatterers. Gamma rays of the following 
energies were observed: from natural iron 0-85-+0-03, 1-22--0-05, 
1-42+0-05 Mev. The 1-4 Mev gamma ray was not present with the °*Fe 
scatterer. A 1-37-+.0-05 Mev line was observed with the ®4Fe scatterer. 
In addition the neutrons scattered from natural iron were studied with a 
stilbene crystal, and groups corresponding to levels at 0-85--0-03 mev and 
1:45-0:05 Mev were found. On the basis of the above the following 
levels can be identified: in **Fe 0-85+0-03 mev, 2:07-0:05 mev, in 
54He 1-40+0-05 mev. 


$1. INTRODUCTION 


THE gamma ray spectrum associated with the scattering of fast neutrons 
by natural iron has been extensively studied (Beghian et al. 1951, Day 
1954, Lafferty et al. 1954, Rothman et al. 1954, Scherrer et al. 1954). 
Gamma rays of 0:85 Mev, 1-25 Mev, and 1-42 Mev have been reported 
(Day 1954). The 0-85 Mev and 1-25 mev lines fit into the level scheme of 
56He established from the 6 decays of *®Co and **Mn. The 1-4 Mev line 
must either be associated with some level of **Fe which is not excited in 
the B decays, or it must originate from one of the other iron isotopes. The 
present work was carried out to decide between these alternatives, to find 
which levels are excited in the various isotopes and to establish their mode 
of decay. 


§ 2. MEASUREMENTS wiTH NaTuRAL IRON 


A 5cm diameter, 2 cm thick natural iron scatterer was employed and 
studies of the neutron groups and gamma rays made by techniques which 
have previously been described (Eliot et al. 1954). As before, pulse height 
distributions were recorded, on a 25 channel kicksorter, with and without 
scatterer, and using both a sodium iodide and stilbene crystal as detectors. 

The gamma ray spectrum was obtained from the sodium iodide measure- 
ments by subtracting (after suitable normalization) the pulse height 


* Communicated by the Authors. 
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distribution taken without scatterer from that made with the scatterer. 
In this way the effect of background neutrons was eliminated. Gamma 
rays of 0:85-L0-03 Mev, 1-22-+0-05 mev and 1-45+-0-05 Mev were observed 
(see fig. 1). In addition there was a small] peak at 2-1 Mev (see § 5). 

The neutron spectrum was determined from the stilbene pulse height 
distributions as before (Eliot et al. 1954) by plotting R as a function of 
pulse height where 


counts per kicksorter channel with the scatterer ) 
a = per kicksorter channel without the scatterer 
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Natural iron scatterer. Gamma ray spectrum. 


The resulting ratio plot will, in general, consist of a series of plateaux 
separated by steps. Each step can be interpreted as a neutron group or 
gamma ray Compton edge. In the present case the Nal measurements 
indicate that all the gamma rays are too energetic to produce Compton 
edges on the ratio plot. Any steps must therefore correspond to neutron 
groups. 

Figure 2 shows the presence of neutron groups exciting levels at 
0:85+0:03 Mev and 1-45-+-0-05 mev. It has been demonstrated in 
subsidiary experiments (details to be published) that gamma rays of 
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all energies produce a distortion of the plateaux in the region of low 
pulse heights. This makes it impracticable to search for neutron groups 
corresponding to levels of more than 1-7 Mey. 


Fig. 2 
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Natural iron scatterer. Ratio plot. 


‘Fig. 3 
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56He scatterer. Gamma ray spectrum. 


§ 3. MEASUREMENTS WITH *°Fe 
The scatterer was 29-62 of 97%, °®Fe powder enclosed in a thin 
(0-1 mm) aluminium can. Gamma ray measurements only were made and 
it is seen from fig. 3 that the 1-4 Mev gamma ray is absent. The 1-22 Mev 
gamma ray was still present and can thus be definitely attributed to °*Fe. 
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§ 4. MEASUREMENTS with °*Fe 


3 g of 54Fe enriched to 67° enclosed in a thin can of aluminium were 
used as scatterer. Again only the gamma spectrum was observed (see 
fig. 4). Itis seen that a 1-37--0-05 Mev gamma ray was present and within 
experimental error this can be identified with the 1-42-0-05 mev line 
observed with natural iron. 


§ 5. CONCLUSIONS 


It has been shown that the 1-4 Mev gamma ray is associated with °*4Fe 
and the neutron groups from natural iron indicate the presence of an 
excited state of this energy. It therefore is reasonable to interpret the 
1-4 Mev gamma ray as a ground state transition in *4Fe. A similar 
conclusion has recently been reached by Sinclair (1955) who studied the 
gamma rays from a 32 g scatterer enriched to 35% °*Fe. The cross 
section for excitation of the 1-4 Mev level in *4Fe, calculated from the 
ratio plot, is 1-0--0-3 barns. 


Fig. 4 
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4Fe scatterer. Gamma ray spectrum. 


In the measurements with natural iron both the neutron group exciting 
the 0-85 Mev level in °*Fe and the ground state gamma ray were observed. 
The isotopic cross section for this level is 1-06-+-0-25 barns. Since there is 
no neutron group corresponding to a 1-22 mev level, the 1-22 Mev gamma 
ray must be regarded as a cascade transition from a 2-07 mev level into 
the 0-85 Mey level. This is consistent with the data on the 6 decay of 
*°Co (Sakai et al. 1954). The 2-1 Mev peak present in the gamma ray 
spectrum can be attributed to an addition process in the Nal crystal of 
the radiations of 1-25 and 0-85 Mey. This interpretation is consistent 
with the observed intensities but the possibility of a crossover transition 
cannot be excluded. 
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The conclusions are summarized in the level schemes of fig. 5. The 
results obtained with °4Fe show that this technique of measuring gamma 
spectra is applicable to only a few grams of scattering material and is 
therefore particularly useful in the case of separated isotopes. 


Fig. 5 
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The level schemes of *4Fe and °°Fe. 
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SUMMARY 


The stresses round a piled-up group of dislocations may either initiate a 
crack or operate a Frank—Read source by pulling it from its locking 
impurities ; the latter is related to the propagation of a Liiders band. 
This interpretation leads to results in accord with experiment. The model 
gives a sharp transition from brittle to ductile behaviour as the temperature 
is increased. 


§ 1. INTRODUCTION 


THE conditions under which the stresses round a piled-up group of dis- 
locations are sufficient to open up a crack have been considered in a 
previous paper (Stroh 1954); to form a crack we require dislocations 
piled up under a shear stress o, such that 


nbo,= 12,0 (S.A 


where b is the Burgers vector and y the surface energy. This equation can 
be put in a form more suitable for comparison with experiment by using 
the result due to Eshelby, Frank and Nabarro (1951) that the pile-up will 
. occupy a length L of the slip plane given by 


L=nbG/r(l=—v).@,.. si. 4 ee 


where G is the rigidity and v Poisson’s ratio. Eliminating n between (1). 
and (2) we obtain 


o,2L=12yG4/n(1—v). iI jes en ge Re 


If the specimen is tested under tension the maximum shear stress is half 
the tensile stress 7, and, supposing the dislocations to pile up on the grain 
boundaries, we must take L equal to half the grain diameter d. Then 
(3) may be written 

rdV2—4{6yQ/r(1—v) #9; 5 www. AY 


thus the cleavage strength should be inversely proportional to the square 
root of the grain size, a result first deduced by Orowan (1933) on the basis 
of Griffith’s theory. This proportionality has been verified experimentally 
by Petch (1953) for ferrite and Greenwood and Quarrel (1954) for zinc. 
The constant on the right hand side of (4) may also be calculated. For 


* Communicated by the Author. 
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ferrite with G=8 x 10 dynes cm~2, y=1-2 x 103 dynes cm~! and v=0-3 
we obtain 7d’2—2-0 108, which may be compared with the value 
2-3 10° from Petch’s graph. For zinc, we have G=3-8x 101, y=103 
and v=0-2, giving 7td¥/2=1-2 x 108, while Greenwood and Quarrel obtain 
(at liquid air temperature) 1-0 x 108. 

This very satisfactory agreement encourages us to consider how far 
these ideas may be applied to other problems. In the present paper we 
discuss the lower yield point of ferrite, and the brittle—ductile transition. 


§ 2. THE Lower YieLtp Pornt 


Petch (1953) has emphasized the great similarity in the behaviour of the 
cleavage strength and the lower yield point of ferrite; the latter also 
varies inversely as the square root of the grain size. If we interpret 
cleavage as resulting from the normal stresses round the pile-up, it is 
natural to suppose that yielding and the propagation of a Liiders band 
arise from the shear stresses round the pile-up. Now at a distance r from 
the head. of the pile-up the shear stress, apart from a factor approximately 
unity but depending on orientation, is 


oo, (L/r)\/. PE cams ey ie Sant (0) 


If slip is to continue, this stress must be great enough to operate the 
sources in the next grain. These sources will be locked by impurity atoms 
(Cottrell 1948, Cottrell and Bilby 1949) and to free the dislocation a stress 
a1, equal to the upper yield point is required. In addition the stress on the 
source must be great enough to bow it out ; if the source has length / the 
stress needed for this is G'b/J. Hence the source will operate if 

o,(L/r)}/?=max [o,, Gb/2], a Rel nc fe a's eer) 
where ‘ max ’ indicates that the greater of the two quantities in brackets is 
to be taken. | 

We suppose, now, that there are a number of sources of varying lengths. 
Those of length / cannot be nearer than a distance / from the pile-up ; 
in general the nearest sources will be a distance r=al away, where « is a 
constant of order, but greater than, unity. With this value ofr, (6) may be 
written 


Glo) == Max jo,)/ 2G bli}. Ca ee Came or FT) 
The right hand side of (7) is least when the two terms are equal, or 
Le Ge Pie oa oo (8) 


and these sources will be the first to operate. With o, equal to the 
observed upper yield point of about 7 x 10° dynes/em (at room temperature), 
this gives =3x10-5 cm. Substituting (8) in (7) we have 
o,(L/«)? = (Gbo,)™, 
or taking L half the grain size d, and o, half the tensile stress r, 
| 7dV2=2(2a0,Gb)¥2—=constant. . . . . . (9) 
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The constancy of 7d1/2 is the relation between lower yield point and grain 
size observed by Hall (1951) and Petch (1953). 
At two different temperatures 7’, and 7’, eqn. (9) gives 


(rd2)_|(rd¥2)p=(0,6)7[(0,G)2*, - «+ (10) 


where the suffixes refer to the temperatures. 

To test this equation we may use the results of Hall at room temperature 
and of Petch at liquid air temperature. From the graphs in their papers 
we find that 7d1/2 has the values 0-7 x 108 at room temperature and 1-9 x 108 
at —196°c, the ratio being 2-7. Both experiment (McAdam and Mebs 
1943) and theory (Cottrell and Bilby 1949) agree in giving a change in the 
upper yield point o, by a factor of 4:1 between the temperatures under 
consideration. The change in G@ will be about 10°. Then the right hand 
side of (10) is equal to 2-1, which is of the right order. It is possible that 
the agreement could be improved by taking into account the rate of 
straining. Hall used the unusually slow rate of 10-4/min ; Petch does not 
state his strain-rate but if it were greater, the effect would be to accentuate 
the difference in temperature. . 

Finally we may consider the value of the constant in eqn. (9). Using 
the value o,=7x 108 dynes cm~ appropriate to room temperature, we 
obtain 3:3 10’ «1/2; this agrees with Hall’s results provided « is about 
4 which seems reasonable. Further the value of 7d’? should be the same 
for all specimens (at one temperature) if o, is the same ; however in all 
but the very purest specimens of ferrite there will be sufficient carbon 
atoms present to.form a saturated atmosphere round the dislocations so 
that the locking forces will have a constant maximum value. In fact 
Petch has found that when the lower yield point is plotted against grain 
size, even for such different materials as mild steel, ingot iron, and spectro- 
graphic iron, the points all lie on a single line with a remarkably small 
amount of scatter. The point is not a trivial one, as if the value of the 
yield point depended on some quantity such as the mean size of the 
dislocation network considerable variation might occur. 

Thus the model proposed accounts for the variation of the lower yield 
point with grain size and temperature, its insensitivity to the nature of the 
specimen, and, in order of magnitude, for the absolute value of the yield 
point. 

§ 3. BrirrLe—DuctTimLE TRANSITION 


We have seen that the stresses round a pile-up may either initiate a 
crack or operate nearby sources; if the former happens the metal is 
brittle, if the latter, ductile. When the sources are locked by impurity 
atoms, thermal fluctuations will help free them, and so high temperatures 
will favour ductile behaviour. At low temperatures no thermal fluctua- 
tion great enough to liberate a source may occur in the time taken for the 
pile-up to grow large enough to initiate a crack ; the material is then brittle. 
In the present section we attempt to estimate the temperature at which 
the transition from brittleness to ductility will take place. 
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Let o be the stress at time ¢ on a source of length / near the pile-up. 
The probability that the source will operate in the next interval of time 
dt, provided it has not done so already, is . 


y(I/b) exp {—u(a)/kT} dt, 


where v is an atomic frequency of vibration, and wu(c) is the activation 
energy for the source to break away from its locking impurities ; J/b is 
the number of points at which the break away can start. If p is the 
probability that the source has not operated at time ¢, the probability 
that it will do so between times ¢ and ¢-+-dt is 


—dp=pv(l/b) exp {—u(o)/kT} dt. . . . . . . (11) 
On integrating (11) we have 
1 yt 
p=exp | - _ exp {— eal ar| Ar nip res Eo) 


We require the probability that the source will not have operated by the 
time the pile-up has grown large enough to initiate a crack ; (12) is now 
taken to refer to this time. The contributions to the integral will come 
mainly from times at which u(c) is least, which will be when o is greatest. 
The integral will be approximately ¢, exp {—u(c)/kT7'} where f, is a time of 
the order of that required to form the pile-up, and o now denotes the final 
value of the stress due to the pile-up on the source. Then (12) may be 
written 

p=exp [—(vli,/b) exp {—u(o)/kT}].  . 2. . . . (18) 


Because of the double exponential p will change fairly abruptly from 1 to 0 
as 7 increases ; this is in accord with the sharp transition from brittle to 
ductile behaviour observed. The transition will occur at a temperature 


u(c) 


eat log (vlt,/b) ’ 


(14) 
To estimate u(c) we require the stress on the sources : from (3) and (5) 

this will be 
o= {12yG/n(1—v)r}”. Rie cre eee nS aig 1D) 


The considerations of § 2 suggest that the first sources to be operated will 
be about 1-2 10-4 cm from the pile-up. With this value of r (15) gives 
o=101° dynes/em?. Taking the yield stress at absolute zero to be 3x zo 
dynes/cm?, curve 4 of Cottrell and Bilby (1949), which they consider to 
be their best, gives u to be about 0-6 ev. Other quantities, which appear 
only in the argument of a logarithm in (14) are not very critical ; taking 
v~1012 sec-1, ty~10-* sec, and 1/b~10* we obtain 7=300°K. This is of 
the right order, though there is clearly considerable uncertainty as to the 
exact value. Under the same conditions we find from (13) that p decreases 
from 0-9 to 0-1 for a rise in temperature of 40°c ; this may be taken as a 
measure of the width of the transition range. 


972 On Brittle Fracture and Yielding 


We may be interested in the spread of a crack instead of in its initiation. 
In Robertson’s experiments (1951) a crack started from a saw cut is 
propagated through a temperature gradient, and the temperature where it 
is arrested observed. The moving crack will produce stresses on any 
sources it passes ; these may or may not then break away from their locking 
impurities and generate dislocation loops. As before we can find a transi- 
tion temperature ; this will be given by (14) with ¢ now the time taken 
for the crack to pass the source, i.e. ty~r/v, where r is the distance of 
closest approach and v the speed of the crack. This temperature will be of 
similar magnitude, but not identical with that obtained previously. 
Below the transition temperature the crack can propagate without 
disturbing any source and so the material will certainly be brittle. Above 
the transition temperature the propagation of the crack will be accompanied 
by plastic flow ; if this absorbs too much of the stored elastic energy the 
crack will be unable to spread. It seems reasonable then to expect the 
arrest temperature of the crack will be closely related to our transition 
temperature. 
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ABSTRACT 


The spin-orbit coupling experienced by a nucleon in a nucleus is inter- 
preted in terms of a neutral two-body spin-orbit potential. An expression 
of the Thomas form is obtained for the effective spin-orbit potential. 
It is shown that the spin-orbit splitting in heavy nuclei can be explained 
satisfactorily if the nuclear density distribution has a finite transition 
region at the surface and if the neutron distribution has a larger radius 
than the proton distribution. 


§ 1. INTRODUCTION 


THE success of the nuclear shell model introduced by Mayer (1949) 
and Haxel, Jensen and Suess (1949) implies that a nucleon moving 
in a nucleus experiences a strong spin-orbit force. Recent experiments 
on the polarization of nucleons scattered by nuclei have served to confirm 
the existence of such a force and it is therefore of interest to investigate 
its origin. 

To explain the force which manifests itself in the shell model two 
types of approach have been considered. Firstly, the spin—orbit potential 
has been interpreted in terms of a sort of Thomas (1926) precession, the 
potential having the form h?/2m?c?r dV/dr L.S where m is the nucleon 
mass. The actual relativistic Thomas precession gives much too small 
an effect, but suggestions (e.g. Inglis 1953) have been made that spin— 
orbit potentials of the above form should arise from meson theories. 
Secondly, the spin-orbit potential has been interpreted in terms of two- 
body internucleon forces. In particular, Keilson (1951) considered the 
contribution of tensor forces to the spin-orbit coupling and showed that, 
subject to the validity of certain approximations, one obtains an effective 
spin-orbit potential having the form —Kp(r)/r dp/dr L.S resulting from 
the interaction of an odd particle with a spin saturated core (here p is 
the nuclear density and K a positive number). Splitting of the right sign 
is obtained if the odd particle spends most of its time in a region of radially 
increasing core density. However, particles of high angular momentum 
which are known to experience the largest splitting spend most of their 
time at the ‘surface ’ of the nucleus where the core density is radially 
decreasing. Hence it seems unlikely that the tensor interaction is 


* Communicated by Professor R. E. Peierls, F.R.S. 
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responsible for the strong spin-orbit coupling. On the other hand, a two- 
body spin-orbit force of the form 
V=J([ry—re |) [(ra—re) A (Pi— Pa) - (61+ F2)] 2(T 1,72) 

has been shown by various authors to give spin-orbit splittings of the 
right order of magnitude and the correct sign (Hughes and Le Couteur 
1950, Blanchard and Avery 1951) provided that the isotopic spin factor 
T (1,75) is suitably chosen. These authors, however, do not attempt to 
make detailed quantitative comparisons of theory with experiment or to 
obtain an explicit form for the effective spin-orbit potential. Further, it 
has been shown by Elliott and Lane (1954) that such an internucleon 
interaction can also account for the observed increase in splitting in light 
nuclei as a given shell or sub-shell is filled (e.g. the 1ps,.—1 pj). splitting 
in }N is about twice the 1p,/,.—1p, 2 splitting in °He). 

The object of the present paper is two-fold. Firstly, starting from the 
above two-body force an expression is derived for the effective spin—orbit 
potential experienced by a nucleon in a nucleus. Secondly, it is shown 
that this spin-orbit potential can then account in detail for the sign and 
magnitude of the splitting observed in nuclei in the vicinity of ?°&Pb 
and is also consistent with the 1p;;.—1p,/, splitting in °He. 


§2. DEDUCTION OF THE Sprin—ORBIT POTENTIAL 

We consider here the spin-orbit potential experienced by an odd nucleon 
interacting with a spin saturated core. Because of the necessary anti- 
symmetrization of the wave functions when the odd nucleon is of the same 
type as those forming the core there will be both direct and exchange 
terms to consider in the interaction. Calculations of Blanchard and 
Avery (1951), however, indicate that the exchange terms are only 
appreciable in the case of light nuclei and for this reason they are neglected 
in the following work. 

For the two-body spin-orbit interaction we take the neutral form 

Vig=J (Lei — yl) e— ry) a (Pi— Py) - (0;-+ 9,)] 

since this accounts satisfactorily for ee increase in splitting as a given 
shell is filled (Eliott and Lane 1954). 

The potential experienced by an odd nucleon at r; in the presence of 
a spin saturated closed shell (specified by the usual quantum numbers 
n, 1) is given by : 

Vj=2 | $*VbdV do 


where ¢ is the wave function for the closed shell. 
Since we are dealing with a spin saturated core the expectation values 
of Xp; and 20; will vanish so that 


Vj;=-2{ $*I ([r;—r;|)(r;—r;) a p;.6;dV do 


=~ [22] eared —eb leeward |av,-< 


and the Density Distribution in Heavy Nuclei 975 


where ¢7i(r) is the wave function for a single particle in the state 
specified by the quantum numbers n, /, m (m=magnetic quantum number). 
By symmetry 
28 | Wi (ree eileen) dV, 


will have the form F'(r;)r; where 
a ath ro; Ee 
Por —22 | deter —es)) [51 | ote, 
m J 


Thus, the effective spin-orbit potential is 
V=Fi(r)rap.co 
=Fir)L.o 


where L and o are the usual orbital angular momentum and Pauli spin 
operators respectively, and the suffix 7 has been omitted. 


§3. EvaLuaTiIon oF F(r) 


In order to obtain an explicit expression for F(r) we use a method of 
evaluation valid for short range forces used by Brink (1954). 
Denoting 


25h" (r 
m 


2 —1 Yni(t’) by —p(r, r’) 
we have F(r)=— | I(|r—r'|)p(r, dV" 


Writing p(r, r’)=p’(s, r) where s=r’—r and expanding p(s, r) about 
s=0 gives 


F(r)=—4r | J(s)s?ds [0'(0, I— F5 | Tls)stds [V2p"(s, rho +- 


The first term in the expansion of (7) vanishes since, for s=0, r=r’ and 
oer 
v2 
and one obtains for F(r) to this approximation 


Poy — FED (| seostas) = FRO 


where f,,,(7) is the radial function for a single particle in the state (nl). 
But if the mean density distribution of particles at a radius r in the closed 
shell is p,,(7), then 


—1=0. The evaluation of the second term presents no difficulty 


py lr) =2f2,(r) (214+ L)/4a 
4 \ bdo. 
so that er J(s stds) = rae 


If we now sum over all particles in the nucleus, it follows that we 


can write F(r) in the form 
F(r)=K/r dp(r)/dr 


where 5 K=—47/3 | a s)s* ds 


22 
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is a positive quantity (assuming J(s) to be of the usual form for nuclear 
forces) and p(7) is the nuclear density distribution. 

Thus it turns out that the radial dependence of the resulting spin-orbit 
potential is of the same form as the Thomas type except that the nuclear 
density distribution replaces the potential V(r) in which the nucleon is 
moving. It is to be expected, of course, that p(r) and V(r) will have a 
similar sort of behaviour. Although one cannot be sure of the effective 
nuclear potential one has a fairly good idea as to the nuclear density 
distribution (e.g. from electron scattering experiments) so that the above 
form of F(r) gives a reasonable guide to the radial dependence of the 
spin-orbit potential which might be used for calculations in scattering 
problems. It is interesting to note that the above result is independent 
of the form of the radial functions used. 


§4. CALCULATION OF THE NUCLEAR SPIN—ORBIT SPLITTING 


From the form of F(7) it is clear that shell model spin-orbit splitting of 
the right sign is obtained in regions of radially decreasing density, i.e. 
at the surface of the nucleus. Calculation of the splitting to be expected, 
however, must necessarily be fairly approximate since some assumption 
has to be made about the radial function of the odd particle. For the 
purpose of the calculation the following crude assumptions are made : 

(i) The nuclear density distribution p(7, t) is taken to have the form 


p(r, t)=po; 0<r<R(1—t) 
p(r, =p R—r/Rt,  RO—t)<r<R 


where ¢(0<t<1) is a parameter measuring the size of the transition region 
in which the nuclear density decreases from p, to 0. 

(ii) The radial function is taken to be that of a particle moving in an 
infinite well of radius R. 

These two assumptions are consistent in that both the density and the 
radial function vanish for values of 7 greater than R. The radius R, 
of the nucleus is conventionally defined as the radius of the uniform 
density distribution which has the same mean square radius as that of the 
distribution under consideration. This implies 


RB) RP 
where 
f2(t)=3(4—6t-+ 4t2—13)/2(6— 15t+ 2002— 1503+ 6t4— #5) 


and for a nucleus of atomic number A we take Ry=r,4”3, Furthermore, 
the normalization of p(r) must be such that 
“R 
4a | p(t)rdr=A. 
0 


This gives at once 
po= 34/7 R3(4—6t+ 40223), 
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The spin-orbit splitting 42= By-1-1/2— Ej ~141/2 for an odd nucleon in 
the state (n, /) is then given by 


AK= —(2/-+1) [i | fF2(7) )PR(r yr? dr 


where f,,,(”) is the appropriate radial eanetign for a particle in an infinite 
well of radius R. 


ee OK N= Kirdpjdr—0 for  0<r<R(1—2) 
z ==K oo htr for RI—t)<r<k 
R 
so that AE=(21+1)K po Rt | | falr)Pr dr 
~ RA —t) 


which, for an infinite well can be written 
A4h=y, i(t) CAS 
where 
Ynailt)= [8(27++ 1 Ne (4) [nil F 143 /2(%n 1) ?4(4— 6t-+ 4t?— #8) 6%(t)] 


1 
ion ee 4 
C= 2 | J (s)s* ds 


5 


(nl 


and NA Je ria) dx. 


%pj(1—t) 
Here «,,, is the nth zero of the Bessel function J ,,/9(z). 
The values of y,,,(¢) for t=0-1 to 0-7 for the different shell model 
states are given in table 1. 


Table 1. Values of y,,(¢) 


State GOs t=0-2 t=0°3 t=0-4 t=0:5 t=0-6 t=0°7 


Ip 0-450 1-489 2-857 4-383 5:758 6-896 7-602 


Id 1-251 3-630 6-716 9-341 10-961 11-790 11-924 
lf 2-056 O91. 11-989 15-171 16-219 16-217 15-787 
2p 1-040 3037 3°989 3504 3°103 4-395 6-963 
lg 4-564 11-527 18-608 21-491 21-381 20-396 19-537 


2d 2-358 5-931 6-600 5-239 5-901 9-495 13-053 
lh 5:737 17-024 25-278 27-279 25-696 23-910 22-731 


2f 4-228 9-409 8-880 7-269 10:377 15-675 18-515 
i 8-400 23-629 32-872 33-296 30-105 27-644 26-207 
3p 1-801 3°657 2-914 Pasian 4-228 4-232 5-141 
2g 6-659 12-941 10-544 10-187 16-136 21-640 23-121 


3d 3°632 6-193 4-533 5-841 6-937 6-788 10-933 


Since U is positive and y,,(¢) is also positive for all values of ¢ it is 

clear that splitting of the sign required for the shell model is obtained. 

We now look for values of U and t which give the correct spin-orbit 
splitting in the region of the double closed shell nucleus ?°*Pb. 

Stripping experiments (Harvey 1953) indicate that the 2g5/2.—287/. 
and 3d,,.—3d,,, splittings for the odd neutron in ?Pb are 2-03 Mev and 
0-98 mev respectively. Further it is possible to make an estimate of the 
neutron 1i,3/2—1i,4/2 splitting (see figure). 
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Hitherto this splitting has been estimated by reference to the semi- 
empirical formula (e.g. Harvey 1951). A more accurate estimate is 
given in the following. At 2°8Pb the break in neutron binding energy is 
3-41 mev (Wapstra 1953). Apart from the pairing energy P(3pj/2) 
of the last neutrons in 2°8Pb this can be regarded as the separation in 


MeV 


Single particle levels near the closed shell of 126 neutrons. On the left are 
shown the levels obtained with the potential well referred to in the text. 
(The 2g5,. level is arbitrarily taken as the zero of energy.) 


energy of the single particle levels 3p,;. and 2g5/,. Now in 2°Pb the 
li41/9 level has excitation energy 0-75 Mev (Harvey 1953) and in 2°7Pb the 
liz3/2 level (obtained by excitation of a neutron from lij3;, to 3p4/9) 
has energy 1:63 Mev (Stelson and Campbell 1955). The 1li,3;.—38p4/. 
single particle level separation is therefore 1-63—P(1i,3/.)+P(3p,/9). 
Thus the lij3;2—1i,,/, separation is 


{0-75+ [3-41 —P(3p,/.)]+[1:63—F (1ij/.)+-P(3p4/2)]} Mev. 
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But the difference in binding energy between Pt and 1°7Au gives a 
measure of P(1i,3,.). According to van Patter and Whaling (1954) this 
difference is about 2 Mev. Thus the li,3,.—1i,,/> splitting can be taken 
to be of the order of 3-8 Mev. Further support for this figure comes from 
**’Pb where the 1i,,/, level has excitation energy 3-60 mev (Stelson and 
Campbell 1955). This level results from the excitation of a 3p,/, particle 
to the 1i,;/ state. Using the figure quoted above for the energy of the 
i43/2 state gives for the li,,,.—1i,,/. splitting 


3°60+ 1-63—P(1i43/2)+P(3p 4/9). 
Now P(3p,/,) can be estimated from the difference in neutron binding 
energies in ?°7Pb and 2°*Pb which, according to Wapstra (1953) is 0-65 Mev. 
Taking P(1i,3/.)=2 Mev then gives the 1li,;/.—1i,,/. splitting as 3-88 Mev. 
Good agreement between theory and experiment is only obtained if 
we take U=5-5 mev and t=0-2 (see table 2). 


Table 2. Experimental and Theoretical Values of the Spin—Orbit 
Splitting Near ?°%Pb 


State Experiment Theory 
2g 2-03 Mev 2-03 Mev 
3d 0-98 Mev 0-97 Mev 
hi 378 Mev 3°71 Mev 


| | 


The postulated single particle levels in the region of ?°°Pb are shown in 
the figure. AJso shown are the single particle levels from which the spin— 
orbit doublets result. These levels have been taken at a distance 
J/(21+-1) x (splitting) from the lower component of the spin-orbit doublet, 
this being the correct position when the spin-orbit force has the form 
Fi(r)L.o. These levels are also compared (as far as spacing is concerned) 
with those obtained in a spherical well having the form V(r)=—V, for 
r<R, V(r)=— V, exp [—ep(r—F)] for r>R# with V)>=34-1 Mev and «—1-83 
with »=1-4x10-8 em. ‘Such a well has been deduced by Malenka 
(1952) from a non-linear meson theory and should be fairly realistic in 
that it lies in form between a square well and an isotropic oscillator well. 
The agreement is very good and indicates that a well of this form is 
reasonable for heavy nuclei. It is also interesting to note that such a 
well gives a nuclear density distribution very similar to the one assumed 
here and with 0:3>t>0-2. 


§5. Discusston 


(i) The Value of U 
To see if the value obtained for U is consistent with other experimental 
data on a two-body spin-orbit force we consider the particular case of 
J(s)=SV, exp (—ps*). This gives U=(x/p*)”°3S Vo/ 219°. Recently Hoch- 
berg et al. (1955) have fitted the data on the scattering of neutrons by 
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alpha-particles by using a two-body spin-orbit force of the type con- 
sidered here. They take V,;=22-5 Mev and p=0-2657 x 10° cm~? so 
that U=648 mev for 77>=1:-45x10-% cm. Thus, in order that this 
value shall agree with experiment we must take S=0-09. The above 
authors find that S=0-1 fits the scattering data on *He. The close agree- 
ment between these two values of S is satisfactory considering the 
approximations made in the calculations. 


(ii) The Nuclear Density Distribution 


It is interesting that the spin-orbit coupling data for heavy nuclei 
can only be fitted if ¢, which measures the transition region in the nuclear 
density distribution, is taken as 0-2. That there should be such a transi- 
tion region is indicated for example by the analysis of high energy electron 
scattering by nuclei (e.g. Brown and Elton 1954). These experiments 
suggest a rather larger value for t(~1/3). It is highly probable, however, 
that if a more realistic smoothed out density distribution had been used 
in the present case a larger value for ¢ would have been obtained since 
the smoothing out at the two ends of the transition region will have small 
dp/dr and will therefore only give a small contribution to the spin-orbit 
coupling. In view of the lack of knowledge as to the form of the radial 
functions it was not thought worthwhile to use such a smoothed out 
distribution. 

In the present calculations 7, has been given the conventional value 
1-45 x 10-18 em which is larger than the value (1-2 x 10- em) obtained 
from electron scattering experiments. This latter value, however, 
refers only to the proton distribution and it is to be expected (e.g. 
Swiatecki 1955) that the neutron distribution extends beyond the proton 
distribution by 10-20% and since p refers to the total density it is the 
larger radius that is relevant. In this connection it is instructive to 
compare neutron and proton spin-orbit splitting. Unfortunately, in the 
region of 2°5Pb such a comparison can only be made for * holes’. Thus 
in 2°7Pb the 3p3/2—3p,/. neutron hole splitting is 0-87 Mev and in 2T] the 
2d5/2—2ds/. proton hole splitting is 0-4 Mev. The values predicted on the 
present model for single particles are 0-57 Mev and 0-94 Mev respectively. 
As pointed out by Elliott and Lane (1954) the hole splitting is expected 
to be larger than the particle splitting so that in the case of the neutron 
hole experiment and theory are consistent. This is not so for the proton 
hole. If it is remembered, however, that the proton radius is smaller 
than the neutron radius, it follows that the proton will spend less of its 
time in a region of decreasing core density. Thus if ry (proton)=0-97, 
(neutron) then for the proton we should take t~0-1 rather than 0-2. 
In this case the predicted splitting is 0-38 Mev. Since the hole splitting is 
expected to be about twice the particle splitting this value is probably 
still too large and suggests that 7) (proton)<0-97, (neutron), i.e. 
9 (proton)<1-3x10-% cm in agreement with the electron scattering 
experiments. 
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To obtain a more satisfactory comparison of neutron and proton radii 
it would be useful to compare the splittings of an odd neutron and an 
odd proton outside the doubly closed 2°°Pb shell. This would require 
information about the levels of ?°°Bi which is not at present available. 
On the present model for ¢=0-1 the predicted 2f,,.—2f,,. proton splitting 
in 2°°Bi is 0-66 Mev. 


(iii) Splitteng of the 1h and 1g States 

It is interesting to find the order of magnitude of the 1h,,;.—1hg/» 
and 1g5/.—1g,/,. splittings predicted by this model for the neutron magic 
numbers N=82 and N=50. For t=0-2 the values obtained are 3-5 Mev 
and 3-0 Mev respectively. Work by Harvey (1951) indicates that at the 
above magic numbers the neutron binding energy drops by about 2 Mev. 
However, the spin-orbit splitting should be larger than this for similar 
reasons to those advanced in the case of the 1li,3;.—1i,,/. splitting (e.g. 
the 51st and 83rd neutrons do not go into the 1g,,. and lhg,,. states but 
into 2d,,, and 2f,,,, the 1g,,. and lhg,, states lying higher than these). 
It therefore seems as if the splittings are of the right order of magnitude. 
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[Received May 11, 1955] 


ABSTRACT 
The y-ray spectrum at the 825 kev resonance in #4Mg has been investi- 
vated and the angular distributions of the three initial y-rays measured. 
The spins and parities of the ground state and first excited states are 
confirmed as 5/2+ and 1/2+ respectively and the second excited state is 
shown to be 3/2+. Evidence is presented which lends a measure of 
support to the view that a level exists in ?°Al at 1-92 Mev. 


$1. [INTRODUCTION 

THE mirror pair 2>Al-?>Mg are known to have the same parity and 
spins for their ground states (5/2+) and first excited states (1/2+). The 
second excited states are both either 3/2+ or 5/2*+ (Goldberg 1953), but 
which of these values is correct has not been determined so far in either 
case. The levels above 3-09 Mev in ?°Al have been studied by the elastic 
scattering of protons from *4Mg and their properties are well established 
(Koester 1952). It should be possible to obtain further information 
about the levels of ?°Al lying below the 3-09 Mev level by studying the 
radiation from one of these higher levels. 


§2. EXPERIMENTAL METHOD 

Two resonances occur in the working range of the Liverpool University 
H.T. Set, at 418 kev and 825 kev. The resonance at 418 kev has a y-ray 
yield which is below the useful working limit, so all the work was carried 
out on the 825 kev resonance level, the spin and parity of which are known 
to be 3/2~ from the elastic scattering experiments. 

A 60 u-amp resolved beam of protons bombarded a target of electro- 
magnetically separated 24Mg which had been evaporated on to a 0-02 in. 
copper backing. An insulated target holder enabled the beam current 
to be measured and integrated. The y-rays were detected in a 2 in. 
diameter by 2 in. long thallium activated sodium iodide crystal mounted 
on a Dumont 6292 photomultiplier. The resultant pulses, after ampli- 
fication and delay-line shaping, were analysed by a Hutchinson—Scarratt 
type pulse analyser. Another scintillation counter at 120° to the beam 
was used to monitor the yield and the pulses from this crystal were fed 
to a single channel pulse analyser set to accept all pulses corresponding 
to energy losses in the crystal between 1-5 Mev and 3-9 Mev. 


*Communicated by H. W. B. Skinner. 
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$3. Tue y-Ray SPECTRUM 


In order to be certain that all the y-rays present were correctly identified 
and to enable a comparison of the relative intensities to be made, the 
pulse height distribution produced by a single y-ray was determined for 
y-rays in the range 669 kev to 8-06 Mev, with the aid of various radio- 
active sources and the proton-y reactions occurring in “B, 20, 130 
and }°F. By interpolating between these pulse height distributions, the 
distributions produced by a y-ray of any energy could be determined. 
The interpolation was aided by constructing curves of the quantity 
A/PH versus energy where A is the total area under the curve, P the 
pulse height of the peak and H the number of counts per channel at the 
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peak, in the manner described by Flack et al. (1954). This method was 
also used for estimating the relative intensities of the y-rays and. has 
given consistently good results in this and other work in progress. it 

Knowing the shapes of the spectral distribution produced, the experi- 
mentally obtained spectrum was then analysed by starting at the high 
energy end and successively subtracting the contributions due to the 
various y-rays in descending order of energy. An example of this is 
shown in fig. 1. 


Counts per channel 
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The spectrum obtained is shown in figs. 1 and 2. Analysis of this 
spectrum shows that y-rays are present as given in table 1, together 
with their relative intensities corrected for their angular distributions. 


Fig. 2 
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y-ray energy 
Table 1 
y-ray Energy Intensity (Arbitrary Units) 
Mev 7 
il. 3:04 40-03 274-2 
De 2-59 +0-03 159+10 
3. 2-11+0-03 19+2 
4. 1-:12+0-02 2-+0:5 
d. 0-96 +0-02 11-1 
6. 0-460 +0-01 158+10 
ie 0-510 -+0-01 200 +15 


The suggested cascade system is shown in fig. 3, where the transitions 
are numbered according to their energy as in table 1. The 0-51 Mev y-ray 
in table 1 is due to annihilation radiation from the positron decay of the 
ground state. The relative intensities of the 2-11 Mev y-ray and the 
0-96 Mev y-ray, which are in cascade, show that the 0-96 mev level decays 
partially to the ground state and partially to the first excited state at 
0:46 Mev. The weak 0-50 Mev y-ray produced by this latter transition is, 


of course, lost underneath the combination of the intense 0-51 Mev and 
0:46 Mev y-rays. 
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The weak 1-12 Mev y-ray is thought to be due to a transition between 
the resonance level at 3-04 Mev and a level at 1-92 mev. Goldberg found 
a very weak neutron group from the *4Mg(d, n)?5Al reaction corresponding 
to a level at 1:94 Mev. He also found a neutron group corresponding to 
a level at 1-81 Mev, but a careful search of the spectrum showed no 
indications of a y-ray transition to this state. De-excitation of the 


Fig. 3 


3.04 
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1-92 mev level could occur to the ground state or either of the first two 
excited states. All the spectra taken at 6 different angles showed an 
inflection at 1-45 Mev indicating that the decay probably proceeds at 
least partially to the 0-46 Mev state. The possibility that the 1-12 Mev 
y-ray might be due to a contaminant is unlikely as the most likely 
contamination would be from #*Mg and ?*Mg. The former has a 
resonance at 820 kev but does not give a y-ray of this energy and ?°Mg 
has no resonance sufficiently close to this energy region. 


$4. ANGULAR DISTRIBUTIONS 


For the angular distribution measurements the pulse amplitude analyser 
was set to include pulse heights corresponding to energy losses in the 
crystal from 1-7 Mev to 3-6 Mev. Measurements were made on either 
side of zero to check the symmetry of the system. Spectra were taken at 
6 angles between 0° and 90° in a random order to avoid errors due to 
slow drifts and deterioration of the target. The spectra were analysed 


986 L. L. Green, J. J. Singh and J. C. Willmott’ on the 


as described above, the total area under the total capture peak of each of 
the three initial y-rays was measured, and the results were then corrected 
for the varying absorption in the target backing at different angles. 

The results obtained are shown in figs. 4, 5 and 6, and in table 2 the 
best fits to the points are compared with the theoretical formula, modified 
by the finite acceptance angle (10°) of the scintillator, for electric 
dipole transitions. 


Fig. 4 
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The results confirm the assignments of 5/2+ to the ground state and 
1/2* to the 0-46 mev state and show that the 0-96 mev state is 3/2+. 
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Table 2 
y-ray Energy Transition Theoretical Experimental 
3-04 Mev 3/2-5/2+ 1—0-14 cos? 6 =1—(0-16 40-025) cos? ¢ 
2-59 Mev 3/2->1/2+ 1—0-58 cos? @ 1 —(0-55 40-03) cos? 6 
oa f 3/2-+5/2+ 1—0-14 cos? 6 
2-11 Mev + (0-72 +0:07) cos? @ 
| 3/2--+3/2+ 1+-0-73 cos? 6 


§5. RapiaTiIon WIDTH 


By knowing the solid angle subtended by the crystal at the source, its 
efficiency and the integrated beam current, it is possible to calculate the 
number of y-ray quanta produced per proton and hence the radiation 


Fig. 6 
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width I, of the level at 3-04 Mev, assuming that the energy lost by the 
protons in traversing the target is several times the total width of the 
level, which is given by Mooring e al. (1954) as 1:5 kev; our targets 
were about 5 kev thick. In this way we obtain a value of 
T,=0-08-6:53 eV. 


§5. Discussion 

The three low lying levels of 25Al are in the order 5/2+, 1/2*, 3/2+ 
which is the order to be expected from excitation of the odd proton on the 
shell model. The resonance level at 3-04 Mev is of odd parity and hence, 
since odd parity levels due to excitation of the d;j. particles would be 
expected to occur at much higher energy, is probably due to excitation 
of the 160 core. This view is confirmed by comparison of the observed 
radiation width with that calculated for single particle transitions by 
Moskowski (1953) which would lead to a value of approximately 50 ev 
in this case. 
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CXIV. Optical Constants of Metals in the I nfra-Red—Conductivity of 
Silver, Copper and Nickel 


By J. R. Brarriz and G. K. T. Conn 
Department of Physics, The University, Sheffield* 


SUMMARY 


The experimental methods used to measure the optical constants of 
metals at infra-red frequencies have already been described. The 
conductivity of evaporated aluminium, silver, copper and nickel is 
reported. Since the penetration depth is small the behaviour is closely 
related to the nature of the surface of the metal; attention is therefore 
paid to the preparation of various representative surfaces including hand 
polished, electrolytically polished and those obtained by evaporating in 
vacuo. 

By giving care to the manner of presenting experimental results, the 
various factors which contribute to the absorption may be discriminated. 
Figures are given of relaxation times, conductivity and the effective 
number of free electrons per atom. Elementary interpretations are given. 


§ 1. INTRODUCTION 


THE depth of penetration of electro-magnetic waves into metals of high 
conductivity is of the order of a micron at 10° cycles per second and 
diminishes as the frequency increases. The behaviour may be described 
using the surface resistance R, defined as the ratio of the electric field 
gradient to the current per unit square at the surface. It can be shown 


that R= 1/(68)= / (270/0c?), lias te etal Mae ele (1) 


where 6 is the skin depth, o is the conductivity and w is the angular 
frequency. At the infra-red frequencies, greater than 10 cycles per 
second, eqn. (1) isinadequate. Since the mean free path of the conducting 
electrons is limited, the surface impedance Z is complex with real and 
imaginary components, R the surface resistance and X the surface 
reactance ; cZ/47 is the ratio of the electric to the magnetic field strengths 
at the surface of the metal (Stratton 1941, Reuter and Sondheimer 
1948). R may be determined from measurements of the ratio, Ro’, of 
the reflected to the incident radiant intensities. To a close approximation 


(en er, eee ee (2) 


The early work of Rubens and Hagen (1904) who measured £,” for many 


metals was interpreted as vindicating eqn. (1) at infra-red frequencies. 
ae 
* Communicated by Professor W. Sucksmith, F.R.S. 
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Doubt has been cast on this view which is in conflict with modern theories 
(Mott and Zener 1934, Dingle 1953). Though R and X are related 
(Robinson 1952) both must be measured if the investigation is limited to a 
particular range of frequencies. Indeed the work here reported shows that 
X is much larger than R so that measurements of R,” alone are of limited 
value. For example fig. 1 (a) shows the surface resistance R and the 
surface reactance X of a film of aluminium plotted against frequency. 
The ‘theoretical points’ are obtained using the dispersion formulae 
discussed below. 

The techniques necessary to examine the characteristics of metals in 
the infra-red are essentially optical and is natural to frame the experiments 
and describe the results in terms of a complex refractive index n—7k. At 
a frequency v the dielectric constant «(v) and the conductivity o(v) are 
then given by (Mott and Jones 1936) 


n?—k?=e(v) and nhyv-e(v)i- «2 a 


Save under certain theoretical conditions discussed by Reuter and 
Sondheimer (1948) the refractive index is equivalent to the reciprocal of the 
surface impedance and 


Z—=RAIX=4nj/{e(n—tk)}: -. yo 


Measurements of the complex refractive index of metals at wavelengths 
greater than 5 » have not hitherto been reported. In two previous papers 
(Beattie and Conn 1955a and Beattie 1955), the difficulties which have 
restricted such investigations were examined, new experimental methods 
were described and preliminary results were presented to illustrate the 
technique adopted. In the present paper observations are reported of 
silver, copper and nickel at various temperatures and over the range of 
2uto12u. Since the behaviour is closely bound up with the nature of the 
surface of the metal, it has been necessary to devote some care and 
attention to the preparation of various representative surfaces including 
hand-polished, electrolytically polished and those formed by evaporation 
in vacuo. 

If means are adopted of distinguishing the various factors, measurements 
of optical constants contribute to studies of relaxation effects, phenomena 
associated with restrictions imposed on the mean free path of conduction 
electrons, surface oxidation and perhaps internal conversion. Attention 
is therefore first given to the manner in which numerical values may best 
be presented to discriminate these phenomena ; in particular to ensure 
that the relaxation and strong absorption of the free electrons do not 
obscure other information. Elementary interpretations are presented. 


§ 2. FRAMEWORK OF INTERPRETATION 


In the simple Drude—Lorentz theory of conduction a relaxation time is 
defined by writing the conductivity in the form 


Coe Ne Tn, ches et ee pee (5) 
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where N is the number of electrons per unit volume, each of charge e and 
massm. The classical dispersion equations are then (Mott and Jones 1936) 


272 
n—k?—e,— — 4 Neer" ? esta ; (6) 
- Meeas eee eee 
477Ne? 
Dipl EE Se ae (7) 
mow (1 + w?7?) Pal 


€, 1s the ‘ static dielectric constant of the atomic cores’. Substitution of 
(5) in (6) and (7) shows that o and 7 suffice to define the dispersion analyti- 
cally. 

This dispersion has been treated in terms of quantum mechanics by a 
semi-classical argument (Krénig 1931) and recently by Wolfe (1954 and 
1955). The model used by Wolfe to calculate the absorption of a metal is 
essentially that of photoelectric transition of conduction electrons to 
higher levels in the same band. Such transitions are forbidden in an 
ideal lattice and Wolfe examined the case in which lattice irregularities 
arise from a random distribution of impurity atoms ; he showed that the 
results of his calculations of the absorption agreed with eqns. (6) and (7) 
over a wide range. The choice of a random distribution of impurity atoms 
to give a model of finite conductivity was dictated by analytical convenience 
and there seems little doubt that formulae of the form (6) and (7) follow 
in more general cases in which for instance the lattice is disturbed by 
thermal agitation. We have therefore used these equations as a basis for 
interpretation of experimental results. 


§ 3. MANNER OF PRESENTATION 


The major contribution to the absorption of high conductivity metals 
‘in the infra-red is made by the conduction electrons and in seeking means 
of discriminating the various factors, the primary question is to distinguish 
this contribution. It is easy to verify that curves of n and & against A 
and v or of n?—k? and 2nk against 2 or v are quite insensitive for this 
purpose. Neglecting ¢,, which is negligible being of the order of unity 
(Mott and Jones 1936), it follows from (5), (6) and (7) that 

e(v)=n?— k? = —4zrog7/(1-+w?7?) and o(v)=nkv=o,|(1-+-w?r?). (8) 
If an Argand diagram is plotted with nkv as ordinate and n?—k? as 
abscissae a straight line passing through the origin is obtained of slope 
—1/(4rr). From such a straight line the relaxation time is easily obtained 
and the value of o, follows from the numerical values of n?—k? or nkv. 

In practice absorption other than that associated with this metallic 
dispersion leads to departure of the experimental points from the straight 
line and to a lateral displacement of the line at long wavelengths. Brief 
attention is given to this elsewhere (Beattie and Conn 1955b). It is 
necessary to take a sufficient number of measurements at long wavelengths 
to ensure that the slope of the line is known. og, is then obtained by 
drawing a parallel line through the origin. Several of these Argand 
diagrams, based on experimental measurements of silver, copper and 
nickel, are presented below. 


8x2 
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Fig. 1 


Argand diagram showing 2nk/A against n?—k? ; experimental results obtained 
with evaporated aluminium. 


Fig. 1 (a) 
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Evaporated Aluminum 
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Surface resistance and surface reactance of evaporated aluminium as a function 
of frequency. 


§ 4. EVAPORATED ALUMINIUM 


The form and use of such Argand diagrams is illustrated in fig. 1 in 
which results obtained with evaporated aluminium are plotted; the 
experimental figures have already been published (Beattie 1955). Circles 
represent the experimental points to which the wavelengths in microns 
are affixed. The deduced straight line representing the contribution of the 
conduction electrons is shown by dashes ; the slope of this line indicates 
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that the relaxation time, +, is 6-1 10-1 seconds and hence that the 
relaxation wavelength, ,, is 11-5. The fitness of eqns. (8) and (9) to 
represent the main features is verified by using A, and 2nk/A and n2—k? at 
9 #, a convenient reference, to calculate values at other wavelengths. 
The crosses indicate the results of such calculations. It is clear that 


absorption other than that of conduction electrons may be present at 
wavelengths less than 8 w. 


Table 1. Characteristics of an Evaporated Aluminium Film 
Thickness of film=1550 +504 u. 
7=6:1X 10-15 sec; A,=11-5 p. 


Ain pw 8-0 9-0 10-0 11-0 12-0 


Gp X 10-1" @.8.u. ae) P13 1-19 1-14 1-23 


Mean value of o)=1-16 x 101” e.s.u. 

d.c. conductivity of film =(1-50 +0-10) x 1027 e.s.u. 
d.c. conductivity of bulk aluminium =3-18 x 1017 e.s.u. 
Effective number of electrons per atom=1-28. 


Table 1 lists various characteristics of this film. The discrepancy 
between oy and the conductivity measured with direct currents is a feature 
of all evaporated films which have been measured ; it will emerge that in 
the case of bulk specimens there is a similar discrepancy. The number of 
electrons per unit volume is obtained from 4, and oo, eqn. (7) ; the number 
per atom, NV, is calculated by assuming that the density of bulk aluminium 
is the same as that of the film the thickness of which was measured by 
multiple-beam, Fizeau fringes (Tolansky 1948). 

It should be emphasized that the construction of an Argand diagram 
from the experimental results serves primarily to clarify the interpretation ; 
the particular diagram chosen presents the relaxation dispersion in its 
simplest form and the nature of the experimental absorption is readily 
grasped. If, as in the present paper, the purpose is restricted to a study 
of the contribution of the ‘ conduction electrons ’, a limited number of 
measurements at long wavelengths is all that is necessary since these 
determine A, and op. 


§ 5. SILVER 


Three specimens were examined. The first preliminary measurements 
were made on a small piece of rolled silver sheet, ground with emery and 
finished with a proprietary polish. The surface was contaminated by 
abrasive from the polish and only one reflection at 84° incidence was used 
so these results, listed in table 2, are of qualitative interest only. 

Two silver mirrors were then prepared by evaporation from an alundum 
crucible in vacuo and condensed on a glass plate at about 40A per second ; 
the final thickness was 5200 A. A study of five reflections showed that 
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the dispersion formula was closely obeyed though there were indications 
of an additional resonance absorption in the neighbourhood of 7 uy. 
Results are listed in table 2. 

Attempts to electro-polish large specimens of annealed sheet were 
only partially successful. Since it is difficult to obtain the critical 
polishing conditions (Gilbertson and Fortner 1942) over a large area, the 
finish was marred by a milky etching. A final polish was given by gentle 
rubbing with chamois leather impregnated with a little magnesium oxide. 
Figure 2 shows the Argand diagram obtained. 


=4900_ 3000. »-2000~ 4000 O 


Figures obtained with annealed polished silver plotted on an Argand diagram. 


Silver does not readily oxidize and the oxide decomposes above 200°c. 
An attempt was therefore made to estimate the temperature coefficient of 
oy by observing the changes at 8 » when the temperature was raised to 
250°c. A section was cut from one of the third, annealed, specimens ; 


Table 2. Optical Characteristics of Silver 


Treatment 7 in sec A,in| o9@.8.u. | o(d.c.), |o(d.c.)/og| WN per 
be €.8.U. atom 
Metal Polish] 2-05 x 10-14 | 38-6 | 12x 1016| 57x 1016 4:8 0-39 
Evaporated | 1-75 x 10-14 | 38-9 | 23-4 x 1016| 33 1016 1-4 0-90 
Annealed 2 i 
and 
polished 1-61 x 10-14 | 30-4 | 19-5 x 1016] 57 x 1026 2-9 0-81 


Temperature coefficient of o) is of the order of 0-001 per °c. 
Temperature coefficient of bulk material is 0-0038 per °c. 
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details of the arrangements for heating specimens are given in the section 
devoted to nickel. The temperature coefficient, quoted in table 2 is much 
smaller than that of the direct current conductivity. 


§ 6. COPPER 


For the study of multiple reflections, pairs of copper specimens were 
prepared in five different ways ; the purity was better than 99-9°%. 
(i) Copper was evaporated from an alundum crucible at a pressure of 
about 5x 10-® mm Hg and condensed on a glass plate at about 164 per 
second. The thickness was 49004. 


Specimens of bulk copper were :— 


(ii) polished with magnesium oxide by hand, see fig. 3 ; 
(ii) buffed, see fig. 3 ; 

(iv) electro-polished without annealing, see fig. 3 ; 

(v) annealed and then electro-polished, see fig. 4. 


Measurements had to be completed within a few hours since copper 
tarnishes readily. Careful and detailed examination is not possible in 
this limited time without using a controlled atmosphere. The deviations 
from a straight line are small in the case of evaporated copper being of the 
same order as those of evaporated aluminium, see fig. 1; in both cases 
these deviations occur at wavelengths less than 8 ». The parameters 
deduced from these specimens are listed in table 3. 


Table 3. Optical Characteristics of Copper 


Treatment T SeC A, ia | oo 68:1. | o(d.¢.) an |o(d.c.)/o9| N- per 
a e.8.U. atom 

Evaporated 

(i) 1-41 x 10-"4 | 26-6 | 1-54 x 101” | 2-50 x 1017 2-0 0-41 
Hand- 
polished (ii) | 1:06 x 10-4 | 20-0 |0-49x 1017} 5-2 10? 10-6 0-21 
Buffed (iii) | 1-80 10-14 | 33-9 | 1-35 1027] 5-2 101? 3°9 0-35 
Annealed 
and electro- 
polished 
(iv) 1-41 x 10-14 | 26-5 | 0-91 x 1017} 5-2 101" 5:7 0-30 


The straight lines drawn in figs. 3 and 4 lead to the values listed in 
table 3. The optical properties are strongly influenced by surface treat- 
ment but the results with the four bulk specimens show common features 
which are believed to be characteristic of copper; departures from the 
straight line predicted by dispersion theory are considered elsewhere. 
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3000 -2500 -2000 “150Q -1000 
Argand diagram using measurements of copper ; upper curve—hand polished : 
lower curve—same specimens buffed. Crosses indicate points obtained 


when the specimen was electro-polished. 


2500 2000 -IS00 -|OOO -500 Oo 
n-k 


Argand diagram using measurements of annealed and electro-polished copper- 
The crosses were calculated using the point at 12 » as reference. 


§ 7. NIcKEL 


The manner of preparing the surface of nickel specimens does not 
appear to have such influence on the measured results as is found with 
copper. Values of p and 4 are less extreme so that a single reflection 
at 84° by one specimen 2 in. x 14in. was adequate. Nickel mirrors do 
not readily tarnish and time is available for detailed study ; they do not 
oxidize at temperatures below 500°c. Measurements were made above 
room temperature using specimens in a pyrophyllite holder fitted with a 
flat nichrome heater and a Pallador thermocouple to measure temperature. 
To ensure no movement of a specimen during heating, the front face 
was firmly located against the rounded ends of three quartz rods. To 
verify that no significant oxidation took place, measurements of the 
optical constants at 24 and 9 were repeated after the specimen had 
cooled again to room temperature. Irreversible changes, when these 
occurred were caused by contamination of the surface with electrolyte 
during polishing ; such cases were disregarded. 
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It is clear from table 4 that the number of electrons per atom, JN, 
remains roughly constant however the surface is prepared. For the 
reasons given in the footnote below, the second evaporated film cannot 
be taken to be representative. Figure 6 shows that when the wavelength is 
less than 5 « the optical constants are largely independent of temperature ; 
the contribution of the conduction electrons to the absorption is small 
and apparently the profile of residuat absorption changes little with 
temperature. 


“1000 —750 »_.* —SOO 250 O 
Argand diagram of hand-polished nickel. 


“1000 750 »2 2 —S00 =o) O 
Argand diagram of annealed, electro-polished nickel at 260°c and 20°c. 


§ 8. DIscussION 


There are few data with which present figures may be compared. 
Ingersoll (1910) extended his measurements of silver, copper and bulk 
nickel to 2-25 while Forsterling and Fréedericksz (1913) examined 
sputtered copper and chemically deposited silver at wavelengths less 
than 5. The figures of Férsterling and Fréedericksz for electro-plated 
copper are not consistent with their other results. The need to maintain 
adequate and detailed control of the treatment of the surface confirms 
that the agreement indicated by typical values listed in table 5 is 
satisfactory. 
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Table 4. Characteristics of Nickel 


11-5 x 1016 


4-6 


1-4 


31 


Fig.| Treatment | Temp. °c T Sec Le Go ©.8.U. 
5 | Hand- 
polished 20°o | 1-14 10-14] 21-5 | 3-72 « 1016 
Buffed 20°c | 0-99 18-7 | 3-80 
yates 5 é Suis 
eb 20°c |0-81 15-2 | 2-68 
polished 270° | 0-70 13-2 | 2:50 
Annealed 20°c | 1-24 23-4 | 3-93 
6 | and | ee |e ee eee | 
electro- 
polished 260°c | 0-84 15-9 | 2-72 
__ | Evaporated 20°c | 1-02 19-2 | 3-16 
(i) 250°o | 0-72 13-6 | 2-36 
Die responce 20°c | 0-50 9-38 | 1-02 
(i) 250°C ‘| 0-40 7-62 | 0-88 


(i) This specimen was evaporated in a vacuum of the order of 5 x 10-® mm Hg from an alw 


crucible at about 504 per second ; the thickness was 30004. 


(ii) This film 30304 was formed by slow evaporation less than 104 per second in a m 
vacuum of about 10-4mm Hg and is thought to have been contaminated by gettering 
residual gas. It was unstable, recrystallized after a few days and peeled from the glass b: 


=I00O0.—-—-750 42 -500.-—- 250 


Argand diagram of nickel evaporated at 504 per second from an alundum 
crucible ; thickness 30004. 
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Table 5. Comparison of Data 
fe a a I SS 


Specimen Authority Ain pw n?2—k? 2nk 
N i (buffed) Present work 2:0 = 63:1 65:8 
Ni (bulk) Ingersoll 2:0 — §9-1 63-1 
Cu (evaporated) Present work 2°25 —154 15:3 
| Cu (electropolished)| Present work 2:25 —12] 22:2 
Cu (bulk) Ingersoll 2:25 —136 24-1 
Cu (sputtered) Forsterling and 
Fréedericksz 2°25 —130 15-5 
Ag (bulk) Present report 2:0 —166 21-7 
Ag (evaporated) Present report 2-0 —230 16-2 
Ag (bulk) Ingersoll 2-0 —187 18-6 
Ag (chemical) Forsterling and 
Fréedericksz 21 —204 28-6 
Cu (evaporated) Present report 5-0 —745 160 
Cu (sputtered) Forsterling and 
Fréedericksz 5-1 —800 169 
Ag (evaporated) Present report 4-0) —710 83-5 
Ag (chemical) Forsterling and 
Fréedericksz 4:0, —821 172 


From the experiments figures of three parameters may be obtained, 
T, oo and N. These are not independent and it is natural to focus 
attention on o, and N since these have an obvious physical significance. 
Caution is required in any general discussion of the present results since 
the body of experimental data is, as yet, not large. 

The effective penetration depth of radiation of wavelength A into an 
absorbing medium is \/47k and with most metals this is of the order of 
10-® cm in the infra-red. It is therefore to be expected that the manner 
of preparing a surface influences o, and the results obtained with silver, 
copper and nickel bear this out. The important factors may be listed 
for convenience as (a) contamination, (b) surface topography, and (c) 
crystal size and strain. To these the influence of “ anomalous skin effect ’ 
must be added (Fuchs 1938, Sondheimer 1952 and Dingle 1953). It is 
not easy to distinguish each since some contribution may be made by 
all but it is possible to offer certain general comments. 

The expressions (6) and (7) leading to (8) which summarize the dispersion 
of ‘free’ or ‘ metallic’ electrons, do appear to describe adequately the 
absorption observed at relatively long wavelengths, A>8p. Figures 1 
and 4 in particular illustrate this. At shorter wavelenghts there are 
significant departures which suggest, particularly in the case of nickel, 
that there are other factors peculiar to nickel which contribute to the 
absorption. 
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The silver surface treated with metal polish yielded a low value of oo 
because of contamination. All low values of o, cannot however be 
ascribed to contamination or surface strain ; even a buffed specimen may 
yield a conductivity as high as that of a specimen which is annealed and 
electro-polished (see tables 3 and 4). oa) may be reduced by surface 
roughness but even with evaporated films in which there can be little 
effect from scratches and relatively little contamination from the conditions 
of condensation, o, is significantly less than the conductivity of the same 
film measured with direct currents or o(d.c.). It does appear that o 
agrees reasonably well in these cases with the best values obtained with 
bulk material. 

In general o, is less than o(d.c.) and with the metals reported here 
o(d.c.)/o, is of the order of 2-5 or more; the lowest value of this ratio 
is furnished by silver (table 2). If measurements of nickel and silver are 
any guide the temperature coefficient of o, is very small compared with 
that using direct currents. A simple calculation shows that it is not 
necessary to assume that the additional contributions to the resistance 
have any dependence on temperature. These might be attributed to 
‘imperfections ’, a word used in its widest sense, but it is necessary to 
extend measurements to low temperatures. It is difficult to resist the 
conclusion that this discrepancy is fundamental and it is natural to seek 
the cause in current theories of ‘ anomalous skin effects ’ since conductivity 
at infra-red frequencies is essentially a surface phenomenon. The 
magnitude of 5/1, / being the electronic mean free path, certainly falls in 
the appropriate range but direct comparison cannot be made with 
microwave measurements at low temperatures since the present experi- 
ments lie in the region of metallic dispersion. Current theories of 
anomalous skin effects predict a different frequency dependence from 
that of classical dispersion. This is not in agreement with experiment. 
More extensive measurements are required before the demands to be 
made of theory can be clearly defined. 

The values of NV are surprisingly low compared with those derived from 
the Hall effect (Wilson 1953). This cannot be due entirely to the use of 
the free electronic mass in the calculation. The figure for nickel, 0-15 per 
atom, presumably refers to the s-electrons since the contribution of the 
d-electrons to relaxation is negligible. With copper and aluminium, 
N is low being 0-35 per atom and 1-3 per atom respectively. The figure 
for silver, 0-9 per atom is more in line with the valency. It may be that 
simple relaxation theory has been interpreted too closely though it may 
be noted that with liquid metals (Kent 1919), N is to be identified with 
the chemical valency to a surprising, and convincing, degree. Measure- 
ments of liquid metals made in this laboratory confirm this to a material 
extent. Whatever the significance of the numerical values, with silver 
and nickel N does not vary over the range of temperature studied. 
If a polished surface of silver or nickel is etched, there is little effect on 
7 but XN is significantly reduced. 
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CXV. Resonance Absorption of Nickel in the Infra-Red Region 
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[Received May 26, 1955] 


SUMMARY 


Means are described whereby the absorption in the infra-red of ‘ free ” 
or ‘ conduction’ electrons may be abstracted from the total which is 
measured experimentally. When this is done, it is clear that there is 
residual absorption in nickel which may be caused by electron transitions 
from the d- to the s-band. The mechanism is being further explored. 


THE conductivity of silver, copper and nickel in the infra-red has been 
discussed in a previous paper (Beattie and Conn 1955) where it is shown 
that the contribution of the conduction electrons, that is those which 
obey simple dispersion theory (Mott and Jones 1936, Wolfe 1954), gives a 
straight line through the origin when 2nk/X is plotted against n?—k? on an 
Argand diagram. n—vik is the complex refractive index and dA is the 
wavelength. With many specimens the experimental points show depar- 
tures from the simple linear relation particularly at shorter wavelengths, 
A<8 » and at least in some cases these must be considered to be significant. 
To assess these additional contributions it is first necessary to dissociate 
from the measured values the absorption arising from the metallic 
dispersion. Caution is required in the interpretation because of the 
influence of the surface examined. It is difficult to avoid the conclusion 
that with the nickel surfaces, however prepared, there is a real absorption 
arising from transitions of d-electrons to vacancies in the s-band. Such 
absorption has been predicted by Mott (1935). The present purpose is to 
present the experimental evidence for selective or ‘ resonance ’ absorption 
by nickel ; confirmation of the interpretation is being sought by examining 
the behaviour of copper-nickel alloys. 

Expressions representing metallic dispersion (Mott and Jones 1936, 
Beattie and Conn 1955) are commonly described in terms of a relaxation 
time 7 with an associated relaxation frequency, v,, and wavelength, d,. 
These are simplified by defining a reduced frequency f so that 


fHv/v,=A,/A=oT=2r7v and v,=1/(277)=Cc/),. 
27v=w is the angular frequency and \ is the wavelength. Then 
n—k?—n,2= —Amroy/(1-+f?) and 2nk=4mro,/{f(1+f2)}. (1) 


* Communicated by Professor W. Sucksmith, F.R.S. 
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No" is negligible being of the order of unity so that the straight line through 
the origin is 


— (n?—k?) = 2nkf Tee gs ee) 


and may be used to determine 7 and o, (Beattie and Conn, loc. cit.). In 
practice it is not difficult to draw the straight line of eqn. (2) through 


Fig. 1 


SPECIMEN - BUFFED NICKEL 


é Ingersoll 


© Present work 
x Classical dispersal 


-I500 “090 2, 0)®) O 


<—§(nrk) 


Argand diagram for a buffed nickel surface with 2nk/A plotted against n® —k?. 
The straight line is taken to represent the contribution of * metallic 
dispersion’; the crosses represent the computed position of the experi- 
mental wavelengths associated with the circles on the experimental 


curve. 


the origin and parallel to the run of points at long wavelengths ; measure- 
ment must be so extended that the linearity is confirmed and the slope 
determined. Values of 2nk/A obtained with buffed nickel are plotted 
against n2—k? in fig. 1; the present results, measured at the wavelengths 
indicated, have been extended to shorter wavelengths by using the data of 
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Ingersoll (1910). The position of points on the straight line which 
correspond to each of the measured wavelengths can be calculated 
knowing X,.; these are represented by crosses. The additional absorption 
at each wavelength is given by the lines or vectors joing each cross to 
the corresponding circle and the components of each vector may be used to 
construct ‘extracted’ curves of 2nk and n?—k? against v or A. The 
graphs so obtained from fig. 1 are shown in fig. 2. If the additional 
absorption is caused by resonance there is no contribution to n?— k? at the 
resonant frequency and the corresponding vector is then parallel to the 
2nk/A axis. Contours of the residual absorption bands so constructed 


Fig. 2 


iO a : 
O tO 20 30 hy in ev. 
Curves showing, as a function of frequency, residual values of 2nk and n?—k? 
after allowance has been made for metallic dispersion. The specimen 
was buffed nickel ; see fig. 1 


depend to some little extent on the line drawn to represent relaxation 
of the conduction electrons but since this must pass through the origin the 
choice is limited. 

The only a priori assumption is that relaxation effects can be represented 
by such a line. The absorption curve is then extracted by calculation 
using A, and the experimental values of n?—k? and 2nk. The errors are 
therefore those associated with measurement of n2—k? and 2nk and with 
the fitting of the straight line which is very rarely less precise. Since the 
contribution of metallic dispersion becomes progressively larger at longer 
wavelengths, small departures from the straight line imply larger changes 
of 2nk. The percentage error associated with ‘extracted’ figures of 
2nk and n?—k? increases therefore ; the experimental difficulty lies in 
measurement of the phase difference 4 between the components of the 
reflected radiation and caution is necessary in interpreting such absorption 
curves. 
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Figure 3 shows dispersion curves obtained with two specimens prepared 
in entirely different ways. For fig. 3 (a) a film of evaporated nickel was 
prepared under good vacuum conditions, the pressure being less than 
5x 10~° mm Hg and condensed on glass at about 50 A per second; the 
thickness was 3000 A. With unstable specimens prepared under poor 
conditions, about 10-* mm Hg, similar absorption bands emerged but the 
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Residual values of 2nk and n?—k? plotted against frequency. The specimen was 
a film of evaporated nickel condensed on glass to a thickness of 3000 A. 


(0) 


Residual values of 2nk and n?—k? plotted against frequency. The specimen was 
of hand-polished nickel. (It should be noted that the scale of abscissae 
in figs. 3 and 4 is quite different from that used in fig. 2.) 


maximum absorption moves to shorter wavelengths and is less than half 
the value shown in fig. 3 (a). Figure 3 (6) illustrates the behaviour of 
hand-polished nickel. In each case, figs. 2, 3 (a) and 3 (6), there is a 
doublet structure near the maximum with a regular fall of 2nk as the 
frequency is increased. 

When the effect of temperature is studied the conductivity is changed 
but the temperature coefficient is much smaller than that of bulk material 
(Beattie and Conn, loc. cit.). Figure 4 was obtained using an annealed, 
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electro-polished specimen. The broad characteristics of the previous 
figures are repeated and, most important, the effect of temperature on the 
band structure is small. At the lowest frequencies it is significant, as is to 
be expected, but at higher frequencies the influence is negligible. 

The present brief communication has two purposes. We wish to report 
and illustrate the practical possibility of distinguishing the contribution of 
the free or metallic electrons to the measured values of 2nk and n?—k?. 


Residual values of 2nk plotted against frequency. The annealed, electro- 
polished nickel was examined at 20°c and 260°c. 


Attention may be then directed to residual absorption if this is present. 
Such absorption may arise from many causes, some associated with the 
nature of the surface, some characteristic of the material. In the case of 
nickel it is difficult to avoid the conclusion that there is evidence of 
resonance absorption. It is natural to ascribe this to transitions between 
the d- and s-bands ; at the moment there is little profit in detailed quanti- 
tative examination of the data. We have preferred to seek confirmation 
by extending the range of experimental study to alloys of copper and 
nickel. This work will be reported in the near future. 
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SUMMARY 


The aim of this paper is to give a simplified calculation of the energy 
required to form a vacancy in a metal, which is applied to the noble and 
alkali metals. Although all such calculations are subject to considerable 
error unless many refinements are introduced, it is felt that such a 
calculation will give a useful indication of the various factors which 
determine the energy and the variation from metal to metal. 

An outline of the method is as follows. The metal is represented as 
a large spherical box in which the positive charge of the ions is uniformly 
distributed and the electrons are free to move. When an ion is removed 
from the centre of the sphere and spread over its surface, the energy 
of the free electrons changes by 4H,,, because the electron waves undergo 
phase shifts in order to screen. the vacancy and the volume of the metal 
box changes. 4H,, is proportional to the Fermi energy of the metal H,, 
and is estimated to be about {#,. For the noble metals another relevant 
contribution to the energy needed to create a vacancy is that of the 
closed-shell repulsion between ions: for Cu this is about —0-:3 ev. 

The energy to move a vacancy is also discussed ; the electronic con- 
tribution is fairly small and in noble metals the term due to the closed-shell 
repulsion predominates. This is probably somewhat underestimated if 
one uses repulsive potentials valid in the neighbourhood of the equilibrium 
interatomic distance. 

The general agreement between the theoretical results and the experi- 
mental values for the energies to form and to move a thermal defect in 
the alkali and noble metals is satisfactory and suggests that vacancies 
are responsible for the transport of matter also in the alkalis. 


§ 1 
In recent years there has been considerable interest in the problem of 
point defects in metals owing to the role that these defects play in the 
transport of matter within solids. Evidence bearing on these defects 


* A preliminary report on this work was given at the Xe Conseil de Physique 
Solvay (Bruxelles, 1954). 

+ Permanent address: Istituto di Scienze Fisiche, Universita di Milano, 
ltaly. 
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has been obtained from widely different experiments, including the study 
of diffusion, radiation damage and the mechanical properties of metals. 
In contrast with the fairly large amount of experimental work, the 
amount of quantitative theoretical work on the problem has been rather 
limited. 

The only detailed calculations of the energies required to form and 
to move point defects in metals are those of Huntington and Seitz (1942, 
1949) and Huntington (1942, 1953) for Cu; rough estimates have been 
made by Bartlett and Dienes (1953) for vacancy pairs in Cu and by 
Paneth (1950) for linearly relaxed interstitials (crowdions) in the alkali 
metals.* The results of Huntington and Seitz are not in a form which 
can easily be generalized to other metals. This has recently led Nachtrieb 
and Handler (1954) to attempt to draw conclusions about the mechanism 
of transport of matter in metals from an empirical relation between the 
activation energy for self-diffusion and the latent heat of melting. The 
underlying assumption of their work is that the vacancy mechanism of 
self-diffusion should have an activation energy simply related to the 
latent heat of sublimation ; we shall see that this is not the case. 

The purpose of this paper has been to devise a simple model which 
would allow some physical insight into the processes of creation and 
motion of vacancies in monovalent metals, and which would make it 
possible to obtain information of a general nature about the energies 
involved. 


§ 2 

To calculate the energy needed to create a vacancy in a monovalent 
metal, we represent the metal as a large spherical box in which the 
positive charge of the ions is uniformly distributed and the electrons are 
free to move. We then compute the change in energy that the electrons 
undergo when we remove an ion from the centre of the sphere and spread 
the charge over its surface. For the noble metals we will have to add 
to the change in the electronic energy computed in this way the change 
in the energy due to the strong non-coulombic forces between the ions. 

With this model the unperturbed wave functions of the free electrons 
in the perfect metal are of the form r—! sin (kyr—4lm) (l=0, 1, 2, ...), 
while the perturbed wave functions in the metal with the vacancy will 
have the asymptotic form r~! sin (k,/r-+-y,—4lz), where 7, is the phase 
shift of the / waves (Mott and Massey 1950 a) due to the vacancy. The 
condition that the perturbed and unperturbed wave functions alike 


* In his estimate Paneth (1950) neglects for instance that the electron gas is 
pushed out of the volume occupied by the interstitial ion and thus his value 
for the energy to form a crowdion is bound to be too small by an ev or two 
(C. Herring, private communication). This contribution to the energy to 
form an interstitial in a metal has been neglected also by Huntington and 
Seitz (1942) and by Huntington (1953). 
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should vanish at the boundary of the spherical box (r=R) yields the 
equations 


k,R=n7-+ila (taal oe): gee em LY 
k, R=nr-+4lr—n, (eee es Ve ee (2) 

It follows that ee, 
R(k,'—k,)+n,=0. SO Can es on se, 


The change in energy of an / electron owing to the change Jk, of its 
momentum is given by #?k,4k,/m. Therefore the energy change of the 
free electrons upon removal of an ion from the centre of the metal sphere 


is equal to H2len,(It) R 
mR 


AE ., ea! D+) dk. 


l 
where ky is the momentum of the electrons at the Fermi level and R/z 
is the lena of states for / electrons. The factor 2 originates from the 
spin. 
In the Born approximation (Mott and Massey 1950 b) 


n(k)=— aan V(r) {J psja(kr) Pr dr Sea a oh C8 


where V(r) is the perturbing potential, that is the potential field of the 
vacancy screened by the conduction electrons. Using the addition 
theorem of Gegenbauer (Watson 1922), 


2kr 
220+ 1) {Sp 41/2(h7) P= TE 
l 
we can rewrite (4) in the following form 
ky? R 
AB y'=s8y| Virjanr@dr, 2. 8) 
3772 J 9 


This is the result given by first-order perturbation theory ; the energy 
change undergone by the electrons is given by the integral of the per- 
turbing potential multiplied by the unperturbed electron density k,3/37?. 

The use of Friedel’s theorem (Friedel 1952) for the phase shifts at the 
Fermi level, still using Born’s approximation, allows us to put 4H,,’ in 
a form independent of V(r). Indeed, following Friedel, we may write 


Dy 
— 3'(21+-1)y, (kp) =—1, pene eI, ees (7) 
EN) 


since a vacancy in a monovalent metal can be considered as an impurity 
atom with an excess valence equal to —1. The left-hand side is equal, 


in Born’s approximation, to —(kym/h?x?) le Vir)4rr? dr. Thus 


oe 
AE / F 


el = “3m = 2H, : . . . : . . (8) 


where E, is the Fermi energy. To improve on this result which is based 
on the Born approximation, we shall have to return to eqn. (4) (see § 3 
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AE,,' is not the entire energy change undergone by the electrons 
upon creation of a vacancy. We still have to compute the energy change 
due to the spreading of the ion we have extracted over the surface of 
the metal sphere ; this causes an expansion of the spherical box in which 
the electrons are enclosed. The energy of the N electrons in the box 
is (3/5)NHy, and the Fermi energy Hy, is proportional to V-?* where V is 
the volume of the box. Thus the change in energy caused by an expansion 
OV of the box is SV 


3 
NEy = 


t 


9 2 
AR,’=—5 2B). ins ee 


OU 


where the number of atoms in the metal has been taken to be equal to 
the number N of free electrons (monovalent metal). The increase in 
surface energy connected with the expansion of the box is negligible 
being of the order of N~1/3. 

Thus for the total change 4H, of the energy of the electrons 

AB j= AE g-+ SE "=~ By. oh rn 

This result neglects the changes in exchange and Coulomb correlation 
energies of the free electrons. However, as Friedel (1954) has noted, 
the result is not altered if one includes in the perturbing potential V(r) 
the change in average exchange energy. Since the volume over which 
the electrons are distributed is essentially the same in the initial and 
final states, it is perhaps permissible to neglect also the change in 
Coulomb correlation. 

$3 

Formula (8), and thus formula (10), have been obtained using the Born 
approximation for the phase shifts », of the free electrons of the metal 
under the action of the repulsive field of the screened vacancy. This 
will give too large a value for 4H,,’, and thus for 4#,,, because the 
Born approximation overestimates appreciably the phase shifts for 
repulsive fields. It is clearly of interest to see how much the value 
of JEH,,’ is affected. 

To calculate 4H,,’ with exact phase shifts we have to choose a definite 
form for the repulsive potential of the screened vacancy. We have 
adopted a constant spherical barrier (Jongenburger 1953, Abelés 1953) 
for which the exact phase shifts can be expressed analytically. To 
estimate the error in 4H,,’ caused by the Born approximation it was 
felt to be sufficient to use exact phase shifts for the s waves which undergo 
considerably larger shifts than the waves of higher /. The calculation 
can be done at once for all the monovalent metals, for which 


key = (97/4) = 1-919 


(r,=radius of the atomic sphere). 
The potential we use is of the form 


Vi(r)=D for r<a, \ a1) 
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We take a equal to r,, and we choose D (=/2k,2/2m) using Friedel’s 


condition for the phase shifts at the Fermi level. This condition reads 
now 


2 
“(21+ 1)inilkr) Jom = { eee) +s) }pomt = = {ng P) fexact= il (12) 
> j 2 ¢ 2 ss a(7,k Vi 
where 72+ fnlkp) }pom=— 5 aoe eee 
6 k a rk)? . 
(rolls) nom 3G pgp (Si Py C08 ryly—reip) 6. (14 
and (Mott and Massey 1950 a) 
aioli) tora = tan" Jee eae / {(rgkp)?— (7k) —r,k 
O\'’P) fexact v/{ (r,kp)2—(r reo)? } V Mn 
For k<k,y 
pal 
{(n0(*) fexact= tan | EE tanh V/{(rsko)?— a fg. (16) 


A#,,’ is then obtained by adding the exact energy change for the s waves 


to the energy change for the higher / waves computed in the Born 
approximation 


og ee? 1 rE sin a COS & 
AE = mre E (7ko)*(Tskp)>— ie (rho)? (vate = | a 3. i) 


0 «x 
Z 3 rk, 
+= {Cee a I. a tan! 
x 
Xx (a= tanh V(trdko)*—0*}) dx 


rokp b v : 
— lee xz tan-1 (oe tan V(02—(rbo)}) dz} | = sa(llieg) 


To obtain 4H, we have again to add to 4E,,’ the energy 4H,,"=—?Ey. 
The value of r,k, which fits Friedel’s condition (12) is found to be 
approximately 1:64. The resulting value of 4H,, is 4H, within the 
accuracy of the numerical calculation. 

The use of exact phase shifts also for J>0 would still yield a 4E,, 
proportional to H, since the exact phase shifts 7,(k,) depend only on the 
product 7,k,; it is difficult to estimate a priori the magnitude of the 
change in the proportionality factor but it appears unlikely that it is 
sizeable. Thus it does not seem warranted to perform the corresponding 
numerical calculations owing to the approximations involved in the model 
itself. 

§ 4 

In the case of the noble metals, the energy to create a vacancy cannot 
be taken to be equal to the energy change undergone by the conduction 
electrons. Important cohesive forces in these metals are the ion—ion 
closed-shell repulsions (Mott and Jones 1936) and the Coulomb correlations 
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between conduction electrons and 3d-core electrons (H. Brooks and 
J. Friedel, private communications) ; while the former contribute only a 
small fraction of an ev to the cohesive energy, the latter are thought to 
contribute more than | ev. 

The Coulomb correlations between conduction and core-electrons 
operate inside individual atomic polyhedra and thus should not contribute 
in a relevant way to the energy needed to transfer an ion from the interior 
of a metal to its surface. 

In the approximation of a rigid lattice, the closed-shell repulsions 
between ions (central forces) would contribute to the energy to create a 
vacancy an amount equal to their contribution to the cohesive energy. 
However, the metal lattice is not rigid and upon removal of a metal ion 
from the interior the nearest neighbours relax to new positions. The 
relaxation occurs in such a way as to minimize the energy and thus the 
repulsive energy decreases more than in a rigid lattice. 

To estimate the contribution of the repulsive forces to the work to 
create a vacancy, allowing for the relaxation of the neighbouring atoms 
only, we note with Huntington and Seitz (1942) that the creation of a 
vacancy in a face-centred metal causes the rupture of 6 bonds and the 
alteration of 84; there are 24 bonds between the nearest neighbours of 
the vacancy, and only 60 of the 84 bonds between nearest and next- 
nearest neighbours are altered to a first approximation. The expression 
for the change in repulsive energy is* 


AE p=12R{d(1+A)}+48R {d(1+A+d2)"/2}-+4 24R{d(1—A)}—9OR(d) (18) 


where A(r) is the repulsive energy per ion pair at distance 7, d is the 
equilibrium interatomic distance and Ad is the inward displacement of the 
12 nearest neighbours of the vacancy. If we use for R(r) the Born—Mayer 
exponential form A exp (—7/p), AE, becomes 


\ 


4HEp=A exp (- *) 
P 


Ir d d 
x {12 exp (- =) +48 exp (- : 5) +24 exp (* a) —90}. (19) 


To maximize this decrease of repulsive energy, we minimize the A-depen- 
dent part. 

We have performed the numerical calculations for copper using the 
two Born—Mayer potentials recently proposed by Huntington (1953) and 
the Born—Mayer potential originally used by Huntington and Seitz (1942) 


* This expression assumes that each of the 12 nearest neighbours of the 
vacancy undergoes an equal displacement toward the site of the vacancy. 
Obviously this is a very special type of relaxation ; angular rearrangements 
may in fact be important (N. Cabrera, private communication). It seems safe 
to say, however, that the order of magnitude of the repulsive energy gain will 
not be affected. 
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in the correct form reported by Huntington (1953).* We found that 
AH, is about —7-5R(d) (with d=2-55 4) rather than —6R(2-55 A) as 
it would be in the rigid lattice approximation : \ ranges between 0-015 
and 0-020. The numerical values of 4H, for the three potentials used 


are respectively { —0-40 ev, —0-28 ev and —0-24 ev. A reasonable value 
to take is —0-3 ev. 


For silver and gold no formulae have yet been proposed for the closed- 
shell repulsion but the order of magnitude of the gain in the repulsive 
energy upon creation of a vacancy is certainly the same as for copper. 


§5 

The calculation of the energy needed to move one of the nearest 
neighbours of a vacancy in a monovalent metal from its equilibrium 
position halfway toward the vacancy, thus reaching the saddle-point 
configuration, is a difficult problem. We will attempt only some semi- 
quantitative considerations on the important contributions to this 
energy. 

We note first with Huntington and Seitz (1942) and Huntington (1942) 
that the saddle-point configuration, in which a metal ion is halfway 
between two neighbouring vacancies, can be envisaged as a perfect 
lattice with two neighbouring half-vacancies in it, if one assumes that 


* It is difficult to say whether any of these potentials is to be preferred. If 
so, it should be the potential adopted by Huntington and Seitz (1942) in the 
form reported by Huntington (1953, footnote 14) since its two empirical 
parameters were determined from the three elastic constants of copper taking 
proper account of the electrostatic contributions to the shear constants com- 
puted by Fuchs (1936). The first repulsive potential proposed by Huntington 
(1953) (H,) R(r)=0-053 exp {—13-9(r—d)/d} ev per ion pair, is based instead 
on the assumption (Zener 1950) that the actual electrostatic contributions to 
the shear constants are appreciably smaller than those calculated by Fuchs 
(1936), owing to a redistribution of the electron gas during the shear distortion. 
The fact that the values of the shear constants computed by Fuchs (1936) for 
the alkali metals are in very good agreement with experiment (Mott 1952) appears 
to indicate that Zener’s assumption is unrealistic, and Dr. Huntington (private 
communication) points out that the calculation that he gave to substantiate it 
(Huntington 1953, appendix A) is incorrect. The second Born—Mayer potential 
proposed by Huntington (1953) (H,) R(7)=0-038 exp {—17-2(r—d)/d} ev per 
ion pair, was simply chosen as an opposite extreme to the first, so as to be 
fairly close to the potential adopted by Huntington and Seitz (1942) (H.S.), 
which can be written as R(r)—0-032 exp {—17(r—d)/d} ev per ion pair. We 
have not considered the repulsive potential used by Zener (1950), 0-078 exp 
{—(19/4/2)(r—d)/d} ev per ion pair, which is derived on the radical assumption 
of a vanishing electrostatic contribution to the shear constants, and without 
using the experimental value of the bulk modulus. It should be mentioned 
that the repulsive potential attributed by Zener (1950, p. 350) to Huntington 
and Seitz (1942) is affected by the misprints which occurred in the original 

er. 
“} The last number is given incorrectly by Huntington and Seitz (1942), 
p. 320, as —0-4 ev. 
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the moving ion deforms in such a fashion as to fill only spaces which do 
not belong to the atomic polyhedra of its neighbouring ions. 

The electronic contribution to the energy to create two half-vacancies 
in a monovalent metal can be estimated along the lines of §2. To a 
first approximation this energy does not differ from the electronic energy 
change (4/15)H, which accompanies the creation of a vacancy ; indeed the 
‘box expansion’ term —3H#y is replaced by two equal terms each 
amounting to —}H, and similarly the ‘ phase-shift ’ term }Hy is replaced 
by two equal terms each amounting to 4H#,. Roughly we should thus 
expect no electronic contribution to the energy of movement for a 
vacancy in a monovalent metal. However, corrections to this first-order 
estimate should be made, mainly owing to the non-spherical symmetry 
of the holes in the electronic distribution corresponding to the two half- 
vacancies (Huntington 1942, section II. 6). These corrections could 
easily give in the closed-packed metals a contribution of a few tenths of 
an ev. 

For the noble metals one must consider also the contribution of the 
closed-shell repulsion to the energy of movement of a vacancy. This can 
easily be estimated following Huntington and Seitz (1942). The saddle- 
point configuration can be obtained from the perfect lattice by removing 
two ions from their adjacent lattice sites and placing one of them at the 
surface and the other at the centre of the rectangle formed by the four 
ions which are nearest neighbours of both vacancies ; upon the intro- 
duction of an ion at the centre of this rectangle, each of the four ions at 
the vertices displaces itself outwards along the line joining it to the 
added ion, by an amount equal to its distance from the added ion 
multiplied by A. The change in repulsive energy can be considered to 
be composed of three parts due: (i) to the rupture of 17 bonds (12 to 
create the first vacancy, 11 to create the second and 6 recovered from 
the atom placed at the surface), (ii) to the introduction of an ion at the 
centre of the rectangle formed by the four common nearest neighbours 
of the two vacancies, allowing only for the relaxation of the ions at the 
vertices of the rectangle, and (iii) to the relaxation of the remaining 
14 neighbours of the two vacancies. The contribution (ii) is given by* 


4[R(fay/10)+-R{fav/6(1-+A)}+ ER{ha(2+-4d)¥2} 4 2B {ha(2-+d)"2} 
t 2K {ta(2—A)/2\+ Rfta(2—2r)/2} 
+ 2K {$a(2—3A)/2}—7-5R(da4/2)). . . . . . . (20) 
For (iii) we take with Huntington and Seitz (1942) a ratio 14/12 of the 
decrease in repulsive energy caused by the relaxation around a single 


* This formula is given incorrectly by Huntington and Seitz (1942), eqn. (6) ; 
instead of 
28 {a(2+4A)1?}+ 2K (fa(2—A)*9}-+ BR {ha(2—2d)2/2} 
they write 
ER {$a(2+2A)1/23+3R {4a(2—A)2/2}. 
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vacancy which we considered in § 4. The contribution of the closed-shell 
repulsion to the energy of movement of a vacancy is obtained by sub- 
tracting from the sum of (i), (ii) and (iii) the repulsive energy gain which 
accompanies the creation of a single vacancy (see § 4). 

We have performed the numerical calculations for copper as we did 
in § 4 and we have again used the two Born—Mayer potentials proposed 
by Huntington (1953) (H, and H,) and the Born—Mayer potential 
adopted by Huntington and Seitz (H.S.) (1942) as reported by Huntington 
(1953). The values of (i), (ii) and (iii) for the three potentials are given 
below (in ev) : 


(i) (il) (iii) 


H, 0-01 1-12 (4=0-047) 0-09 
H, —0-65 1-06 (A=0-052) - —0-07 
H.S.* —0-55 0:90 (A=0-052)  —0-06 


The final result turns out to be fairly insensitive to the particular form 
of potential used, ranging from 0-52 ev for the H, potential to 0-62 ev 
for the H, potential. 

This value for the contribution of the closed-shell repulsion to the energy 
of movement of a vacancy in copper may be somewhat too small. The 
Born—Mayer potentials that we have used to calculate it are really valid 
only in the neighbourhood of the equilibrium interatomic distance. 
A rather approximate analysis of the (27°c) pressure-volume isotherm 
for copper up to pressures of 0-5 megabar (Walsh and Christian 1955), 
following the lines of Zener (1950), suggests that the repulsive energy per 
ion pair at distances appreciably smaller than the equilibrium distance 
(such as the distance r=a,/6/4 which appears in the repulsive energy 
contribution to the energy of the vacancy saddle point) may be larger 
than that estimated from the equilibrium forms of the repulsive potential 
by a factor of a few units. This implies that the contribution (ii) may 
be underestimated, while the same will not be true for the contributions 
(i) and (iii) in which one deals respectively with the equilibrium distance 
and with distances differing only slightly from it. It is tempting to 
stress the analogy with the case of NaCl where the values that Mott and 
Littleton (1938) computed for the energies of movement of Na* and Cl~ 
vacancies in NaCl (0-5 ev for Nat vacancies ; 0-6 ev for Cl” vacancies) are 
well below the experimental values (Etzel and Maurer 1950, Aschner 1954, 
Seitz 1954) (0-85 ev for Nat vacancies, >1 ev for Cl- vacancies) while the 
values they computed for the energies of formation are in fair agreement 
with experiment. It is natural to attribute the inaccuracy of the 
theoretical energies of movement to the repulsive energy as Mott and 
Littleton did, since the remaining contributions should be computed 
with considerable accuracy in an ionic solid. 


* These numbers are incorrectly given by Huntington and Seitz (1942), 
Table II as —0-69, 0-51 and —0-19 respectively. 
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The experimental values for the energies to create and to move point 
defects in monovalent metals are at present uncertain, particularly for 
the case of the noble metals, and it should be stressed that theoretical 
results cannot be a decisive element in these disputes owing to the 
approximations which are made to perform the calculations. 

Accurate values for the activation energy for self-diffusion are available 
for lithium (Holcomb and Norberg, private communication) and for 
sodium (Nachtrieb, Catalano and Weil 1952), while an approximate value 
is available for rubidium (Holcomb and Norberg, private communication). 
Values have been reported also for the energies needed to create thermal 
defects in the alkali metals (Carpenter 1953, MacDonald 1953, 1955) 
and the reasonable assumption has been made that the defects 
created are the defects responsible for self-diffusion ; in this way, the 
activation energies for self-diffusion have been separated in an energy of 
formation and an energy of movement. The available results are 
summarized in table 1. The more probable values appear to be those 


Table 1. Experimental Values for the Activation Energy for Self- 
diffusion (U) and for the Energy to form (U,) and to move (U,) 
a Thermal Defect in the Alkali Metals (in ev) 


' 
Li Na K Rb Cs 


U 0-57 0-45 — 0-41 - 
U,(MacDonald 1953) 0-40 0-39 0-39 — — 
(MacDonald 1955) —_ _ 0:39 ae — 
(Carpenter 1953) “= 0-26 0-24 — = 
U,(MacDeonald 1953) 0-17 0-06 — —- wae 
(Carpenter 1953) — 0-19 — a — 


which MacDonald has obtained by measurements of electrical resistance 
(MacDonald 1953) and of specific heat (MacDonald 1955) at high 
temperatures. Thus the thermal defects in the alkali metals appear to 
have a very small energy of movement. This checks well with our 
discussion of the energy of movement of a vacancy in the alkali metals 
(see §5) which gave a zero value. The theoretical values for the energy 
of formation of a vacancy in the alkali metals that we computed in 
§§ 2 and 3 are shown in table 2; for Li, the only alkali metal for which 
the effective mass m* of the conduction electrons is known to be 
appreciably different from the free electron mass m, we have introduced 
in Hy the correcting factor m/m*=0-69 (Brooks 1954). The values 42, 
agree well with the more probable experimental values for the energy 
of formation of a thermal defect, even if they are still slightly too large. 
It seems reasonable to conclude that the thermal defects responsible 
for matter transport in the alkali metals are vacancies. 
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For the noble metals, accurate values for the activation energy for 
self-diffusion are available (Kuper, Letaw, Slifkin, Sonder and Tomizuka 
1954—Copper ; Slifkin, Lazarus and Tomizuka 1952—Silver ; Gatos and 
Kurtz 1954—Gold) and vacancies are normally accepted as the defect 
responsible for the transport of matter, owing to the theoretical results 
of Huntington and Seitz (1942, 1949) and the experimental observation 
of the Kirkendall effect (Seitz 1950, Zener 1950). However, alternatives 
have been proposed for the way in which the total activation energies 
should be divided into the energies of formation and of movement of 
a vacancy (Manintveld 1952, 1953, Meechan and Eggleston 1954, Kauff- 
man and Koehler 1955). These experimental values are given in table 3. 
On the ground of the available experimental evidence the more probable 
values appear to be those which give an energy of movement for the 
vacancy larger than the energy of formation (Seeger 1955). The theoretical 
values for the energy of formation of a vacancy are reported in table 4 
and favour such an assignment. The value we estimated in § 5 for the 
energy of motion of a vacancy in Cu (c. 0-6 ev) appears then to be too 
small by a factor 2; as we noted in § 5, this is not unreasonable. 


Table 2. Theoretical Values for the Energy (U,) to form a Vacancy 
in the Alkali Metals (in ev) 


Li Na K Rb Cs 
(m/m* =0-69) 

Hgt 3-30 320 4 = 2188. 8 1°58 

U Bon —(4/15)Ey (§2) 0-88 0-85 0-57 0-49 0:42 

U,=(1/6)Ep (§ 3) 0-55 0-53 036 031 0-26 


For copper, our theoretical results are in good agreement with the 
results of Huntington and Seitz (1942) and Huntington (1942). Indeed 
Huntington and Seitz (1942) obtain an energy of formation for a vacancy 
ranging from 1 to 1-5 ev, and the more detailed calculations by Huntington 
(1942) effectively give an energy of about lev. Huntington’s calculations 
(1942) give also an upper limit of about 1-8 ev for the energy of movement 
of a vacancy. It is frequently stated in the literature that the results 
of Huntington (1942) give an energy of formation definitely larger than 
the energy of movement and that the latter is computed to be about | ev: 
these statements, however, should be corrected. The high value (1-8 ev) 
for the energy of formation that Huntington (1942) reports as the result 
of his calculations is faulty since it includes twice the electrostatic inter- 
action of the shielded vacancy with its nearest neighbours : the interaction 
is included already in the calculation leading to the value of 1-2 ev 


+ The values of the Fermi energy Hy for the alkali metals quoted in the 
literature differ somewhat (cf. Kittel, C., 1953, Introduction to Solid State 
Physics (Wiley), p. 240; and Wilson A. H., 1953, Theory of Metals (Cambridge), 
p- 270). The values reported in this table have been computed anew with the 
equation rk» =(97/4)1/3 adopting the following values for the lattice parameter : 
Li, 3-484; Na, 4-254; K, 5-204; Rb, 5-624; and Cs, 6:05 A. 
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(Huntington 1942, p. 332), and should not be added again as Huntington 
(1942) does in Table II of his paper. The actual value for the energy of 
formation which follows from Huntington’s calculation is about | ev, 
since one must subtract from the quoted value of 1-2 ev the gain in 
repulsive energy which accompanies the breakage of six bonds : according 
to Huntington (1942), there is no relaxation around a vacancy in copper 
because the net electrostatic field of the shielded vacancy at its nearest 
neighbours is repulsive, a result which is likely due to the neglect of the 
phase shifts for / greater than 2. For the energy of movement, it should 
be stressed that the value of about 1 ev which can be obtained from 
Huntington’s calculation (1942) by subtracting the energy to form a 
vacancy (1 ev) from the estimate of the energy of the saddle point 
(1-6-+-0-4) ev* is a lower limit, since a relevant contribution to the energy 
of the saddle point has been neglected as Huntington (1942, section II. 6) 
emphasizes. 


Table 3. Experimental Values for the Activation Energy for Self- 
diffusion (U) and for the Energy to form (U,) and to move (U,) 
a Vacancy in the Noble Metals (in ev) 


F Cu Ag Au 

U 2-03 1-96 1-96 
U, (Meechan and Eggleston 1954) 0-90 — 0-67 
(Seeger 1955) 0-90 c. 0-8 0-67 
(Manintweld 1952, 1953) 1-15 1-27 1-27 
(Kauffman and Koehler 1955) — _ 1-28 

U, (Meechan and Eggleston 1953) 1-13 — 1-29 
(Seeger 1955) 1-13 c. 1-2 1-29 
(Manintweld 1952, 1953) 0-88 0-69 0-69 
(Kauffman and Koehler 1955) — 0-68 


Table 4. Theoretical Values for the Energy (U,) to form a Vacancy 
in the Noble Metals (in ev) 


Cu Ag Au 
Ky 7:04 5-51 5-54 
AE 480m —(4/15)Hyp (§ 2) 1-88 1:47 1-48 
AE ,=(1/6)Ey (§ 3) 1-17 0-92 0-92 
AEp (§ 4) —0:3 c. —0°3 c. —0°3 
U,=4E + AER 0-9 c. 06 c. 0-6 


* The value 1-6 ey stands for the value 1:57 ev which actually follows from 
Huntington’s calculation (1942, section IT. 5) (through a misprint 1-57 is given 
as 1:55). The value 0-4 ev is an estimate of the repulsive energy contribution 
to the energy of the saddle point which is obtained by considering only the 
contributions (i) and (ii) of §5: according to Huntington (1942) there is no 
relaxation around a vacancy and thus (iii) is zero. The values of the nearest 
neighbour contribution quoted by Huntington (1942, Table III) suffer from the 
same fault discussed in the text just above. 
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A final comment should be made on the widespread opinion that the 
energy to form a vacancy in a metal would be related to the heat of 
sublimation (see, for instance, Nachtrieb and Handler 1954). This is, of 
course, true for a rigid lattice bound by central forces but there is no 
reason to expect it to be true for a metal, where the main contribution 
to the energy to form a vacancy comes from the change in energy of 
the conduction electrons ; as we have seen, this change in energy is 
related to the Fermi energy of the metal and not to the heat of sublimation. 
In this connection it is interesting to consider the case of the noble metals 
where the latent heat of sublimation and the latent heat of melting have 
different trends ; while the heat of sublimation of copper is smaller than 
that of gold, the opposite is true for the heat of melting (Slater 1939). 
The energy of formation of a vacancy, related to Hy, follows the trend 
of the heat of melting, and not that of the heat of sublimation; of 
course the correspondence between the energy of formation of a vacancy 
and the heat of melting is purely qualitative. 
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CXVII. CORRESPONDENCE 


Effect of Heat Treatment on the Internal Friction 
of Aluminium Crystals 


By M. Levy 
The Haloid Company 
and M. MrtzcEr 
University of Illinois 


Tue strengthening of metal single crystals by increasing the rate of 
cooling after an elevated temperature anneal has recently been observed 
in simple shear tests of zine crystals (Li, Washburn and Parker 1953) 
and. in tensile tests on aluminium crystals over a wide range of purity 
(Maddin and Cottrell 1955). Maddin and Cottrell proposed that the 
strengthening is due to the restraint on the motion of dislocations resulting 
from the condensation of vacancies on the dislocation lines to form jogs 
which move with difficulty and increase the applied stress necessary for 
the nucleation of slip. The purpose of this note is to cite additional 
evidence in support of this proposal based on measurements of the 
damping of aluminium crystals in forced longitudinal oscillation in the 
kilocycle range. The strain amplitude dependent damping observed in 
these studies, which results from the periodic motion of dislocations 
(Read 1940), is related to the earliest stage of plastic flow before disloca- 
tion multiplication has begun. A low damping corresponds in a general 
way to a high yield stress. In this sense, the results of these studies are 
in good agreement with the results of Maddin and Cottrell. 

The measurements were made using the composite oscillator method 
described by Read (1940). Principal attention was devoted to maximum 
strain amplitudes from 10-® to something more than 10~° (the upper 
part of the obtainable range) since the effects of heat treatment were 
more pronounced in this range. The specimens were 4mm diameter 
single crystal rods about 8 cm long of 99-992 or 99-996°% aluminium 
grown from the melt. The decrement (the fraction of the vibrational 
energy lost per half cycle) was measured at room temperature after 
various heat treatments carried out with the specimen supported in an 
alundum boat. It was not necessary to remove the oxide film since the 
decrement was found not to be affected by its presence. 

A marked effect on the decrement of the rate of cooling after an anneal 
had been observed in preliminary work by J. Intrater and 8S. Shapiro 
(unpublished), and this was investigated in some detail for a one day 
anneal at 640°c. Although in each case there was some decrease in the 
decrement with time at room temperature, the decrements characteristic 
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of each cooling rate remained in the same relation to one another. In 
fig. 1 are shown measurements made after one to three days at room 
temperature when the changes with time had become very slow. These 
data are typical of those obtained for several crystals as regards both the 
effect of cooling rate and the magnitude of the decrement. Both the 
magnitude of the decrement and the rate of its dependence on strain 
amplitude are seen to be greatest for furnace cooling, intermediate for 
air cooling and least for water quenching. The broken curve in fig. 1, 
which represents an earlier water-quenching test, is given to provide an 
indication of the reproducibility of the data. The reproducibility was 
always satisfactory, even after a number of intervening annealing 
treatments, so that the structure of a crystal after annealing at 640°C 
was independent of thermal history so far as these measurements were 


concerned. 
Fig. 1 


Decrement x 105 


Maximum strain amplitude x 106 


Decrement at 34 ke of a single crystal of 99-992°% aluminium as affected by 
the rate of cooling after annealing 1 day at 640°c. FC=furnace cooled 
at 50°c hr-} (test No. 6); AC=cooled in still air (test No. 5); WQ= 
water quenched (test No. 7). The broken curve represents an earlier 
water quenching test (test No. 3). 


The condition of a specimen after furnace cooling may be taken for 
the moment as reference since a slow cool after annealing has been the 
traditional practice in damping studies. Now, the effect of increasing 
cooling rate to depress the damping is not a trivial one due to the strains 
produced by rapid cooling since it has been observed that the damping 
of a specimen is lower when it is furnace cooled at 50°c per hour than 
when it is furnace cooled at a much slower rate. Nonetheless, a word 
of justification is required with regard to the use of a quenching treatment 
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in this work since, even though it was performed with care, this is a 
drastic treatment where measurements of a highly structure-sensitive 
property are concerned. It is to be noted that the thermal and mechanical 
shock to which the crystal is subjected during the quench do not appear 
responsible for the very low decrement observed after quenching since, 
in the vicinity of room temperature, both these factors markedly increase 
the decrement of these crystals. It is believed that the effect of these 
factors is obscured by the opposite and much stronger effect which 
depresses the damping, and thus the low decrement after quenching is 
believed to result from the same phenomenon which produces a decrease 
in the decrement with increased rate of cooling in the range of low 
cooling rates. 

The measurements made in the range where amplitude dependence is 
substantial show hysteresis, i.e., the decrement remeasured at a given 
strain amplitude immediately after an intervening exposure to a higher 
strain amplitude is now higher than the first reading. This is attributed 
to the breakaway of the dislocation line at some of the sites where it 
was pinned and is consistent with the conditions treated by Koehler 
(1952). He considered the dependence of the decrement on the average 
free length of dislocation line and showed the decrement to vary inversely 
as a power of the linear density of pinned sites. Then, since the properties 

of the aluminium crystals studied in this investigation did not change 
- non-reversibly in the course of successive anneals, the general dislocation 
configuration is considered to have remained unaltered and the effects 
of heat treatment are considered to have resulted from reversible changes 
in the density of pinned sites. 

The observed effects of cooling rate cannot be ascribed to the redis- 
tribution of sites pinned by impurity atoms since the Cottrell or Suzuki 
atmospheres of the dislocations would be dispersed at 640°c and would 
be most efficiently reformed to give the lowest decrement at the slowest 
cooling rate, contrary to observation. It is believed, in agreement with 

the proposal of Maddin and Cottrell, that the major part of the effects 
observed can be attributed to the condensation on dislocations of lattice 
vacancies which are unable to reach the crystal surface during the time 
the crystal is cooling to room temperature after the anneal. The 
argument, briefly stated, is that a vacancy has an associated dilatational 
stress and would be attracted to the edge part of an extended dislocation, 
and the jog produced would move with difficulty thus serving to pin the 
dislocation at low stresses ; the stable vacancy pairs formed after quench- 
ing (Bartlett and Dienes 1953) have associated shear stresses and could 
interact with the screw part with similar results. 

The number of vacancies available was estimated by using Nowick’s 
(1951) estimate of the volume self-diffusion coefficient of aluminium as 
0-45 [exp (—1-43/k7')] em? sec? and by taking 65% of the activation 
energy as the energy required for vacancy formation in accordance with 
the results of Kauffman and Koehler (1955) for gold. The energy to 
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form a vacancy in aluminium is thus estimated as 0-9 ev and the fraction 
of vacant sites at 640°c as 10-5. A network of intersecting dislocations 
is postulated having a density of the order of 108 lines per cm? corres- 
ponding to a value of 10-4 cm as the possible free length of the dislocation 
lines, as required from consideration of the plastic properties of soft 
single crystals (Cottrell 1953). Thus, the ratio of the length of dislocation 
line per unit volume to the number of vacancies per unit volume at 
640°C is estimated to be of the order of 10-®cm per vacancy. The 
number of vacancies originally present is therefore so large that, even if 


Effect of Final Annealing Temperature before quenching on the Decrement 
at 35 ke of a 99-996% Aluminium Crystal. In each case the Crystal 
was initially annealed one day at 640°C 


Decrement x 10° at a maximum 
strain amplitude of 5 x 10-° 


Test Subsequent heat treatment SS 
No. 8-11 min. | 40-42 min.| 1-3 days 
after after after 
quench quench quench 
1 | quenched directly 5:0 3°8 37 
2 | cooled in 5 days to 550°c—quenched 13 9-2 
3 | cooled in 6 days to 450°c—quenched 19 13-3 
4 | quenched directly 4-6 3°6 3°8 
5 | cooled in 10 days to 350°c—quenched 32 21 iN 


many escape to the crystal surface during the quench or combine with 
one another to form larger groups, there will still be a sufficient number 
migrating to the dislocations and producing jogs to reduce the average 
free length of the dislocation lines to a value much below 10-4cm and 
thus to account for the low decrement observed after quenching. 

If, instead of quenching from 640°c, the crystal were quenched from a 
lower final annealing temperature where the equilibrium vacancy concen- 
tration is smaller, the decrement would be expected to be higher. This | 
was observed to be the case, as indicated by the data given in the table, 
both soon after quenching and after one to three days when the decrement 
was no longer changing at an observable rate. For each test in this 
series, the crystal was first annealed at 640°c, then it was cooled very 
slowly in the furnace in order to permit the excess vacancies to migrate 
to the crystal surface and the equilibrium concentration of vacancies 
to be approximately maintained, and finally it was quenched from the 
desired temperature. The magnitude of the decrement after quenching 
is believed to have been determined by the relative extent to which damp- 
ing was decreased due to vacancy condensation and increased as a result 
of thermal and mechanical shock. The large decrement observed after 
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quenching from 350°c suggests that the latter effect was then predominant. 
The estimated concentration of vacancies at 350°c (which is 1/200 of the 
concentration at 640°c) still appears high enough for the vacancies to 
have been able to exert some influence. 

A further test of the vacancy hypothesis is afforded by observations 
of the rate of decrease in the decrement with time at room temperature 
after quenching. Maddin and Cottrell had noted that the critical shear 
stress measured soon after quenching was lower than that obtained after 
the crystal had remained at room temperature for a few days. They 
attributed this to the gradual migration to the dislocations of those 
vacancies which were still dispersed immediately after quenching. In the 
present work, the decrement after quenching fell rapidly with time, as 
indicated in the table by the values obtained in three of the tests 40-42 
minutes after quenching. The major part of the change was accom- 
plished within a few hours in all crystals quenched from elevated 
temperatures. Now, the diffusion coefficient for vacancies can be 
obtained by dividing the volume self-diffusion coefficient by the equili- 
brium fraction of vacancies (Bardeen and Herring 1952). The present 
estimates of these quantities lead to an estimate of the vacancy-diffusion 
coefficient at room temperature as greater than 10-19 cm? sec~!. This is 
large enough for diffusion over distances of the order of 10-*cm in a 
few hours to be nearly complete, i.e., the diffusion parameter a/2D1/2t1/2 
is about 0:1, and is thus consistent with the view that the decrease in the 
decrement with time is due to the diffusion of vacancies. This effect 
cannot be due to the diffusion of impurity atoms to the dislocations ; the 
volume chemical diffusion coefficients for solute atoms in aluminium 
range from 10-2 to 10-24 at room temperature as calculated from 
Nowick’s (1951) estimates, and even if the diffusion of the solute atoms 
were enhanced by the excess vacancies present after quenching, the 
solute atoms could not have effective diffusion coefficients even within 
several orders of magnitude of 107 1°. 


This work was done at Columbia University and was supported in 
part by the Office of Ordnance Research. 
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Impurity Diffusion in Polar Crystals 


By J. TELTOW 
Institut fiir Kristallphysik, Deutsche Akademie der Wissenschaften, Berlin 


[Received June 3, 1955] 


In a previous paper submitted to this journal, A. B. Lidiard* has given 
a very lucid and elegant statistic analysis of the diffusion of divalent 
ions in monovalent ionic crystals, e.g. Cd++ in AgBr. The same problem 
was treated by Schéne, Stasiw and the author a few years ago, starting 
from general thermodynamic transport equations. The purpose of this 
note is to show that, with a slight simplification, these calculations still 
remain valid and, in Lidiard’s special case, lead to the same result. 

Lidiard’s principal objection concerns our use of a diffusion coefficient 
De for unassociated Cd++ ions, of effective charge +1, as distinct from 
the diffusion coefficient D,, for associated pairs of cadmium ions and silver 
ion vacancies. Owing to their neutrality, the diffusion of the latter is 
not influenced by the Nernst diffusion potential, in contrast to the 
behaviour of the unassociated Cd++. Assuming that the divalent 
impurity ions can move only by jumps into adjacent vacancies, this kind 
of double-line description of an essentially uniform diffusion mechanism 
seems somewhat unsatisfactory. 

The Dg term was introduced by us to account for the ideal case of 
“no association’ when only accidental meetings of cadmium ions and 
silver ion vacancies occur. Diffusion of the cadmium will then still be 
possible. Meanwhile it has become clear that instead of putting the 
number per cm? of pairs n=0, this ideal case is described correctly by the 
vanishing of the free energy of association 4G. According to the associa- 
tion equilibrium (Stasiw and Teltow 1947, Teltow 1949) 


Nin 
(V,—n)N, 


this is equivalent to K,=12 or n=12 N,N,/(N,+12N,). Here V,, 
N, and NV, are the numbers of cation sites, totally added Cd++ ions and 
silver ion vacancies per cm, respectively, and 12 is the number of nearest 
cation neighbours of a cation site in a face centred cubic (NaCl type) 
lattice. For small degrees of disorder, 12V,<N, and n=12N,N,/N, 
which agrees with direct statistic computation. Thus the D,, term alone 


gives rise to a non-vanishing macroscopic diffusion coefficient D in this 
case too. 


=K,=l2.exp (AGP); > ee 


* The author is very much indebted to Dr. Lidiard for receiving a preprint 
of this paper. 
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This ideal case should not be confounded with the extreme case ‘ no 
pairs at all’, characterized by 4G—+— oo (strong repulsive forces between 
the partners), and therefore K,—0, n=0 and D=0. 

We therefore agree with this part of Lidiard’s criticism and shall 
obtain a correct description of the AgBr—Cd case by putting Dg=0 in our 
calculations (see below). This procedure automatically rules out the 
influence of the diffusion potential. It must be emphasized, however, 
that the general equations with Dg+0 are still valuable in cases where 
diffusion of the impurity ions proceeds by two different mechanisms. 
This has been demonstrated by Schéne, Stasiw and the author for the 
diffusion of cuprous ions in AgBr, which partially move via interstitial 
sites. In such cases, the diffusion potential has an important bearing on 
the migration of the impurity ions. It will even, to a small amount, 
modify the probabilities of the orientation jumps of the pairs in the 
AgBr-—Cd case and thus involve the diffusion coefficient. This effect, 
though presumably negligible, should be examined in a refined analysis. 

From eqn. (12) of our paper (Schone, Stasiw and Teltow 1951), then, 
we obtain by putting Dp=0, in Lidiard’s notation 


2nN , 

a, 2N,N,—(N,—n)2 SOS gece Whe eg te ers (2) 

where D,=D,, by a reasoning analogous to eqns. (2), (3) and (5) of our 
paper, is given by 

Ov Main Nita |e eee eis 28 1 | AD) 


Here a is the nearest anion—cation separation (half lattice constant), 
and w, is the number of jumps of the cadmium ion into an adjacent 
vacancy per sec. If w, is not small compared with the jump probability 
w, of the vacancy to another attached cation position, then in (3) ws, 
has to be replaced by w,w,/(w,+w,.). This has been stated already by 
Schottky (1950). 

Equation (2) includes the more general case of the number of thermal 
defects being comparable with the number of impurity ions (e.g. silver 
halides above 300°C with additions of about 0-1%). To obtain Lidiard’s 
result, we only have to specialize by neglecting the thermal disorder, 
ie. V,=N,—n. Evaluation of n from (1) and insertion in (2) then leads 
immediately to 

De DN ea(leea KIN IN Jat) 2 ae oe Sey. (4) 


which is Lidiard’s final equation. If, on the other hand, the thermal 
defects have to be accounted for, the calculation of V,, 2 and N,, (number 
of positive lattice defects, e.g. interstitial cations (Frenkel type) or anion 
vacancies (Schottky type)) by means of the eqn. (1), 


N~Np=Ni-” (neutrality condition), . .91i-.. 2, «+ (5) 
and N NN p=A1 (law of mass action for thermal disorder) (6) 
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yields, upon insertion in (2), no simple relation. Nevertheless, this 
procedure is sometimes indispensable. As an example we take the 
approximate proportionality of D with N, resulting from (4) if 
4K,N,/N,<1. It follows from (1), (2), (5) and (6) that this is no more 
valid when V,/N, becomes as small as K,~1/?. For still smaller additions 
N, the diffusion coefficient is, expanding to the linear term 


D=D,(1—L+K,I4N JN ,+ ...); 


where L=(1+K,K,~/?)-1. It should be kept in mind that the whole 
association model is only an approximation for very small (thermal or 
foreign ion induced) disorder of the lattice. Escaping towards the range 
of higher impurity concentrations in order to simplify the equations may 
be sometimes of doubtful practical value. 
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Strain-Hardening of a Metal by Laminar or Turbulent Flow 


By W. A. Woop 
Metallurgy Research Department, University of Melbourne 


[Received June 3, 1955] | 


THE theories of strain-hardening due to G. I. Taylor and Mott suggest 
that the more a metallic structure is disordered the more is it likely to 
develop barriers to free dislocations and the more rapidly therefore is it 
likely to strain-harden. Cottrell (1952) has generalized this view by 
distinguishing between a turbulent flow that disorders the structure 
during deformation and a laminar flow which does not ; turbulent flow 
hardens more effectively. 

It is therefore of interest to test by experiment how effective different 
forms of disordering may be in practice. One form particularly amenable 
to experiment is produced in polyerystalline metal by uni-directional 
deformation. This is the disorientation of the grains that gives rise to 
X-ray asterisms. 

Tests on this particular disordering, summarized in the present note, 
show rather surprisingly, however, that it has little effect. All it affects 
is the stability of the hardened structure. 


W. A. WOOD Phil. Mag. Ser. 7, Vol. 46, Pl. 19. 


- -reflection (400) rings from copper. (1) As annealed ; (2) after 5° 
Saal are aoe 5 Vv PN nl) ae cyclic stress, hardness 67 VEN c 
(4) after heating stretched specimen, softening by reerystallization ; 

(5) a specimen hardened cyclically to 63 VPN ; (6) same specimen after 
heating, softened to 44 VPN without recrystallization (same spots as 


before heating). 


i) 
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These points may be demonstrated by x-ray back-reflection photo- 
graphs from polycrystalline copper. 

Figure | (figs. 1-6, Plate 19) is typical of the structure in the undisordered 
state after annealing and before deformation. The perfection of the grains 
is shown by the sharpness of the reflection spots. The hardness in this 
state was 35 VPN. 

Figure 2, from the same material after 5°/ extension under a tensile 
stress of 7 tons/m?, is typical of disordering by uni-directional deformation. 
The grains reflect over an extended angular range, causing the reflections 
to diffuse into the continuous rings. Each grain is therefore reduced to 
differently oriented domains. This flow, evidently turbulent, raised the 
hardness to 66 PNV. 

Figure 3 is from the same material hardened to the same degree, 
67 VPN, but by 10° cycles of alternating tension—-compression at 
+7 tons/in.? applied in a Haigh fatigue-testing machine. This is typical 
of the alternating deformation about a mean zero stress that is commonly 
used in fatigue testing. Clearly this hardening occurs without significant 
disorientation of the grains, or in fact any detectable disordering, for the 
reflection spots are virtually as sharp as they were initially. Thus we have 
the same hardening by flow that relative to the foregoing is laminar. 

The comparison, then, demonstrates the first point that disordering by 
disorientation is not essential to strain-hardening. 

Its influence on stability was examined by heating specimens such as 
the above in stages of one hour at 200°, 250°, 300°C... and subjecting 
them after each stage to hardness and x-ray tests. The disordered and 
undisordered states then responded differently as follows. 

Disordered state. Here (i) the disorientation is only removed by 
recrystallization ; (ii) the strain-hardening also is removed only by 
recrystallization ; (iii) this recrystallization sets in and is completed at 
a relatively low temperature. 

Thus fig. 4 from the disordered specimen after heating at 350°C is 
typical of the first signs of softening. New reflection spots from re- 
crystallized grains appear in the still continuous rings. The spots 
multiply at the expense of the rings until recrystallization is complete, 
in this case at 400°c. 

Undisordered state. Here (i) the strain-hardening may disappear 
largely before recrystallization ; (ii) recrystallization is postponed to a 
higher temperature. 

These features may be illustrated by a fatigue specimen first hardened 
to 683 VPN at -L7 tons/in.2.. This specimen was so held that the same 
areas could be x-rayed before and after heat-treatment and the same 
reflection spots thereby kept under observation. In this way a sensitive 
and non-destructive test for recrystallization is obtained, for new grains 
and therefore new reflection spots are produced as soon as recrystallization 
sets in. Figure 5 shows the reflection spots in the initially hardened 
state of 63 VPN. Figure 6 shows that the same spots persist after 
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heating to 600°c, when the hardness fell to 44 VPN. Hence this softening 
has occurred before recrystallization. In this instance recrystallization 
did not occur until the temperature was raised to the relatively high value 
of 650°C. 

These further tests show, then, that the disordered state is the less 
stable, since it recrystallizes first. 

It is believed that these observations are of interest in showing that 
hardening of the same pure metal can occur by different dislocation- 
arrays according to the mode of deformation. 

Thus, on the one hand, the array responsible for strain-hardening under 
the turbulent uni-directional deformation might be interpreted on 
Mott’s theory as a relatively coarse array of piled-up groups, each 
containing many unit dislocations ; for the main slip bands produced by 
such deformation are spaced widely, and the slip jumps are avalanches of 
many unit slip movements. Further it would appear that these piled-up 
groups are located on the network of sub-boundaries delineating the 
disoriented domains; for both disorientation and strain-hardening, 
removal of which implies a dispersing of the pile-ups, disappear together 
on heating. 

On the other hand the array produced by the alternating deformation 
is more likely to be a finer distribution of smaller pile-ups, on the following 
grounds. The criterion for the comparatively laminar flow in the 
fatigue specimens is that the stress cycle should initially be built up by 
alternating increments that are small, preferably producing not more 
than about --0-1°% strain, so that the stress never produces any large 
uni-directional strain (Wood and Davies 1953). A strain of 0-1°% in 
a grain of say 10-? cm means an extension of 10-* cm. Thus the criterion 
is that the plastic strain should be restricted to slip jumps of less than some 
10-8 cm, that is, to what may be classed as ‘fine’ slip. In this case, 
therefore, the pile-ups must consist of fewer unit dislocations; and 
since they are capable of hardening as effectively as the coarser pile-ups, 
they must be more closely spaced. Whether they also are located on 
sub-boundaries it is not possible on the evidence to say, since such fine 
slip does not significantly disorient the grain. 

On this view the disoriented structure becomes more unstable because 
of the higher stress concentrations associated with the larger pile-ups. 
Recrystallization occurs before they can disperse in any other way. 

It is suggested that experiments of this kind, dealing with pure metal 
hardened to the same degree but by different internal changes, may help 
in a much-needed physical interpretation of the dislocation-arrays 
postulated by the theories of strain-hardening. 
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Recent beta-spectrometer studies of the 4°Sc decay (Keister and Schmidt 
1954) give an indication of the order of emission of the two gamma-rays ; 
the decay scheme is shown in fig. 1. The spins and parities have been 
determined by angular correlation and polarization experiments (Brady 
and Deutsch 1948, 1950). Both gamma-rays are expected to be electric 
quadrupole. Weisskopf’s (1951) theoretical relation for the lifetimes of 
gamma-ray transitions yield lifetimes of 110-1! sec and 7x 10-™ sec 
for the 1-12 Mev and 0-89 Mev transitions respectively. The results of 
delayed coincidence experiments (Murdoch and Webb 1954, Whittle and 
Porter 1953) place an upper limit of 1-5 x 10~9 sec on the lifetimes of the 


tates in 4°Ti. 
states in 4°Ti Fig. 1 


oo Tia 


O0.36MEV 99.9% 


Decay scheme of *Sc. 


Experiments have been carried out in this laboratory to examine more 
closely the lifetimes of these states in *°Ti. A fast coincidence apparatus 
similar to that described by Bell, Graham and Petch (1952) with a resolving 
time 7)2-5x10-* sec was used. In this case we have used separate 
side channels with a display on a cathode ray tube (see Ferguson and 
Lewis 1953, Azuma and Lewis 1955) to effect the necessary pulse height 
selection. Three coincidence experiments were performed: (1) a y-y 
coincidence experiment with **Sc using ®Co to provide prompt coinci- 
dences for comparison, (2) a f-y, y-y ‘ self-comparison ’ experiment 
involving only the **Sce source, (3) a f-y coincidence experiment in 
which the 46Sc decay was compared with that of Co. 

The first experiment determines the lifetime of the intermediate state of 
46Ti by ‘ comparing ’ the y~y coincidence resolution curve with that of 
609. The two counters used plastic phosphors (from Pilot Chemical 
Inc.) in the form of right circular cylinders, 1 in. diameter by 1 in. long, 
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surrounded by magnesium oxide and contained in aluminium cans thick 
enough to stop all beta particles, and mounted on E.M.I. photomultipliers 
Type 6262. The side channel kicksorters accepted energy bands 100 kev 
wide with one side accepting only pulses due to the 1-12 Mev gamma-ray 
from 4°Ti, the other side accepting pulses from both gamma-rays. The 
positions of these bands were set to include regions of similar pulse height 
distribution from each source, and also to exclude coincidences due to 
Compton back-scattered quanta. The Co lifetime has been found to be 
less than 10-11 sec (Bay et al. 1955) ; but even supposing the lifetime to 
be as large as 5X 10-" sec this source provides a ‘ prompt > comparison 
since either counter can accept the first gamma-ray, and the two y-rays are 
detected with almost equal efficiencies. For these reasons a slight broad- 
ening of the resolution curve rather than a shift of the centroid is produced. 
A series of coincidence resolution curves were obtained, the curves for 
each source being taken alternatively, with one run consisting of a group of 
four curves. The shifts of the centroids give directly a measurement of the 
lifetime (Bay 1950). The result of three runs gave a lifetime of less than 
5x 10-11 sec, which is regarded as the limit of stability of the apparatus. 

Since the lifetime of the intermediate state of *°Ti is less than 
5x 10-11 see, a ‘ self-comparison ’ method was used to determine an upper 
limit for the lifetime of the second excited state. The plastic phosphor on 
one counter was replaced by a 1 cm cube of trans-stilbene forthe detection 
of both beta and gamma rays. The output from the other counter, the 
gamma-counter, was passed through a discriminator rather than a kick- 
sorter, and was biased to accept pulses from both gamma-rays. Thus a 
f-y coincidence resolution curve modified by about 20° y—-y prompt 
coincidences was obtained. Since the gamma-counter accepts pulses 
due to both gamma-rays, the B-y resolution curve contains contributions 
from the lifetimes of both excited states of *°Ti. The y-y prompt 
curve was obtained by inserting a 270 mg/cm? aluminium absorber 
between the source and the f-counter. The side channel kicksorter 
‘ windows ’ were again set in positions where the shapes of the pulse height 
distributions were similar for each source. A lead absorber, 0-5 em 
thick, was placed between the source and the gamma-counter to attenuate 
back-scattered quanta. Several runs, again in groups of four, were 
obtained and showed that the excited state lifetimes of 4°Ti are both less 
than 5 x 10-1" sec. 

To confirm the two above results, a third measurement was made by 
comparing the 6~y cascade of **Se with that of Co, which is known to 
have lifetimes of less than 10~™ sec (Bay 1955). A thin stilbene crystal, 
1-5 mm thick, was used on the beta-counter, and to reduce further the 
effects of back-scattering a 0-5 cm lead absorber was placed between the 
source and the gamma-counter. A comparison of the coincidence 
resolution curves for the two sources shows that the lifetimes of the two 
excited states of 4°Ti are both less than 5x 10-™ see, thus confirming the 
above results. Figure 2 shows a set of coincidence resolution curves for 
this experiment. 
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Observed coincidence counts 


ts =2 fe} +2 +4 +6 


Delay x 10-° sec. 
B-y coincidence resolution curves for *6Sc and ®Co. 
O 2*S¢ —— + — — Co 
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Tue lifetime of the 280 kev state of 2°TI, produced following the P-decay 
of 23Hg, was shown by Deutsch and Wright (1950) to be less than 
310-9 sec. Conversion coefficient data suggested the transition was an 
M1, E2 mixture (cf., for instance, Wilson and Curran 1951, Rose et al. 
1951). Lifetimes of 2x 10- sec and 5x 10~® see are expected from the 
Weisskopf (1951) formula for pure M1 and pure E2 transitions respectively. 
The resonance scattering experiments of Metzger and Todd (1954), 
Metzger (1954), involving ?Tl, indicated a lifetime (7'y) of the state 
against y-emission, of between (10+4) x 10-1 sec and (9+4) x 10~1° sec. 
This heightens the interest of a direct determination which was in progress 
here and which we report in this communication. Comparison by 
results of the two methods has hitherto only been possible for 1°8Au and 
then only in a limited manner, the region of overlap being in general a 
difficult one for both methods. 

A coincidence unit of resolving time ~2 x 10~® see was employed, using 
two selected E.M.I. photomultipliers type 6262 (ef. fig. 1). The fast 
mixing unit with delay cable was similar to that of Bell, Graham and 
Petch (1952), but the two sets of pulses, contributing to the coincidences 
finally observed, were selected to lie in narrow pulse height bands by 
kicksorting after mixing (cf. Ferguson and Lewis 1953). At a fast 
coincidence, the cathode ray tube beam was brightened, and lengthened 
pulses from the dynodes of the two counters were applied to the X and Y 
deflecting plates respectively causing brightened spots to appear. Spots 
appearing on a small selected area of the display activated an auxiliary 
photomultiplier, operating a scaler. This apparatus allows comparisons 
to be made of the coincidence resolution curves involving radiations of 
different kinds. 

A delayed coincidence curve was obtained for the B-y rays of Hg. 
It was compared with a prompt curve for coincidences between the two 
simultaneous quanta from the annihilation of positrons, generated by a 
22Na source enclosed in aluminium. A stilbene crystal, 1 em cube, was 
used on one counter, A plastic scintillator, 1 in. long, 1 in. diam. (from 
Pilot Chemical Inc.) surrounded by a 4 in. thickness of magnesium oxide 
in a 35 in. aluminium container, on the other counter, detected y-radiations 
only. The source lay between the two scintillators. The narrow energy 
channels were positioned to eliminate coincidences involving Compton 
back scattering or x-rays, at settings where the respective pulse height 
distributions of the radiations of ?°°Hg and *Na were similar in shape. 
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C, photomultiplier collector ; B, dynode; L, limiter; A, amplifier; d, main 
delay cable; h, stub; F, fast crystal mixer (~210-% sec); D, discri- 
minator; T, slow coincidence unit (~1 sec); P, pulse former; e, 1, sec 
delay ; G, gate and lengthener; XY, cathode ray tube; M, auxiliary 
multiplier ; 8, scaler. 

Fig. 2 
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Coincidence counts as a function of delay. 
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A lead sheet 7 in. thick with a hole } in. diam. was inserted between the 
source and the stilbene as an additional precaution against Compton 
back scattering. Following pile irradiation of mercury the radioactive 
mercury had been allowed to decay ubefore se to remove isomers of 1°’Hg. 
For a prompt coincidence run in which ?*Na was substituted for 7°Hg 
the counters were each moved 5 ecm in line from the ?2Na source, so that the 
1-28 Mev y-ray contributed: only very slightly to the coincidence count ; 
the small delay times associated with it (cf. Ferguson and Lewis 1953) 
were then of no account. Several pairs of runs were taken in succession. 
Figure 2 shows results for one such set of 4 runs. The observed centroid 
shifts, giving the lifetime, lay near 3 x 10~1° sec. 

As a check on these results, comparison of the delayed coincidence 
curve for the B-y rays of 2°Hg was made with the delayed coincidence 
curve for the B-y rays of 18Au. 1%Au yields a state of half-life 
(1-0+1-7) x 10-1 sec (Graham and Bell 1951), which can be considered 
prompt for this experiment. A stilbene crystal, 1:5mm thick, was 
employed on the f-side. The narrow energy channels were used as 
previously also the lead screen. The results for many pairs of runs in 
succession showed similar centroid shifts to those found above. 

The lifetime of the excited state of 2°°T] is, on the basis of all the 
measurements here (3-4-1) x 10-1 sec. 

Because of internal conversion the value of the lifetime (7'y) against 
y-emission, of the 280 kev state, determined by the resonance experiments 
quoted earlier, will exceed by 20°, the corresponding figure for the 
lifetime of the excited state. The lifetime of the excited state from 
the resonance method (8+3-2, 7-2+3-2) x 10-7” sec, is therefore in a range 
above that found directly here, the limits of error just adjoin. In 
conclusion, comparison with the Weisskopf formula shows a substantial 
electric quadrupole component in the transition. 


One of us (R. E. A.) is indebted to the Nuffield Foundation for a 
Studentship. 
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